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Ml.  AUTPUCT 


-ae  resu3.ts  of  a  three-year  research  program  are  described.  The  overall  objective 
of  the  program  was  to  carry  out  a  fundamental  study  of  nucleation  and  film  growth 
mechanisms  in  heteroepitaxial  semiconductor  thin  films.  The  specific  technical 
objectives  were  l)  investigation  of  various  aspects  of  the  mechanisms  of  hetero¬ 
epitaxial  film  growth,  to  establish  technical  guidelines  for  the  preparation  of 
better  films;  2)  preparation  of  improved,  high-quality,  device-grade  heteroepitaxia] 
firfSJ°f  \  aDd  GaAs  °n  lnsulatin6  substrates  by  chemical  vapor  deposition  (CVD) 
methods;  3)  development  of  new  or  improved  methods  of  characterizing  heteroepitaxia! 
r.  in  conductor  films;  and  M  design  and  fabrication  of  selected  thin- film  devices 
■c  ing  Advantage  of  the  unique  properties  of  heteroepitaxial  films.  The  program 
involved  both  theoretical  and  experimental  investigations  of  nucleation  and  growth 
mechanisms  and  development  of  improved  techniques  for  film  and  substrate  preparatioi 
end  characterization.  The  CVD  method  of  growing  semiconductor  films  was  emphasized 
because  of  its  importance  in  the  semiconductor  device  industry.  Main  emphasis  was 
on  the  Si-on-Al203  system,  with  attention  also  given  to  the  Si-on-MgAlo0i  and 
C^As-on-Al203  systems.  The  work  was  divided  among  seven  subtasks:  l)  theory  of 
ap  taxy  and  heteroepitaxial  interfaces;  2)  deposition  studies  and  film  growth* 

3j  analysis  and  purification  of  CVD  reactants;  1*)  preparation  and  characterization 
of  substrates;  b)  studies  of  in  situ  CVD  film  growth  in  the  electron  microscope; 

6)  evaluation  of  film  properties;  and  7)  design  and  fabrication  of  special  devices. 
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Principal  technical  accomplishments  include  1)  experimental  discovery  and 

n  JT?Cterl:Zatl0n  f  Slgnificant  anisotropy  in  the  electrical  properties 
i/A^2°3  of  several  crystallographic  orientations,  and  theoretical 

explanation  of  the  effects  in  terms  of  a  model  that  combines  thermally-induced 
stresses  due  to  expansion  coefficient  differences  with  the  piezoresistance  effect  in 

film  lro^inra^°n  /  T**  substrate  orientations,  deposition  temperatures, 

film  growth  rates,  and  carrier  gas  atmosphere  for  optimized  film  properties 

in  the  Si/Al203  and  Si/MgAl^  systems;  3)  identification  of  the  influence  of 
surface-state  conduction  on  measured  electrical  properties  in  Si  films  with 
carrier  concentrations  <10lbcm--j;  4)  identification  of  the  role  of  the  Alo0o 

Pyrolytic  ^composition  of  SiH4  in  the  formation 

l  techniques  for  preparing  surfaces  of  A1203  and  MgAlpO,  substrates  for 

0f  semiconductor  films  i  h  development  of  an  ion-beam 
sputtering. technique  for  preparing  ultrathin  (<500fl)  regions  in  Alo0„  substrates 
for  transmission  electron  microscepy;  7)  observation  of  the  in  sitd  growth  of 
Si  by  physical  vapor  deposition  (PVD)  and  CVD  on  A1203  and  a^s^s  carbon 
substrates  in  the  electron  microscope;  8)  observation  and  characterization  of 
!£n+t«JnSP°r!1  of  P^?inJe^ted  electrons  through  single-crystal  A120,j  9)  measure- 
^  ^ft?arrief.?;lfet^rae3  in  Si  heteroepitaxial  films  on  A120  and2correlation 
of  lifetimes  with  various  experimental  parameters;  10)  determination  of  surface- 
state  density  distributions  in  Si  heteroepitaxial  films  on  Alo0-j;  and  11)  the 
successful  fabrication  in  Si/AlpO,  of  charge-coupled  devices  exhibiting  good 
charge-transfer  efficiency  at  high  frequencies.  ■  Details  of  these  and  other 
results  and  investigations  are  given  by  means  of  data  tabulations,  graphs,  photo- 
t?arratiVe*  extensive  bibliography  of  electron  microscope  in  situ 
aerate  parts"  ®  3tUdle3  13  3130  lnolu,ied-  The  report  is  boundTn-U^ee 
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objective  'of  °the  iVoIZZs'TlIrry ZITf**!  *“***'  *he  overall 
«nd  film  growth  mechanisms  in  heteroepitaxia/se^"^1  8tUdy  °f  nucleati°n 
The  specific  technical  objectives  1  !  semiconductor  thin  films, 

of  the  mechanisms  of  heteroepitaxial  mi  1 nvestigation  of  various  aspects 
guidelines  for  the  preparation  of  hpt^  growth,  to  establitjh  technical 
high-quality,  device-grade  het  nam  ^  f ±^Tf '  2  PreParatl°n  of  improved, 
ing  substrates  by  chemical  vapor  deD0^^  <rZ,°f  Si  3nd  GaAs  on  *»'**tL 
°f  new  or  improved  methods  of  character! ,! °  ^V°  raethods;  3)  development 
films;  and  4)  desig  r^d?.bMci““tjr«2cttdes2imi*i  5ric™d“ct»^ 
advantage  of  the  unf^  proper tie.  of  het«r“pIta»aTf nL 
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film  and  substrate  preparation  *  nd  develoPment  °f  improved  techniques  for 
growing  aenlco^^tor  a“a  was  elha^LCHtehriZaC1°"-t  Ihe  CTD  "eth“d  °£ 
semiconductor  device  industry  Main  emn^  becau8e  of  its  importance  in  the 
«*t|»  attention  also  given  ihel  o  ^  W"0n  the  Si^-^2°3  W-, 

T^e  work  was  divided  among  se^en  ILtZl* 2?4  aad  CaAs-on-Al  0  2ys3tems. 
epitaxial  Interfaces;  2)  deposition  studies  a  c£  op^-taxy  and  hetero¬ 

purification  of  CVD  reactants-  4)  nrei  ar^t-t  ^  8rowth;  3)  analysis  and 
Strates;  3)  studies  of  ir  situ  CTOfi  lm  ^  ^  oW«tci'Izatlon  of  sub- 

6)  evaluation  of  film  p'EFpFFtles-  and  7^°”^  ln  the  electron  microscope; 
devices.  Properties,  and  7)  design  and  fabrication  of  special 

Principal  technical  accomplishments  inrlnrfp  i\  . 
empirical  characteriration  of  alg“fleajt  *  ^  aXpCrl”e"tal  dlecovery  end 
properties  of  Si/Al.O  films  0f  feveral  *  ands°JroPy  in  the  electrical 
theoretical  explen.hJn  of  the efflcL £e™ T a"d 
thermally-induced  stresses  due  fn  !  <  of  3  model  that  combines 

the  piezoresistance  ef fee  in  Si  diff«ences  with 

orientations,  deposition temoer^uiL  m*3^011  °f  preferred  substrate 
atmosphere  for  optimized  film  properties  in  the^! /I^n8’  3nd  Carrier  8as 
ayatems;  3)  identification  of  the  influence  of  £  2°3  and  Si/MgAl  0 

0ea?firedielect:rical  Properties  in  Si  fiill  jlu surfa^e-state  conduction  on 
4)  ldentl?lcatlon  of  the  role  of ^he  r“nCratl°ns 

for  the  pyrolytic  decomposltlon  of  SiH  in  th^  fAl2°2<SUrf?Ce  aa  a  cataly3t 
5)  development  of  gas-phase  etchlno  $  -t  the  forma^ion  Si  films  by  CVD; 
Qlquee  for  prep.rlS8\^fa“  a  ofS  f  P°Ushl"«  «a»- 

epitaxial  growth  of  semiconductor  f  As  •  6)  8ubstrates  for  ^  hetero- 

sputtering  technique  for  prepar^g  ultrlth  n  f^nnir^  f  an4ion~b^ 
strates  for  transmission  elect-onmir—Ir  ^<500A)  regions  in  AJ  0  sub¬ 
growth  of  Si  by  physical  vann/d»  5y!  P  observation  of  the  In  situ 

carbon  substrates  in  the  electron^nicrosr  CVD  °D  A12°3  and  am°rphous 

tion  of  the  transport  of  photoiniected  e]er^’  observation  anl  characteriza- 

9)  measurement  of  carrie r  lifetimes  in  Si  bp^°nS  tbrough  slngle-crystal  A1  0  ; 

iiretlmes  in  Si  heteroepitaxial  films  on  AljO  a£d3 
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correlation  of  lifetimes  with  various  experimental  parameters;  10)  determina¬ 
tion  of  surface-state  density  distributions  in  Si  heteroepitaxial  films  on 
Al^O^;  and  11)  the  successful  fabrication  in  Si/A^O-  of  charge-coupled 
devices  exhibiting  good  charge-transfer  efficiency  ac  high  frequencies . 

Details  of  these  and  other  results  and  investigations  are  given  by  means  of 
data  tabulations,  graphs,  photographs,  and  narrative.  An  extensive  biblio¬ 
graphy  of  electron  microscope  in  situ  film  nucleatioi  and  growth  studies 
is  also  included.  The  report  is  bound  in  three  separate  parts. 


PROGRAM  SUMMARY 


The  overall  objective  of  this  three-year  program  was  to  carry  out  a  funda¬ 
mental  study  of  nucleatlon  and  film  growth  mechanisms  in  heteroepitaxiai 
semiconductor  thin-film  systems  which  would  lead  to  new  knowledge  and  under¬ 
standing  of  these  processes,  and  then  to  apply  the  results  to  the  preparation 
of  improved  semiconductor  thin  films  and  thin-film  devices  on  insulating  sub¬ 
strates  . 

The  specific  technical  objectives  were  the  following:  1)  investigation  of 
various  aspects  of  the  mechanisms  of  heteroepitaxiai  film  growth,  to  establish 
technical  guidelines  for  the  preparation  of  better  films  which  could  be  applied 
to  real  situations;  2)  preparation  of  improved,  high-quality,  device-grade 
heteroepitaxiai  films  of  Si  and  GaAs  on  insulating  substrates  by  chemical  vapor 
deposition  (CVD)  methods;  3)  development  of  methods  of  characterizing  hetero¬ 
epitaxiai  films  as  to  their  suitability  for  subsequent  device  fabrication;  and 
4)  design  and  fabrication  of  selected  thin-film  devices  which  take  advantage 
of  the  unique  properties  of  such  films. 

The  plan  for  accomplishing  these  objectives  involved  the  study  of  the  funda¬ 
mentals  of  heteroepitaxiai  semiconductor  film  growth  on  insulating  substrates 
as  the  primary  activity,  with  specialized  device  fabrication  used  as  a  means 
both  of  evaluating  film  properties  and  of  exploiting  certain  unique  properties 
of  heteroepitaxiai  semiconductor- insulator  systems.  Both  theoretical  and 
experimental  investigations  were  involved.  The  theoretical  studies  consisted 
of  two  types:  1)  direct  interaction  with  the  experimental  program,  involving 
data  analyses,  suggestion  of  definitive  experiments,  and  postulation  of  specific 
models  to  explain  experimental  observations;  2)  development  of  original  contri¬ 
butions  to  the  theory  of  heteroepitaxiai  growth.  The  experimental  investigations 
were  also  of  two  types:  1)  fundamental  explorations  to  delineate  mechanisms  and 
general  empirical  principles  of  the  heteroepitaxiai  growth  process;  2)  practical 
device  studies  accompanying  the  fundamental  investigations,  so  that  new  develop¬ 
ments  could  be  applied  to  the  improvement  of  films  and  thin-film  devices. 

The  chemical  vapor  deposition  (CVD)  method  of  growing  semiconductor  films  was 
emphasized  because  of  its  importance  in  the  semiconductor  device  industry. 

Main  emphasis  was  on  the  Si-on-A^J.  system,  with  attention  also  given  to  the 
Si-on-MgAljO^  and  GaAs-on-A^O-  systems.  The  work  was  divided  among  seven  sub¬ 
tasks:  1)  theory  of  epitaxy  ana  heteroepitaxiai  interfaces;  2)  deposition  studies 
and  film  growth;  3)  analysis  and  purification  of  CVD  reactants;  4)  preparation 
and  characterization  of  substrates;  5)  studies  of  in  situ  CVD  film  growth  in 
the  electron  microscope;  6)  evaluation  of  film  properties;  and  7)  design  and 
fabrication  of  special  devices. 

The  program  was  carried  out  primarily  at  facilities  of  the  Electronics  Group 
of  the  Rockwell  International  Corporation,  by  Rockwell  personnel.  Parts  of 
three  of  the  specific  subtasks  were  performed  by  personnel  of  the  University 
of  California  at  Los  Angeles  (UCLA),  in  the  Department  of  Electrical  Sciences 
and  Engineering  and  the  Chemistry  Department.  Work  on  another  subtask  was  done 
in  part  in  the  Department  of  Chemistry  of  California  State  University,  San  Diego 
(CSUSD).  Both  the  UCLA  and  the  CSUSD  programs  were  supported  by  subcontracts 
from  Rockwall. 
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films  on  A^O-.  Excellent  agreement  was  obtained  between  theory  and  experi¬ 
ment  in  the  first  two  cases;  the  amount  of  anisotropy  and  the  directions  (in 
the  plane  of  the  film)  of  maximum  and  minimum  carrier  mobility  were  correctly 
predicted.  In  (lll)Si  (on  two  different  A1„0„  substrate  orientations),  however, 
i he  theoretically  predicted  effects  were  too  small  to  account  for  the  experi¬ 
mental  observations,  indicating  that  other  (not  yet  identified)  phenomena  are 
more  important  in  determining  anisotropy  in  the  (lll)Si/Al203  system. 

Extensive  calculations  were  made  for  the  general  (xxl)Si  orientation  (0<x<°°) , 
which  defines  all  Si  orientations  along  the  zone  that  includes  (001),  (111), 

(221),  and  (110)Si  and  thus  all  of  the  Si/Al_0„  epitaxy  inodes  of  major  interest. 

For  n-type  Si/Al20  ,  the  predicted  stress  effects  ranged  from  a  -30%  reduction 
Fnr  (001) Si  to  a  -25%  enhancement  in  maximum  mobility  near  (110) Si.  For  p-type 
films,  the  predicted  stress  effect  is  always  a  mobility  enhancement,  with 
the  maximum  mobility  ranging  from  ~l.()7p  for  (001)Si  to  ~2.4p  for  (110)Si, 

/here  Pq  is  the  zero-stress  mobility.  I?  is  significant  that  ?he  Si  film  orienta¬ 
tion  most  used  in  commercial  devices  -  (001)  -  is  the  one  showing  the  lowest 
mobilities  of  all  those  investigated.  Thus,  considerable  improvement  in  perfor¬ 
ce  of  certain  types  of  devices  could  be  realized  by  exploiting  these  predicted 
^  acts,  which  have  major  significance  for  heteroepitaxial  device  technology. 

Stih task  2:  Deposition  Studies  and  Film  Preparation.  A  major  part  of  the  work 
of  this  subtask  consisted  of  preparing  a  variety  of  Si  and  GaAs  heteroepitaxial 
'  xlm  samples  for  use  in  other  parts  of  the  program.  Experimental  investigations 
of  the  effects  of  various  deposition  parameters  upon  the  properties  of  Si  and 
CcAs  films  continued  throughout  the  program,  including  examination  of  the  follow¬ 
ing:  1)  dependence  of  electrical  properties  on  growth  temperature,  growth  rate, 
and  crystallographic  orientation  of  substrate  (including  the  (1120)  orientation, 
not  previously  used  for  heteroepitaxy  studies);  2)  variations  in  Si  film  prop¬ 
erties  with  thickness;  3)  formation  and  properties  of  p-type  Si/Al„0„  films; 

4)  effects  of  autodoping  (at  temperatures  above  -1050C)  in  Si  films  on  A1„0  and 
MgAl20. ,  5)  effects  of  reactor  configuration  on  film  properties;  6)  characteristics 
of  early-stage  growth  of  Si  films  on  A1  0  ;  7)  growth  of  Si  films  by  SiH  pyrolysis 
at  reduced  pressures  (1  to  10  torr)  ;  0)  growth  of  Si  films  in  gaseous  atmospheres 
other  than  pure  H2;  9)  effects  of  annealing  during  growth  on  properties  of  Si 
films  (no  significant  improvement  in  film  properties  observed) ;  and  10)  growth 
of  GaAs  films  on  and  MgAl20^.  These  studies  revealed  the  strong  inter¬ 

relationships  that  exist  among  the  various  parameters  involved  in  optimiiing  Si 
growth  on  insulators.  Evaluation  of  the  electrical  properties  of  Si  films  on 
A1„0  demonstrated  that  growth  conditions  (ll  must  be  optimized  for  the  particular 
substrate  orientation  chosen;  (2)  differ  for  those  A1203  orientations  which  pro¬ 
duce  the  same  Si  orientation;  (3)  are  dependent  upon  reactor  geometry  and  gaseous 
atmosphere;  and  (4)  must  be  optimized  for  the  particular  film  thickness  desired. 

Based  on  electrical  properties  of  the  films  grown,  the  preferred  substrate  for 
Si  heteroepitaxy  was  identified  as  one  that  produces  (111) Si  growth  rather  than 
(100) Si  growth,  namely  ~4l20)Al20  or  (lll)MgAl  0,,  the  highest  mobilities 
being  obtained  on  (lll)MgAl20 , .  Optimum  growth  temperatures  varied  with  sub¬ 
strate  orientation:  (100)Si_on  (0112)A170,,  1050-1075C;  (lll)Si  on  (10l4)Alo0„ , 
1075-1100C;  (lll)Si  on  ~(1120)A1203,  1075-1100C;  and  (lll)Si  on  (lll)MgAl„0? ,3 
-1025C.  The  preferred  growth  atmosphere  appeared  tc  be  either  H  or  a  He^i 


"MfU'e.f°r  Sr°"th  “n  (0li2)Ai2°3  ot  ~(1120)A1  0,,  but  both  (10l4)Al,O  and 
hiah'n^??4  p??arentJ;y  need  a  He_rich  atmosphere  (-90%  He-10%  H  )  foAbtaining 
w’V  h  /fLl!/Si  F°r  essentially  evefy  orientation  studied,  growth  rates  from 

on  (0ll2)M  04M:  m  n  f°Und  Safficient>  although  good  growth  was  also  achieved 

(0112)A120  at  lower  rates.  These  growth  conditions  are  optimum  for  the  reactor 

system  used,  tfut  should  be  useful  guides  for  CVD  film  growth  in  systems  with  other 
Si/Ai°^  8^etries‘  ?hese  observations  should  lead  to  definite  improvements  in 
C./V72^3  £  ,  ProPe.rties  and  thus  to  improved  performance  and  reliability  in 

bi/AJ^O^  devices  and  circuits.  1 


j_ubtask  3;  Analysis  and  Purification  of  CVD  Reactants.  Since  the  impurity  content 
of  the  various  reactants  used  for  CVD  of  Si  and  GaAs  heteroepitaxial  films  problblv 

selve^^  I  I631  llmil:atl0a  °n  the  achievable  impurity  levels  in  the  films  them- 
ti^nF Uadfrtaken  to  attempt  to  identify  and  establish  the  concentra- 

pioperties  ifPthe  Pa  fn  the  reactanta  that  might  influence  the  film 

P  Pities  if  they  were  to  become  incorporated  into  the  films  during  growth.  Dur- 

ng  the  first  year,  techniques  of  gas  chromatography  were  developed  for  analysis 
the  reactants  used  for  Si  and  GaAs  heteroepitaxy  by  CVD,  with  silicone  oil  and 

PeremohsC°lTS  ^  ^  the  chromatoSraPby •  Several  extraneous  impurity  peaks 

tivelv  idenHf|na  e  ams  o£  SiH4  samPlea»  aad  diborane  (B_H,)  was  tenta- 

identified  as  a  significant  impurity  (-10  ppm),  although  noAonfirmed  by 

spectrometer  techniques.  Small  quantities  of  purified  SiH.,  free  of  diborane 

injections  in  the  chronograph,  bu^the 

were  too  small  for  use  in  laboratory  CVD  experiments.  Beginning  in  the  second  year 

CVD  e6  P^a8ra^‘,  samPles  of  SiHA  and  oE  t  rime  thy  lgallium  (TMG)  used  for  Si  and  Jas 
CVD  experiments  were  analyzed  Tor  impurity  content  by  sensitive  mass  spectrometric 
techniques.  DisLLane  and  trimethylsilane,  together  with  several  othe/impurities 
of  less  concern,  were  found  in  the  SiH4  samples.  impurities 


Significant  Impurity  concentrations  in  some  of  the  reactants  (especially  SiH  ) 

on  s^HI  fc?e  accuracy  of  the  study  the  effects  of  deposition  parameters 

on  Si  film  properties.  Cooperative  efforts  with  vendors  for  preparation^ 

improved-purity  reactants  continued  throughout  the  program,  as  did  analyses  of 
reactants  by  mass  spectrometer  techniques.  It  was  made  Aar  that  significant 

from'  C  D  °ccurred  £n  r-he  purity  of  reactant  materials;  not  only  did  they  vary 
m  supplier  to  supplier  but  also  different  tanks  of  the  same  material  from 
the  same  supplier  were  not  consistent  in  purity.  Another  problem  was  the  lack 

Iheafairr\  f"^  analyS6S  suPPliad  by  different  analytical  laboratories  for 
th  ame  tank  of  reactant  and  the  differences  detection  limits  for  the  same 
f  -umfoui.u  that  different  laboratories  possessed,  making  it  difficult 
to  determine  which  results  were  the  most  reliable  and  which  supplier  of  gases 
was  best.  A  use  test  still  appears  to  be  the  most  reasonable  way  of  evalu- 
ating  materials  tor  the  intended  application.  The  analyses  demonstrated  the  ulti- 

ind  tlTt  <  UnderStandinS  of  the  role  of  impurities  in  epitaxial  film  growth 
3"  I  minimum  requirements  for  the  analytical  methods  that  should  be  used  in 
detecting  these  impurities.  An  extensive  collaborative  study  involving  the 
reactant  supplier,  the  analytical  laboratory,  and  the  ultimate  user  is  needed. 


6^ine  S°me  °f  the  fundamentals  of  the  chemistry  and  reaction  kinetics  of 
fhe  rvn  processes  used  for  growing  heteroepitaxial  films  of  Si  and  GaAs,  Investi- 
gations  of  the  reactions  involved  in  the  formation  of  Si  by  SiH.  pyrolysis  and 
“  tbe  trimethyl  gallium  <TM5)-AsH3  reaction  were  uLert^Se  i^u- 


ence  of  the  Al^  ^rjrfa-e  u  t'^  iituehanism  uf  decomposition  of  Sia"’ as Tf Unction 
mpera  ure,  and  the  decomposition  modes  and  reaction  products  of  TUG  and 
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AsH3  as  functions  of  temperature,  concentration  ratio,  and  H_  partial  pressure 
(to  attempt  to  determine  the  requirements  for  formation  of  optimum-quality  GaAs) 
were  examined.  It  was  found  that  the  surface  of  single-crystal  A1„0„  does 
catalyze  the  thermal  decomposition  of  SiH  ,  whereas  the  pyrolysis  is3homogeneous 
over  a  Si  surface.  No  effect  of  preheating  the  A1„0-  to  temperatures  as  high 
as  1200C  was  found.  A  variety  of  observations  made  on  the  TMG-AsH-  system 
included  the  following,  among  others:  1)  GaAs,  does  result  from  thereaction  of 
IMG  with  AsH3,  and  an  excess  of  AsH  should  be  used  in  any  application  of  the 
process;  and  2)  both  CH^  (methane)  and  U  are  produced  in  the  series  of  step¬ 
wise  reactions  leading  to  the  final  product.  Further  study  of  the  reaction 
at  high  temperatures  is  recommended. 

Sub  task  4: — Preparation  and  Characterization  of  Substrates.  It  was  demonstrated 
that  A1  03  surfaces  prepared  by  mechanical  polishing  techniques  and  used  rou¬ 
tinely  for  semiconductor  heteroepitaxy  typically  had  severe  surface  and  subsur¬ 
face  damage,  with  many  scratches  often  several  microns  deep  yet  often  rendered 
invisible  to  close  inspection  because  of  amorphous  or  fine-grained  debris 
embedded  in  the  scratches  in  the  final  polishing  stages.  Early  in  the  second 
year  a  much  improved  technique  for  mechanical  polishing  of  (1014)A190„  was 
developed,  and  very  good  surfaces  in  this  previously  troublesome  orientation 
were  then  obtained.  Gas-phase  etching/polishing  procedures  using  SF,  and  various 
fluorinated  ha^ocarbons  in  the  1350  to  1500C  temperature  range  were  found  to 
produce  essentially  scratch-free  surfaces  on  (0112)  and  near-(1120)Al  0  sub¬ 
strates.  Extensive  gas-phase  etch-rate  data  were  obtained  as  a  function  of 
crystallographic  orientation  in  this  temperature  range.  The  technique  was 
further  developed  for  (1)  thinning  Al^O  substrates;  (2)  evaluating  the  effects 
of  prolonged  etching  on  (0112),  (0001J,  and  ~(1120)A1  0  ;  and  (3)  assessing 
the  subsurface  damage  caused  by  various  mechanical  polishing  procedures. 
Evaluation  of  mechanical  polishing  methods  for  MgAl  0.  surfaces  indicated  that 
surface  fill-in  occurs  for  this  material,  just  as  forVuO,.  Some  exploratory 
gas-phase  etching  experiments  with  MgAl  0,  surfaces  were  also  carried  out  during 
the  second  year  of  the  program.  z  ^ 

Ion-beam  sputtering  techniques  were  developed  for  preparing  ultra-thin 
(~200A)A1203  wafers  for  use  as  substrates  in  the  in  situ  CVD  experiments  with 
Si.  Wafers  successfully  thinned  to  ~50pm  or  less  by  mechanical  polishing  tech¬ 
niques  were  subsequently  thinned  by  ion  etching  to  the  point  of  perforation 
in  some  areas,  resulting  in  adjoining  regions  of  thicknesses  suitable  for 
transmission  electron  microscopy  as  applied  in  the  in  situ  experiments.  Three 
different  Alo0  orientations  were  successfully  thinned  by  this  method  -  (0001), 
(1014),  and  (0112) .  Considerable  study  of  properties  of  the  resulting  thinned 
substrates  was  carried  out,  and  improvements  in  the  ion-thinning  process  were 
realized  during  the  final  year  of  the  program,  when  the  thinned  substrates  were 
used  in  the  in  situ  CVD  experiments  (SubtasU  5). 

Mechanical  lapping  and  polishing  methods  that  produce  good  quality  surfaces 
suitable  for  use  as  substrates  for  epitaxy  were  developed  during  the  contract 
for_several  orientations  of  Al.Oj  -  (0001),  (0112),  (1014),  (1120),  ~6  deg  off 
(1120)  and  (1122).  By  means  of  etch-rate  techniques  developed  for  this 
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ma5e!lai  Xt  *as  possible  to  determine  the  apparent  depth  of  damage  in  A1„0 
substrate  wafers  at  various  stages  of  preparation.  At  a  given  stage  in  £he 
processing,  the  damage  depth  was  found  to  increase  for  various  orientations 
m  the  following  order:  (1014),  (1120),  (0ll2) ,  and  (0001). 

A  technique  utilizing  the  photoelectric  process  in  a  metal- insulator- 
semiconductor  structure,  consisting  of  an  A1  0  substrate  with  a  Si  or  CaAs 
iim  grown  on  one  surface  and  a  semitransparent  A1  film  on  the  other,  demon- 
strated  that  photoexcited  electrons  from  either  the  semiconductor  or  the 
metal  film  can  be  transported  through  single-crystal  A1_0,  of  several  mils 

MWk  StUdy  °f  tHe  chai‘8e  transport  proceIs3also  established 

that  (1)  the  A1203  used  for  substrates  for  growth  of  heteroepitaxial  semi- 
c  nductor  films  has  trap  levels  approximately  0.18  eV  below  the  conduction 
band  (2)  the  transport  of  photoinj ec ted  electrons  occurs  through  the  Poole- 
Frenkel  conduction  mechanism;  and  (3)  the  quantum  efficiency  for  the  photo¬ 
electric  process  is  quite  low,  approximately  10“5. 

Routine  characterization  of  substrate  surfaces  at  various  stages  of  prepara¬ 
tion  continued  throughout  the  program,  utilizing  various  standard  techniques 

scanning,  Microscopy?  “lffraCtl°"  ^  tZZ  Ing 


Subtask  5: 


Studies  of  in  situ  Film  Growth  in  the  Electron  Microscope. 


In  tae. first  year  ot  the  program  many  of  the  modifications  required  in  the 
transmission  electron  microscope  for  in  situ  observation  of  the  nucleation 
and  early-stage  growth  of  CVD  semiconductor  films  on  insulating  substrates 
were  completed.  Early  in  the  second  year  a  series  of  electron  microscope 
modifications  and  tests  was  finished,  culminating  in  the  first  series  of 
successful  PVD  experiments  inside  the  electron  microscope.  A1  was  deposited 
onto  a  heated  carbon  substrate  and  a  sequence  of  micrographs  was  taken  during 
the  growth  process,  demonstrating  the  feasibility  of  performing  in  situ 
nucleation  and  growth  studies  in  the  equipment.  Additional  in  situ  PVD 
experiments  were  carried  out  in  the  second  year,  with  both  A1  and  Au  deposited 
onto  amorphous  carbon  substrates  to  delineate  further  the  required  techniques 
and  experimental  problems  to  be  encountered  in  the  later  CVD  experiments, 
a  culations  and  design  for  the  CVD  microchamber  were  also  completed  during 

the  second  year,  and  the  fabrication  of  the  microchamber  and  associated 
hardware  was  begun. 


During  the  final  year  of  the  program  the  fabrication  of  the  CVD  microchamber 
and  its  mounting  flange  was  completed,  and  a  gas-handling  manifold  was 
installed  on  the  electron  microscope.  Gas  flow  experiments  were  performed 
to  determine  the  flow  rate  of  gas  through  the  microchamber  as  a  function  of 
pressure  and  to  determine  the  maximum  pressure  attainable  in  the  microchamber 

In  addition,  a  number  of  iji  situ  Si  CVD  experiments  were  performed  resulting 
m  the  successful  growth  of  Si  films  in  the  electron  microscope  by  Hie 
pyrolysis  of  SiH^  The  nucleation  and  early  growth  of  Si  on  both  amorphous 
carbon  and  single-crystal  A^O-  substrates  was  observed,  leading  to  the 
following  conclusions:  1)  the  pyrolysis  of  SiH,  to  form  Si  films  by  CVD 
inside  the  electron  microscope  is  feasible;  2)^the  in  situ  study  by  trans- 
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mission  electron  microscopy  of  the  post-nucleation  and  early  growth  stages 
of  a  semiconductor  film  grown  by  CVD  is  feasible;  3)  the  nucleation  and8 
early  growth  processes  for  CVD  Si  are  fundamentally  similar  to  those  of 
metal  films  grown  by  PVD,  although  some  specific  differences  exist;  4) 

nucleroTf31  f1  grrVn  (0112)A1203  results  P^arily  from  the  growth  of 

and  not  T  °rief ation  at  ««  “P^se  of  randomly  oriented  nuclei, 

and  not  from  the  large  nucleation  rate  of  these  favorably  oriented  nuclei; 

IrM  ±CTbeam  sPuttering  process  can  produce  electron-transparent  A1„0„ 


...I,  r  w  ‘  -  Li.  uu-Li  duaudi  ent 

suitable  for  use  as  substrates  for  in  situ  CVD  film  growth  experiments?  3 


— — -task  6: — Evaluation  of  Film  Properties.  From  the  beginning  of  the 
program,  routine  evaluation  of  film  properties  was  carried  out  by  established 
and  electriXiry  electron  diffraction  analysis,  metallographic  analysis, 
evaluating  th  “®asuremeJtl 3  °f  transport  properties.  A  new  technique  for 

the,Cha,raC  eriStiGS  °f  the  ^terfacial  region  of  heteroepitaxial 
film,  was  developed,  involving  measurements  of  photoelectron  emission  from 
mono chromatically-illuminated  films  in  the  MIS  configuration  on  insulating 
U  es  (see  Subtask  4).  Photocurrents  due  to  electron  transport  through 


the  single-crystal  A1  O3  substrates  were  measured  as  a  func^  h l?* 

P.nprcrv.  nnrmUHt,n  — j x  j  r-  .  >  LWU 


/  j  - - «  X  UUL  L  i-Ull  ui  Ull 

nergy,  permitting  determination  of  various  parameters  in  the  Si/Al  0  and 

M  Q1tfl«>nnin  mi _ _  ^  2  3  °l*U 


GaAs/Al  0  systems, 
affinity  of  A1  0 


These  measurements  gave  values  of  l.OeV  for  the  electron 

4  50eV  for  f°rcthe  bafrier  height  at  an  A1  0  -A1  interface, 

.50eV  for  the  Sl-Al  03  interface  barrier  height,  ~0.37eV2f3r  the  band- 

bending  an  Si  near  tfie  Si-Al  interface,  -O.lOeV  for  the  band-bending  in 

12^!  f  5  t12°3  ln^crface,  and  electron  escape  lengths  of  at  least 

pm  in  Si  and  23pm  in  GaAs .  Values  for  the  work  function  of  various 
metals  were  also  determined  by  these  measurements. 


Determination  of  the  energy  spectrum  of  back-scattered  proton  or  alpha- 
particle  aeams  injected  in  channeling  directions  in  heteroepitaxial  semi¬ 
conductor  films  was  investigated  as  a  means  of  measuring  the  density  and  the 

r  SMTtUwral  d6feCtS  in  the  filmS‘  Experiments  indicated'that 
i/insulator  fiims  have  less  imperfect  interfacial  regions  than  do  GaAs/ 

films  £n\oiTpUl  n  fSt  structures  of  those  examined  were  found  in  (100)Si 
films  on  (0112)A1  0  substrates  and  in  (lll)Si  films  grown  on  near-(1120)Al  0 

I"*"-?'1™  early  growth  stages  of  SI  on  ALG,  was  also2  3 

btained  by  conventional  transmission  electron  microscopy  of  very3thin  Si 
films  grown  on  ion-thinned  A1203  substrates  in  a  conventional  (aLospheric- 
pessure)  vertical-flow  CVD  reactor.  These  experiments  showed  that  the  growth 
of  a  single-crystal  Si  film  by  CVD  is  the  result  of  coalescence  processes  ^n 
filmsaw^h8r°Wt?  fnd  n0t  °f  nucleation  Phenomena  alone,  producing 

easier  mobiHti::  ^  °f  defeCtS  and  relativel*  lower 


Rapid  acquisition  of  data  on  electrical  properties  of  the  films  was  very 

and°third  vLrs-  C°ndUCt  °f  ^  exPeriments-  especially  in  the  second 

and  third  years,  measurements  of  film  conductivity  type,  resistivity,  carrier 

eDiJr.  l  n’iand  carrier  mobility  were  made  routinely  on  a  majority  of  the 

methnHlal  h3™13163  PrePared  on  the  program,  utilizing  either  the  van  der  Pauw 
thod  or  the  more  accurate  and  conventional  Hall-effect  bridge  method. 

lino/  eSSent\ial  for  the  study  of  the  effects  of  changes  in  depo¬ 

sition  parameters  on  Si/Al^  and  Si/MgAl^  film  properties,  and  provided 
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considerable  insight  into  the  factors  which  most  strongly  influence  film 
quality  so  that  identification  of  the  conditions  for  optimized  film  growth 
could  be  made. 

It  was  determined  that  the  electrical  properties  of  undoped  n-type  hetero- 
epitaxial  Si  films  grown  on  various  orientations  of  A1„0  (and  also  MgAl  0.) 
by  the  pyrolysis  of  SiH^  are  dominated  by  surface-state  conduction  for  2  ^ 
carrier  concentrations  of  ~10-*-°cm“  or  below.  It  was  found  that  there  were 
inhomogeneities  in  the  donor  concentration  of  typical  CVD  Si/Al  0  films 
over  the  film  area,  and  that  a  concentration  gradient  existed  from  the  center 
of  the  susceptor  radially  outward;  films  reflected  this  variation  depending 
upon  the  placement  of  the  substrate  on  the  susceptor  during  CVD  growth.  Gas 
flow  characteristics  or  a  non-uniform  temperature  of  the  rf-heated  pedestal 
(susceptor)  were  thought  to  account  for  the  effect.  Measurements  were  made 
of  the  variation  of  electrical  properties  of  Si/Al  0  with  temperature, 
and  some  of  the  observed  effects  were  attributed  to  high  defect  densities 
(e.g.,  deep-lying  donor  levels)  or  inhomogeneous  strains  in  the  films. 

The  most  significant  development  to  come  from  the  film  evaluation  procedures 
was  the  observation  of  the  anisotropy  in  electrical  properties  in  Si/Al,  0 
films.  Mobility  measurements  as  a  function  of  azimuthal  direction  (every3 
18  deg)  in  the  film  plane  indicated  a  maximum  mobility  in  two  directions  and 
a  minimum  mobility  in  two  directions,  the  latter  displaced  by  90  deg  from 
the  former.  The  mobility  anisotropy  factor  A,  defined  as  the  ratio  of  the 
difference  between  the  maximum  and  minimum  values  of  mobility  in  the  plane 
of  the  film  to  the  average  value  of  the  mobility  in  that  plane,  was  found 
to  be  about  40 %  for  (22 1) Si/ ( 1122) Al^CL  and  about  9%  for  (001)Si/0112)Al  0  . 
Results  of  theoretical  calculations  (Subtask  1)  agreed  well  with  the  explr?- 
mental  data.  The  calculations  and  the  experimental  results  indicated  that 
(221) Si  exhibits  higher  electron  mobilities  than  other  more  commonly  used 
orientations.  Measurements  of  anisotropy  at  77K  were  also  consistent  with 
the  corresponding  increases  in  piezoresistance  coefficients  at  that  tempera¬ 
ture.  Data  analyses  predicted  zero-stress  mobilities  significantly  below 
bulk  crystal  values,  however,  indicating  mechanisms  other  than  thermally- 
induced  stresses  were  dominant  in  reducing  carrier  mobilities  in  hetero- 
ep itaxial  films.  Extensive  studies  of  (lll)Si  on  (1120)A1  0  and  on  (10l4)Al.,O 
gave  experimentai  anisotropies  averaging  16%  and  30%,  respectively,  much  2  3 
larger  than  theoretical  predictions,  again  indicating  the  presence  of  other 
major  influencing  factors.  There  appeared  to  be  an  inverse  relationship 
between  anisotropy  in  (lll)Si  and  the  minimum  or  the  average  mobility,  higher 
anisotropy  corresponding  to  lower  mobility.  Attempts  were  made  to  correlate 
these  results  with  reactor  configuration  (i.e.,  horizontal  or  vertical)  and 
various  deposition  conditions.  Measurements  were  also  begun  to  examine 

the  possibility  of  mobility  anisotropy  being  present  in  the  Si/MgAlo0,  system. 

2  4 

The  surface-state  density  of  thermally  oxidized  Si  films  on  A1  0  was  deter¬ 
mined  late  in  the  program  using  the  MOS  C-V  technique.  Evidenceof  both 
donor-  and  acceptor-type  surface  states  was  found;  a  peak  in  the  acceptor- 
state  density  appeared  at  ~0.l6eV  below  the  conduction  band,  but  the  exact 
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location  of  the  donor-state  density  peak  could  not  be  determined.  Measure¬ 
ments  of  high-field  transport  properties  of  Si  and  GaAs  heteroepitaxial  films 
on  A1  03  were  also  undertaken,  early  in  the  program,  to  obtain  drift  mobility 
aata  for  some  of  the  films. 

?-ul)task  ?!  Design  and  Fabrication  of  Devices.  In  the  first  year  of  the 
contract  apparatus  for  determining  minority  carrier  lifetime  by  pulsed  C-V 
measurements  in  MOS  structures  was  designed  and  constructed  and  tests  were 
begun.  A  special  MOS  structure  was  designed  for  measurement  of  channel 
conductance,  high-  and  low-field  transport  properties,  and  various  inter¬ 
face  characteristics  of  heteroepitaxial  films.  Initial  attempts  to  fabricate 
Schottky  barrier  diodes  in  Si/Al „0_  films  as  a  means  of  evaluating  their 
electrical  properties  were  not  successful  and  were  not  pursued  further.  In 
the  second  year  the  design  of  a  Schottky-barrier  type  of  FET  was  completed 
or  use  in  fabricating  experimental  FET  structures  in  GaAs/insulator  films 
for  operation  at  1  GHz.  Most  of  the  device-oriented  effort  centered  about 
the  determination  of  carrier  lifetimes  using  the  MOS  pulsed  C-V  technique 
(work  which  extended  to  the  end  of  the  program)  and  attempts  to  fabricate 
Schottky-barrier  FET's  in  GaAs/Al^. 

Recent  device  efforts  produced  Schottky-barrier  diodes  (in  n-type  Si/Al  0 
samples)  having  good  reverse  but  unsatisfactory  forward  characteristics?  3 
The  Schottky-barrier  FET  structures  were  still  not  satisfactory.  Preliminarv 
work  on  fabricating  and  evaluating  Schottky-barrier  photovoltaic  cells  using 
illumination  from  the  back  side  was  begun,  and  charge-coupled  devices  (CCD's) 
in  Si/Al^O^  composites  were  successfully  designed,  fabricated,  and  tested. 

The  lifetime  measurement  method  used  has  the  important  advantage  that  the 
actual  carrier  lifetime  is  magnified  by  the  factor  N/n.,  where  N  is  the 
impurity  concentration  in  the  semiconductor  and  n  is  the  intrinsic  concen¬ 
tration,  Sq  that  Very  Short  llfetlmes  typical  of  Heteroepitaxial  systems 
(10  -10  •’sec)  could  be  measured.  Carrier  lifetimes  and  values  of  surface 

recombination  velocity  were  obtained  for  As-doped  n-type  Si/Al samples 
and  for  As-doped  n/n  -type  Si/Al20„  samples  grown  by  CVD  on  thiscontract, 
as  well  as  for  some  commercially-obtained  Si/Al  0  samples  that  were  P-doped. 
Lifetimes  for  As-doped  (lOO)-oriented  Si  varied  from  ~10~10  sec  for  films 
l-3pm  thick  to  -5  x  10“9  sec  for  films  ~10ym  thick;  the  particular  P-doped 
films  measured  exhibited  lifetimes  nearly  an  order  of  magnitude  longer. 

dePendfnce  of  lifetime  on  As  doping  concentration  or  on  Si  orientation 
((100)  and  (111)  were_s tudied)  was  detected.  It  was  found  that  an  underlying 
n  layer  (~3  x  101  cm  )  significantly  enhanced  the  lifetime  in  a  3ym  top 
layer  of  As-doped  (-lO^cnf  )  Si  for  a  given  total  film  thickness,  the  enhance¬ 
ment  being  greater  the  greater  the  n+  layer  thickness  (T~lysec  for  18ym  total 
tuckness);  the  n  sublayer  evidently  acts  as  a  "getter"  for  the  trapping 
centers  that  tend  to  lower  the  lifetime  in  the  n-type  material. 

Since  CCD's  had  previously  been  fabricated  in  Si  and  showed  good  charge  transfer 
and  since  other  Si/Al ,0  devices  had  exhibited  good  high-frequency  and  radiation- 
resistant  characteristics  it  was  determined  that  fabrication  of  CCD's  should 
be  undertaken  in  Si/Al^O^.  The  devices  were  successful,  with  good  transfer 
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efficiency  at  high  frequencies;  four-phase  8-mil-per-cell  CCD's  were 
at  2MHz  with  0.99  efficiency.  Low-frequency  operation,  however,  was 
to  be  limited  by  the  short  carrier  lifetime  which  allowed  charge-up 
potential  well.  6  v 
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The  model  described  in  this  paper  uses  an  electronic  wave  function  which  is  defined  to  be 
nonzero  only  along  the  lines  connecting  first  nearest  neighbors  in  the  metafile  lattice.  The 
electrons  arc  assumed  to  move  freely  along  the  lines  between  nearest  neighbors.  No 
electron -electron  or  electron-nucleus  force  is  included  in  the  model  calculations  (except 
for  forces  arising  from  the  Pauli  exclusion  principle).  The  work  function  is  defined  as  the 
amount  of  energy  required  to  move  an  electron  from  a  point  slightly  inside  the  crystal  to 
a  point  slightly  outside.  The  contribution  or  the  electronic  double  layer  is  included  in  the 
calculation  of  the  work  function  as  well  as  the  dependence  of  the  double  layer  potential 
on  the  surface  geometry.  Surface  states,  where  the  electron  is  localized  in  the  neighborhood 
of  the  face  of  the  crystal,  are  found  to  have  energies  sufficiently  above  the  Fermi  level  to 
eliminate  the  possibility  that  they  make  any  contribution  to  the  double  layer  potential  for 
the  case  of  the  (100)  crystal  plane.  Consequently,  surface  states  have  been  ignored  in  all 
the  calculations.  The  surface  double  layer  is  assumed  to  be  caused  by  the  presence  of  a 
finite  potential  barrier  at  the  surface  of  the  crystal.  Bulk  electronic  wave  functions  can 
pene'rate  this  barrier  and  decay  exponentially  outside  the  crystal.  The  only  parameters 
required  by  the  model  arc  the  nearest  neighbor  distance  for  the  lattice  and  the  height  of  the 
potential  barrier  at  the  surface.  The  former  quantity  is  fixed  by  the  lattice  structure 
(body  centered  cubic  for  the  alkali  metals)  and  by  the  density,  while  the  latter  quantity  can 
be  adjusted  to  give  the  best  agreement  between  the  mo  .el  calculations  and  experiment. 
For  the  alkali  metals,  lithium  through  sodium,  the  best  value  of  the  barrier  height  is 
about  50  %  of  the  sum  of  the  ionizatioi:  potential  energy,  the  heat  of  vaporization,  and  the 
calculated  Fermi  level  for  the  corresponding  metal  In  addition,  the  value  of  the  double 
layer  potential  for  sodium  agrees  very  well  with  a  more  sophisticated  calculation  by 
Bardeen  and  is  reasonably  close  to  the  experimental  measurement. 


1.  Introduction 

The  free  electron  molecular  orbital  method  (FEMO)  has  been  successfully 
applied  by  Kuhn1)  and  by  Rudenberg  and  Scherr2)  to  the  treatment  of 
7t-electron  systems  in  conjugated  molecular  systems.  As  in  the  case  of  Hiickel 
theory3),  the  free  electron  orbital  method  is  directed  towards  the  develop- 

*  Research  supported  under  NSF  Grant  #GP-28722X  and  North  American  Rockwell 
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merit  of  a  “relative”  internally  consistent  theory  for  these  complicated 
systems.  In  the  case  of  Hiickel  theory,  parameters  such  as  overlap  integrals 
are  determined  by  experiment  in  one  molecule  and  are  assumed  to  have  the 
same  magnitude  in  others.  Strictly  speaking,  in  the  FEMO  method  there  are 
no  adjustable  parameters  (unless  the  distances  between  atoms  are  allowed  to 
depart  from  the  values  actually  observed  in  the  molecule).  Nevertheless, 
both  methods  provide  a  reasonable  semi-quantitative  picture  of  the  electronic 
state  of  affairs  within  rc-systems,  and  have  been  useful  for  making  semi- 
quantitative  estimates  concerning  the  chemical  and  spectroscopic  properties 
of  n-systems. 

Recently  Montroll  and  his  coworkers  ')  have  expanded  the  FEMO  method 
to  include  the  treatment  of  local  potentials  due  to  individual  atoms.  The 
theory  of  course  remains  highly  schematic,  as  well  as  parameterized,  but  it 
does  allow  relatively  simple  analytical  solutions  to  be  obtained  for  previously 
complicated  problems.  In  this  way.  Montroll  and  his  collaborators  have 
been  able  to  explore  the  qualitative  features  of  many  interesting  solid  state 
phenomena  including  the  electronic  properties  of  defects  and  surfaces. 
Since  the  original  FEMO  method  proved  useful  in  the  semi-quantitative 
discussion  of  the  electronic  properties  of  molecules,  it  is  interesting  to  see 
whether  both  the  FEMO  and  Montroll’s  modification  of  it  can  be  used  in  the 
semi-quantitative  study  of  solids.  We  have  in  mind  not  only  the  examination 
of  schematic  relationships  (Montroll  and  his  collaborators  have  already 
done  that),  but  also  the  development  of  an  internally  consistent  parameter¬ 
ized  theory  in  the  same  spirit  as  Hiickel  theory. 

In  order  to  study  this  question,  we  have  elected  to  apply  the  network  theory 
to  the  estimation  of  the  work  functions  of  monovalent  metals,  especially 
the  alkalis.  Our  principal  objective  has  been  to  determine  the  work  function 
and  its  variation  with  crystal  orientation.  The  approach  is  as  simple  as 
possible,  involving  a  single  adjustable  parameter,  namely  the  difference 
between  the  energy  of  a  free  electron  in  vacuum  and  the  lowest  level  of  the 
conduction  band.  This  parameter  is  invariant  to  crystal  orientation  and 
could  in  principle  be  determined  from  measurements  on  a  single  crystal  plane. 
All  refinements,  such  as  exchange  ai  1  correlation  effects,  are  explicitly 
ignored  but  should  be  implicitly  involved  within  the  barrier  height  which 
constitutes  the  adjustable  parameter  mentioned  above.  The  detailed  explicit 
treatment  of  these  additional  effects  would  be  inconsistent  with  the  severe 
approximations  already  contained  within  the  theory.  One  would  hope, 
however,  that  the  network  model  would  give  a  reasonable  semi-quantitative 
account  of  the  variation  of  work  function  with  crystal  orientation,  once  the 
adjustable  parameter  is  determined,  by  using  the  experimental  data  available 
for  one  crystal  plane.  Furthermore,  one  would  hope  that  there  would  be  an 
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orderly  relationship  between  the  predictions  of  the  theory  for  one  alkali 
metal  as  compared  to  another. 


2.  The  network  model 

Our  treatment  actually  makes  use  of  the  FEMO  method  coupled  with  the 
Green’s  function  approach  developed  by  Montroll  for  the  treatment  of 
defects  and  surfaces  in  connection  with  his  more  elaborate  network  model. 
The  one  electron  wave  function  is  chosen  to  be  nonzero  only  along  the  lines 
connecting  the  first  nearest  neighbors  in  the  metallic  lattice.  This  is  equivalent 
to  requiring  that  the  potential  be  infinite  everywhere  within  the  crystal  except 
on  these  lines.  The  electrons  are  assumed  to  move  freely,  i.e„  we  do  not 
associate  potential  wells  (as  in  the  Montroll  modification)  with  the  individual 
atoms.  To  this  extent,  our  model  is  more  like  the  original  FEMO  method 
However,  we  shall  have  to  make  use  of  the  Bloch  theorem  *)  in  order  to  deal 
with  the  crystal  and  the  Green’s  function  approach  will  have  to  be  used  to 
treat  the  extended  defect  represented  by  the  surface.  The  crystal  symmetry 
enters  the  problem  through  the  nodal  conditions  eqs.  (3a)  and  (3b)  below 
and  through  the  Bloch  theorem. 

The  work  function  is  defined  as  the  amount  of  energy  required  to  move  an 
electron  from  the  point  slightly  inside  the  crystal  to  a  point  slightly  outside 
and  is  effected  by  the  electronic  double  layers)  which  in  turn  depends  on 

crystal  surface.  One  of  our  main  tasks  will  involve  calculating  the  strength  of 
this  doubL;  layer. 

As  long  as  an  electron  remains  on  a  line  joining  two  nearest  neighbor 
atoms  ins»  ,e  the  lattice,  it  has  no  force  acting  on  it,  and  its  potential  energy 
may  be  set  equal  to  zero.  We  assume  that  an  abrupt  change  in  potential 
occurs  as  the  electron  moves  on  to  a  dangling  surface  bond  where  it  ex¬ 
periences  a  constant  potential  VB  (|/„>0).  The  surface  barrier  VB  is  the  only 
parameter  in  the  theory  and  will  have  to  be  determined  through  reference  to 
some  single  experiment.  In  fact,  the  form  of  the  potential  barrier  in  the  real 
crystal  is  not  abrupt,  but  since  we  are  searching  for  the  simplest  possible 
parameterized  theory  (and  also  because  of  the  severe  approximation  already 
contained  within  the  network  model  itself),  we  assume  that  it  is  abrupt 
The  dangling  bonds  for  surface  atoms  have  the  same  direction  as  if  they  were 
joined  to  another  lattice  site;  but  instead  of  being  of  length  R  (where  R  is  the 
nearest  neighbor  distance),  they  extend  all  the  way  to  infinity.  The  wave 
equation  for  this  system  may  be  expressed  in  the  following  form: 

VV  =  (a2-A2)i/^,  (|) 

in  which  «2  =  2 mV!h\  /. 2  =  2mKEjh\  and  ^  is  the  electronic  wave  function, 
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in  the  mass  of  the  electron,  KE  the  energy  of  the  state  described  by  i \ji,  V  the 
potential,  and  h  is  Planck  s  constant  divided  by  2k.  Clearly  V  (and  ot)  equals 
zero  lor  positions  along  bonds  in  the  bulk,  while  V=  K„  on  the  dangling 
bonds.  The  solutions  ofeq.  (I)  for  the  surface  and  bulk  bonds,  respectively, 
are; 

•As  =  «s  exp  [-  (a2  —  A'2)1  .v} ,  Os$v<  +  co,  (2a) 

•An  =  «B  cos  \kx  +  ()„} ,  0<v<«,  (2b) 

Each  ip  is  defined  only  on  a  particular  bond  and  is  zero  elsewhere.  Where 
several  bonds  meet  at  a  lattice  point,  the  wave  functions  associated  with  them 
must  be  equal.  In  addition,  there  are  requirements  associated  with  the  con¬ 
servation  of  momentum  These  restrictions  may  be  expressed  as  follows2): 

*Ai  (0)  =  <M0)  =  ipAO)  =  ...  =  (MO),  (3a) 


i  =  I 


In  these  equations  n  is  the  number  of  bonds  meeting  at  a  particular  lattice  site. 

If  the  quantity  </>(/,  m,  n)  is  the  value  of  the  wave  function  at  the  lattice 
site  (/,  in,  n),  eqs.  (3a)  and  (3b)  require,  for  a  body  centered  cubic  lattice,  the 
following  relationship: 

8  cos(kR)  </>(/,  in,  n)  =  <p(l  +  I ,  m  +  I ,  n  +  I )  +  <p  (I  +  I ,  m  -  |,  «  +  |) 

+  <p  (I  -  I ,  in  +  I ,  n  +  1 )  +  rp  (I  -  l,iu-  I ,  n  +  | ) 

+  (p  (I  +  I,iii  +  I ,  n  —  \  )  +  <p  (I  +  I,//;  —  I ,  n  —  | ) 

+  '/>(/-  I,  m  +  I,  //  -  \  )  +  ip(l  -  I,  />/  -  I,  n  -  I). 

(4) 

In  addition,  for  the  bulk  crystal  we  make  use  of  periodic  boundary  con¬ 
ditions: 


<P(I,  m,  n )  =  (l>(l  +  N,  in,  n )  =  (p  (I,  in  +  N,  //)  =  cP(l,  „  +  N).  (5) 
A  possible  solution  to  eqs.  (4)  and  (5),  which  also  satisfies  Bloch's  theorem,  is: 

J2 


</>(/,  in,  //)  =  A  exp  j2*'  (/.v,  +  ms2  +  //jj)j, 


(6) 


where  v,,  s2,  s2  =  I,  2,  3 . N.  Substitution  of  eq.  (6)  into  (4),  with  appro¬ 

priate  collection  of  terms,  yields: 


cos (kR)  - 


(7) 
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Since  the  energy  of  an  electron  depends  on  *,  and  j„  s2,  and  j,  arc  indices 
denoting  possible  states  of  the  system,  eq.  (7)  determines  the  density  of  states 
for  network.  The  density  of  sets  (j„  s2,  j3),  satisfying  eq.  (7)  in  the  range  d E, 
has  been  derived,  in  another  context,  by  Jelitto7)  (see  tig.  I).  The  density-of- 
states  function  is  symmetric  about  cos  (A /?)=(),  and  it  follows  that  the  Fermi 


Fig  1.  The  density  of  stales,  d/Vn/df,  for  a  body  centered  cubic  lattice  plotted  as  a 
function  of  E,  where  E  is  defined  to  be  the  cosine  of  the  reduced  momentum  wave  vector, 
kR.  This  graph  is  given  by  Jelitto7). 


level  of  the  network  corresponds  to  an  energy  level  such  that  kR  =  \n. 
Therefore,  the  Fermi  level  is  given  by: 

/<,:  =  /f2/32  HlR2  .  (g) 

In  order  to  deal  with  surface  states,  we  assume  that  the  network  is  infinite 
in  both  the  .v  and  y  directions,  but  has  a  finite  thickness  in  the  z  direction. 
It  is  possible  to  write  a  general  relationship  [analogous  to  ?q.  (4)]  which 
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applies  to  either  bulk  or  surface  sites.  Thus,  we  have: 

8cos(AR)  (f>(l,  in,  n)  -  <f>(l  +  I,  in  +  I,  «  +  I)  -  4>(l  +  I,  in  -  I,  n  +  I) 

—  <t>  (I  —  I ,  in  4-  I ,  n  4-  I )  —  (p  (I  —  I ,  hi  —  I ,  n  +  I ) 

-  <i>  (/  4-  I  ,m-\-  I ,  n  -  !)-</>(/  +  I ,  n  -  I ) 

-(/>(/-  I ,  m  +  I ,  ;i  —  \  )  —  (p  (I  —  I ,  in  -  I ,  n  ~  I ) 

=  H  (l,  in,  n),  (9a) 

</>(/+  N,  in,  /;)  =  <p(l,  in  +  N,  n)  =  <J>(I,  in,  n ),  (9b) 

where 

H (I,  m,  n)  =  [<5„  ,  +  <)„,  ,v.]  [4  cos  (A  R)~  4A(kR)]  c/)(l,  in,  n) 

-<V  ,[</>(/  +  I,  m  +  I,  N.)  +  <t>(l  +  I  .in  -  I  ,N.) 

+  <p(l  -  I,  ill  +  I,  (V.)  +  <f>(l  —  I,  ill  —  I,  /V.)] 

—  ^n.  A,  [</*(/  +  I  ,  HI  +  1,1  )  +  </>(/+  I  ,  III  —  I  ,  I  ) 

+  </>(/-  \,  III  +  I  ,  !)  +  </»(/-  I ,  III  ~  1  ,  I  )], 

and  the  <5's  are  Kronecker  delta  functions.  An  equation  of  this  general  form 
was  introduced  by  Montroll  >)  although  he  did  not  specilically  treat  the  case 
in  which  an  abrupt  barrier  of  finite  height  Vn  was  located  at  the  surface. 
We  assume  that  the  z-direction  is  the  (001)  direction  and  that  the  exposed 
surfaces  are  (001)  surfaces  having  coordinates  z=  1  and  z  =  N:,  respectively. 

Because  of  the  finite  extent  of  the  crystal  in  the  z  direction,  we  can  no 
longer  demand  that  <j)  satisfy  the  Bloch  theorem  in  that  direction.  However, 
this  requirement  may  still  be  applied  to  the  .v  and  y  directions.  A  solution  to 
eq.  (9a),  which  satisfies  the  Bloch  theorem  in  these  directions,  is: 

4>{l,  in,  n)  =  A  exp  j2*'  [s,/  +  j2//i]j  0(/i), 

st,s2  =  1,2,3,..  ,N.  (10) 

If  we  substitute  this  into  eq.  (9)  and  collect  terms,  the  resuit  is 

2Fn<P(n)  -  4>(//+  I  )-<P(n-  I  )  =  //„(«),  (II) 

in  which 

W B(«)  =  [<V  i  +  ‘Vn,]  »?^(«)  “‘Vi  ^(^r)  “  <V,v,  ^(1), 

F0  =  cos(AR)/[cos(27ti,/A')  cos(2iti-2/iV)]  , 

//  =  [cos (kR)  -  A  (kR)]l[cos(2nsl/N)  cos  (27ts'2/A  )]  . 
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A  difference  equation  having  the  same  form  as  eq.  (II)  has  been  solved  by 
Montroll4)  using  an  appropriate  Green's  function,  g(n).  In  terms  of  g,  this 
solution  may  be  expressed  as: 

*(»)=  I.  aU-n)Hn(j),  (12) 

j=  i 

where  the  Green's  function  is  defined  as: 


9(») 


3  =  I 


exp(2jm/i/AL) 
f'„  -  cos  (2  tti'/X)’ 


(12a) 


For  the  |Fn|>  I,  </>(/»)  decays  exponentially  as  n  refers  to  sites  further  from 
the  surface,  i.e.,  deeper  in  the  bulk  of  the  crystal.  Therefore,  this  range  of  Fn 
is  the  range  for  surface  states;  and  by  substitution  of  the  solution  in  this  range 
back  into  eq.  (1 1 ),  it  may  be  shown  that  the  following  relation  holds: 

[eos2(£R)  -  A1(kR)]i  =  cos(2tw,/JV)  cos(2iw2/7V).  (13) 

This  equation  is,  for  surface  states,  the  equivalent  of  eq.  (7)  for  bulk  states 
and,  like  eq.  (7),  determines  the  density  of  states. 

The  bulk  wave  functions  may  be  normalized  as  follows.  First,  Hie  square  of 
the  wave  function  (square  of  the  modulus)  is  integrated  over  each  individual 
bond,  yielding  a  quantity  proportional  to  the  fraction  of  electrons  residing 
on  the  bond.  This  fraction  is  then  summed  over  all  bonds  and  the  sum  is  set 
equal  to  unity.  The  fact  that  each  state  is  occupied  by  two  electrons  is  then 
taken  into  account  at  a  later  step  by  multiplying  the  density  of  states  function, 
eq.  (16),  by  2.  The  integral  for  the  square  of  the  wave  (unction  on  a  dangling 
surface  bond  is,  according  to  eq.  (2a): 


o,  =  l(MO)|2  J cxp(-  2/(.v)d.v  =  (14a) 

0 

In  the  above  equation,  /(  =  (x2-a2)*  and  ;  \(0)|  is  the  value  of  the  wave 
function  at  the  surface  node.  The  integral  for  the  square  of  the  wave  function 
on  a  bond  connecting  two  nearest  neighbor  lattice  sites  is,  from  eq.  (26): 

d. = j  cos*  + s.)  J  {« + si"  i2  (t  *  - sin  <“•>} ,  „4W 

where  R  is  again  the  distance  between  nearest  neighbors.  It  should  be  noted 
that  although  both  the  above  integrals  depend  on  the  kinetic  energy  of  the 
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electron,  Ds  is  also  a  function  of  the  surface  barrier  height,  while  Z)„  is  a 
function  of  the  distance  R.  In  a  body  centered  cubic  lattice,  each  site  has 
eight  first  nearest  neighbors  and  therefore  eight  radiating  bonds.  However, 
for  interior  sites,  each  of  these  eight  bonds  is  shared  by  one  other  lattice  site 
and  thus  each  interior  site  has  only  four  bonds  on  the  average. 

In  accordance  with  the  prescription  laid  out  above,  we  therefore  write  the 
following  expression. 

4t>B|aB|2  ^  +  "+  0(N2)=  I ,  (15a) 

where  N3  is  the  number  of  sites  in  the  lattice  and  N2  is  proportional  to  the 
number  of  surface  sites.  All  terms  on  the  left  of  eq.  (15a),  except  the  first,  are 
of  order  N2  or  less;  and  since  the  size  of  the  network  can  be  thought  of  as 
being  arbitrarily  large,  only  the  cubic  term  needs  to  be  considered.  Con¬ 
sequently,  the  following  expression  is  obtained: 

|aB|2  =  1/4Dd/V3.  (1 5b) 

It  turns  out  that  for  reasonable  values  of  the  barrier  height,  K„,  surface 
states  on  the  (0,  0,  I )  plane  have  energies  above  the  Fermi  level  and  are 
consequently  unoccupied.  As  a  result,  surface  states  do  not  influence  the 
surface  double  layer  potential  or  the  work  function,  a'  least  not  within  the 
requirements  of  the  network  model.  Thus,  the  properties  of  the  surface  states 
need  not  be  included  in  the  calculations  for  the  double  layer  and  work 
function  which  follow.  We  can  therefore  make  full  use  of  the  density  of 
states  derivable  by  the  method  of  Jelitto7). 

This  density  of  states  in  reduced  momentum  space  is: 

dIV 

p{kR)  =  Cn$\n{kR)  B,  (16) 

a  L 

where  d/VB/d£  has  been  expressed  rigorously  by  Jelitto  in  terms  of  elliptic 
integrals.  It  may  also  be  approximated  using  more  elementary  functions 
as  follows: 


d£B*2[l  ~  |£|]*[ln(“|r)]  C16  679  +  3.6364  |E|  +  2.4880 1£|2]  ; 

an  expression  good  to  four  significant  figures  over  the  range  0.005  <|£|<  1.0; 
E = cos  (kR)  and  CN  is  a  normalizing  constant  such  that  the  total  number  of 
quantities  in  the  region  O^kR^n  is  N3,  i.e., 

n 

J  p(kR)dkR  =  N3. 

0 
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Because  the  expression  for  p(kR)  is  rather  complicated,  the  above  integral  is 
evaluated  most  conveniently  by  numerical  means.  Unfortunately,  the  value 
of  djVB/d£  is  infinite  for  |£|=0,  and  the  numerical  integration  cannot  be 
performed  over  the  entire  integral  even  though  the  integral  converges.  The 
density  of  states  function,  p(kR),  is  symmetric  on  the  integral  about  the 
point  kR  =  {n,  and  consequently  the  numerical  integrations  need  only  be 
performed  from  k R  =  0  to  k  R  =  u,  where  u  < ^ir.  There  will  be  some  uncertain¬ 
ty  about  the  exact  value  of  CN,  depending  on  the  choice  of  upper  limit  for  the 
integral.  For  the  calculations  in  this  paper,  u  has  been  set  equal  to  0.4985  it. 
This  particular  value  was  chosen  for  two  reasons:  (1)  the  approximation  for 
diVu/d£  is  known  to  be  reasonably  accurate  over  this  range  of  integration, 
and  (2)  the  value  of  CN  differs  by  only  a  few  percent  from  the  value  assumed 
when  the  upper  limit  on  the  integral  is  changed  from  0.485  ji  to  0.49999  n. 


3.  Estimation  of  surface  charge 

The  charge  on  each  dangling  surface  bond  may  be  calculated  in  a  straight¬ 
forward  manner,  being  proportional  to  £)„•  Multiplying  this  ratio  by  2  to 
account  for  electron  spin  and  finally  integrating  over  all  states  (all  possible 
values  of  momentum)  gives  the  fraction  of  the  total  number  of  electrons  re¬ 
siding  on  a  surface  bond.  We  call  this  quantity  G(oiR): 

in  ^  COS^b) 

<?(«*)“  f  2p(kR)  Ds  d  (kR)  =  f -[M^ 2jV/  d  (A :R). 

J  J  ,  ,  sin  [2  (A/?  +  <5„)]  —  sin  (2<5B) 


,  +  sin  [2  (A-K +  <?„)] 
2kR 


sin  (2<5B) 


In  this  equation,  Ds  is  defined  by  eq.  (14a)  and: 


tan  (<5b)  = 


cos(AK)  - 
sin  (kR) 


I'/'s  v^)|2  =  lflB|2  cos2  (<5b)  . 


These  two  relations  are  derived  from  the  requirement  that  the  square  of  the 
modulus  of  the  standing  waves  represented  by  eqs.  (2a)  and  (2b)  must  be  the 
same  at  all  lattice  sites  for  wave  functions  associated  with  the  bulk  energy 
states. 

Because  of  the  difficulties  already  discussed  in  connection  with  p(kR), 
the  actual  upper  limit  of  integration  used  in  these  calculations  was  0.4985  ji 
instead  of  ]ir.  The  bulk  wave  functions  are  assumed  to  differ  only  by  a  phase 
factor  from  one  lattice  site  to  another  even  at  the  surface.  Because  G(aR)  is 
normally  a  few  tenths  of  an  electron  or  less,  and  because  the  conductivity  of 
the  metal  is  high,  all  the  positive  charge,  balancing  the  negative  charge 
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located  outside  the  lattice,  is  assumed  to  reside  on  the  crystal  surface.  The 
total  positive  surface  charge  is  therefore : 

ffs  =  w5fi(aR)^^(?,  (18) 

where  ns  is  the  number  of  dangling  surface  bonds  and  ( NJA  J  is  the  number 
of  sites  per  unit  area,  while  e  is  the  charge  on  an  electron. 


4.  Calculation  of  the  double  layer  potential 


The  work  required  to  move  an  electron  front  the  surface  to  some  point 
outside  the  crystal  may  be  calculated  using  Gauss'  theorem.  Since  the 
distribution  of  electric  charge  is  uniform  parallel  to  the  crystal  surface,  the 
field  acting  on  any  electron  is  perpendicular  to  the  surface  and  has  an  inten¬ 
sity: 

F(x)  =  ff(.v)/e0,  (19) 

where  e0  =  8.85x  10“12  C2/(Nm)2  and  ff(.v)  is  the  net  amount  of  positive 
charge  lying  between  the  surface  and  the  point  ,v  on  the  dangling  bonds. 
The  charge  <r(.v)  is  evaluated  from  the  expression: 


cos2(dn)vXp(-  2/?.v) 

_  [W-A/G2]r~ 

,  +  sin  [2JA-R  +  <>„)]  -  sin (2<>„) 
2k  R  ~ 


/>(*«) 

2  N1 


d  (kR) 


(20) 


in  which  we  have  combined  eq.  (17)  with  (2a).  The  work  performed  against 
the  double  layer  is  then: 


W'd=K0J  F(/)d/.  (21) 

0 

In  this  equation,  d /  represents  distance  along  the  direction  of  the  surface 
bond,  and  Ka  is  a  factor  such  that  ATt;d/  is  the  distance  normal  to  the  surface. 
The  integral  in  eq,  (21)  is  most  conveniently  evaluated  in  terms  of  the  reduced 
distance  ///?,  and  the  exact  value  of  the  upper  limit  is  not  very  important 
because  essentially  all  of  the  contribution  to  the  double  layer  potential  occurs 
within  less  than  five  times  the  nearest  neighbor  spacing  from  the  surface. 
Thus,  within  the  limitations  of  this  model,  changes  in  the  double  layer 
potential  from  one  surface  to  another  are  determined  by  the  variations  in 
three  factors:  (1)  the  number  of  lattice  sites  per  unit  area  of  surface,  (2)  the 
number  of  dangling  surface  bonds  at  each  site,  and  (3)  the  geometric  factor 
*0. 
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S.  The  work  function 

As  indicated  earlier,  the  work  function  is  defined  as  work  necessary  to 
move  an  electron  from  the  point  just  inside  the  lattice  to  a  point  just  outside. 
The  energy  required  for  this  process  depends  on  the  geometry  of  the  surface 
through  which  the  electron  passes  because  of  the  contribution  of  the  double 
layer.  In  keeping  with  the  simplicity  of  the  network  model,  the  work  function 
is  defined  as: 

=  vb  ~  Rr  ~  K  :  (22) 

whereas  before  VB  is  the  surface  barrier,  /iF  is  the  Fermi  level,  and  now  we 
denote  the  double  layer  potential  by  )FD.  The  first  two  quantities  on  the  right 
of  eq.  (22)  are  independent  of  surface  orientation,  and  the  variation  of  work 
functions  with  surface  is  therefore  due  to  I VD. 

The  double  layer  potential,  eq.  (21),  can  be  expressed  most  conveniently 
in  the  following  form: 

W0(otR)  =  ALl(aR),  (23) 

where 


and  is  only  a  function  of  the  nearest  neighbor  distance  and  surface  geometry, 
while  /(«/?)  is  an  integral  which  depends  only  on  the  product  txR.  The  values 
of  I(cxR)  for  various  <xR  have  been  listed  in  table  I .  In  table  2,  the  densities  of 


Table  I 


txR 

H<xR) 

aR 

/(<!/?) 

0.6 

0.10553 

1.2 

0.01447 

0.7 

0.05890 

1.3 

0.01207 

0.8 

0.03944 

1.4 

0.01025 

0.9 

002879 

1.5 

0.00882 

1.0 

0.02215 

1.6 

0.00767 

l.l 

0.01765 

Element 


Lithium 

Sodium 

Potassium 

Rubidium 

Cesium 


Table  2 


Density 

(g/cm3) 


0.534 

1.007 

0.870 

1.532 

1.873 


Nearest  neighbor 

A  i.(001) 

distance 

Ah  (Oil) 

(A) 

(V) 

(V) 

3.038 

25.82 

18.26 

3.666 

21  39 

15.12 

4.594 

17.08 

12.07 

4.937 

15.88 

11.23 

5.349 

14.67 

10.37 

A-ll 


. *■-.  ■  ■■  -- 


592 


I’.  K.  RAWLINGS  AND  H.  RIISS 


the  alkali  metals  and  their  first  nearest  neighbor  distanees  are  given,  as  well  as 
the  values  of  AL  for  the  (001)  and  (01 1)  planes. 

The  surface  barrier  Vh  may  be  estimated  by  the  following  expression: 

1  b  ~  /  +  II  +  /<F.  (24) 

Here  /  is  the  ionization  potential  for  an  atom  in  gas  phase,  //  is  the  heat  of 
sublimation  per  atom,  and  /iF  the  Fermi  level  which  can  be  estimated  from 
eq.  (8).  Eq.  (24)  approximates  an  equation  given  by  Wigner  and  Bardeen6) 
and  has  its  origin  in  comparing  certain  electronic  integrals  which  occur  in  the 
expression  lor  1  H  with  integrals  which  also  occur  in  expressions  for  /  and  //. 
The  barrier  height  may  also  be  estimated  by  adjusting  J'B  so  that  the  value 
predicted  by  the  network  theory  agrees  with  experiment.  This  was  done  Tor 
the  alkali  metals  using  the  work  function  calculated  for  the  (001)  plane. 
The  experimental  heats  of  vaporization,  ionization  potentials,  and  work 
functions  have  been  compiled  from  various  sources  by  Wigner  and  Bardeen6). 
In  table  3,  we  compare  the  barrier  heights  estimated  from  eq.  (24)  with  those 
obtained  through  adjustment  of  F„  to  make  the  work  function  calculated  by 
the  network  model  agree  with  experiment.  The  best  fit  value  obtained  from 
the  network  model  is  always  about  one-half  the  value  derived  from  eq.  (24). 
On  the  other  hand,  eq.  (24)  is  itself  by  no  means  precise.  Considering  the 
s  mplicity  of  the  network  model,  the  results  arc  not  very  disparate.  In  Tact, 
in  the  case  of  both  Hiickcl  theory  and  the  FEMO  method  applied  to  mole¬ 
cules,  the  same  order  of  magnitude  is  observed  for  the  mismatch  between 
theory  and  experiment. 

Also  included  in  table  3  are  double  layer  potentials  for  two  different  sur¬ 
faces.  Bardeen8),  using  the  free  electron  model,  calculates  that  the  surface 
double  layer  for  the  (110)  plane  of  sodium  is  0.4  V.  The  network  model 
predicts  a  value  of  0.37  V  (using  the  value  of  Vn  which  brings  the  theoretical 
and  experimental  work  functions  into  register)  as  compared  with  an  experi¬ 
mental  value8)  of  0.15  V. 

From  these  results,  it  appears  as  though  the  simple  network  model  posses¬ 
ses  about  the  same  quantitative  degree  of  validity  with  respect  to  surface 
state  problems  that  the  F-EMO  model  possesses  with  respect  to  molecules. 
This  suggests  that  it  may  be  used  to  make  ballpark  estimates  in  a  variety  of 
solid  state  situations. 
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Anisotropy  in  electrical  properties  of  fOOlj  Si/{0112j  Al203  f 

A.  J.  Hughes  and  A.  C.  Thorsen 

iC^T^^rn,,,on- 3J70M,ra,oma  Awnue- Anaheim ■ caii/ornia  mo3 


moiTSmm^Ai  n"  m"*  ?"  has  betn  camed  out  011  a  type 

°?l  S./|0II2|  AIjO,  films.  A  specially  designed  Hall  bridge  pattern  has  been  used  to  obtain  ^ 

lnhd;r,»ent!  me“Uremen,S  0f  mob‘ll‘*  as  a  function  of  current  direction  ,n  the  plane  of  the  film 
The  data  show  an  anisotropy  in  the  mobility  of  approximately  9%.  with  a  maximum  in  mobility 
(Kcurmg  along  the  <KX)>  St  direction  that  is  parallel  to  the  (2110)  Al203  direction  in  the  plane 

1  ani«1^  a'!h  ",  °r  “  f°Und  10  beaCOnSe<luence' lbrou^  Pie^oreststance  effect. of 
he  amsotropic  thermal  contraction  of  AJjO,  on  cooling  from  the  deposition  temperature,  which 
leads  to  an  anisotropic  thermally  induced  stress  in  the  Si 


INTRODUCTION 

The  Haii  mobility  in  epitaxiai  semiconductor  films  is 
usually  considered  to  be  a  good  criterion  of  film  quality 
and  has  been  used  extensively  in  the  optimization  of  the 
film  growth  process.  It  is  therefore  important  that  the 
properties  of  the  film  be  homogeneous  and  uniform  over 
the  piane  of  the  film  if  a  realistic  electrical  parameter 
is  to  be  deduced  from  a  particular  measurement.  In 
r.iost  cases  il  is  tacitly  assumed  that  the  mobility  is 
isotropic  in  the  piane  of  the  film  and  does  not  depend 
upon  the  direction  in  which  the  measui  'ment  is  made. 

We  present  results  in  this  paper  of  a  detailed  investiga¬ 
tion  of  the  electrical  properties  of  n-type  {001 }  Si/ 

{0ll"2}  Ai2Oa  and  show  that  for  this  system  an  assumption 
of  isotropy  is  not  justified. 1  The  mobility  is  found  to 
vary  with  current  direction  in  the  plane  of  the  film  and 
has  a  maximum  value  along  the  particular  (100)  Si  direc¬ 
tion  that  is  parallel  to  the  (2H0)  Ai203  direction  in  the 
plane  of  the  substrate. 

The  experimental  anisotropy  in  mobility  reported  in  this 
paper  wiii  be  explained  theoreticaiiy  in  terms  of  the 
piezoresistance  effect  and  anisotropic  substrate-induced 
thermal  stresses.  While  we  believe  this  to  be  the  prin¬ 
cipal  mechanism  behind  the  observed  mobility  aniso¬ 
tropy,  phenomena  associated  with  surface  eiectric  fieids 
and  surface  quantization  can  often  play  a  roie  in  current 
conduction  in  some  other  measurement  or  device 
situations. 

Surface  quantization  and  surface  transport  in  semicon¬ 
ductor  inversion  and  accumulation  layers  are  reviewed 
briefly  by  Stern2  and  a  number  of  recent  references  are 
listed.  Sato  el  al . 3  have  investigated  the  mobility  aniso¬ 
tropy  of  carriers  in  p- type  and  in  n-type  inversion  lay¬ 
ers  on  oxidized  Si  surfaces  in  MOSFET  device  configu¬ 
rations.  These  references  compare  field-effect  mobili¬ 
ties  in  each  of  two  perpendicular  orientations  for  a  num¬ 
ber  of  different  surface  piane  orientations.  For  some 
orientations,  sizabie  anisotropies  were  reported.  How¬ 
ever,  for  the  (001)  Si  orientation,  the  field-effect  mo¬ 
bility  is  isotropic  because  of  symmetry  considerations. 

The  work  of  Sato  et  al.  relates  specifically  to  bulk  Si 
and  therefore  did  not  consider  any  effects  of  substrate- 
induced  stresses  on  the  mobility.  If  there  were  a  surface 
perturbation  in  the  electrical  conductivity  of  heteroepi- 
taxiai  (001)  Si/AlaO,  films,  induced  by  surface  charge 
and/or  oxide  layer  impurities  and  taking  the  form  of  an 
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accumulation  or  inversion  layer,  the  resulting  surface 
component  of  mobility  would  also  be  isotropic,  except 
for  the  fact  that  there  is  an  anisotropy  in  the  film 
stresses.  The  mobility  cl'  a  surface  layer  would  there¬ 
fore  exhibit  a  small  anisotropic  effect  which  would  be 
similar  in  relative  magnitude  to  that  found  in  bulk  Si 
under  similarly  stressed  conditions.  It  is  unlikely  there¬ 
fore,  that  a  small  component  of  surface  conductivity 
could  contribute  substantially  to  the  total  mobility  aniso¬ 
tropy  that  would  be  measured  in  (001 )  Si/ALjO,  films. 

The  mobility  anisotropy  to  be  discussed  in  the  present 
paper  is  believed  not  to  be  associated  with  either  sur¬ 
face  quantization  or  self-accumuiation  or  Inversion  ef¬ 
fects.  There  is  no  applied  electric  field  normal  to  the 
surface  and  the  oxide  film  formed  during  the  after¬ 
growth  annealing  sequence  is  removed  prior  to  the  mea¬ 
surement  of  film  properties.  In  addition,  the  high  donor 
concentrations  in  the  films  tend  to  mask  contributions  to 
the  electrical  conduction  from  surface  effects.  To  nu¬ 
merically  assess  the  magnitude  of  any  possible  surface 
effects  present  under  these  conditions,  however,  would 
require  the  measurement  of  eiectrical  properties  in  a 
MOS  device  configuration,  and  this  has  not  been  done  in 
this  work. 

Coilectively,  the  above  arguments  form  the  basis  for 
our  assumption  that  the  mobility  anisotropy  reported 
here  is  a  bulk  effect  and  is  due  primarily  to  substrate- 
induced  thermal  stresses  acting  through  the  piezoresis¬ 
tance  effect.  The  agreement  between  theory  and  experi¬ 
ment  lends  confidence  to  this  interpretation. 

THEORY 

In  order  to  explain  the  experimental  observations  dis¬ 
cussed  below,  we  have  examined  theoreticaiiy  the  effect 
of  stress  in  the  Si  film  on  the  electron  transport  prop¬ 
erties  of  the  film.  As  is  well  known,  transport  proper¬ 
ties  such  as  resistivity  and  mobility  can  be  reiated  to 
applied  stresses  through  the  piezoresistance  effect.  Nu- 
meronn  sfudUs  cf  tl'AJ  pitioftsista.iCe  effect  in  bulk  Sl 
have  been  reported  in  the  literature  over  a  number  of 
years.  Much  less  attention  has  been  paid  to  stress  and 
its  effect  on  the  transport  properties  of  epitaxial  films. 

Stress  in  Si  films  on  sapphire  (A1203)  has  been  examined 
by  Dumin,4  and  in  Si  films  on  spinel  substrates  both  by 
Schldtterer5  and  by  Robinson  and  Dumin. 6  Schldtterer 
presents  formulae  for  the  fractional  change  in  resistiv- 
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of  the  current  flow.  Then  £',A  =  p0(l  +  nltTk),  where  p„ 
is  the  zero-stress  Si  resistivity  and  k  is  summed  from 
1  to  6.  Defining  Ex/Jx  as  p,,  the  iongitudinai  piezoresis¬ 
tance  is  then  given  by 

^P\  f Pi)  fpi  ™  Po^/ Pet  =  [E[  /Jx pfx)  —  \  “  It x (1) 

The  iongitudinai  piezoresistance  effect  can  thus  be  cal¬ 
culated  if  the  rr  coefficients  and  the  Si  fiim  stresses  Tk 
are  known.  For  a  Cartesian  coordinate  system  with 
axes  aiong  the  Si  crystai  axes,  the  rru  are  particularly 
simple  and  the  sum  involves  oniy  rrn7i  +  nxAT,  +  F,).  For 
an  arbitrary  orientation  of  axes,  which  is  of  primary  in¬ 
terest  to  us  here,  "rotated”  n  coefficients  are  given  by 
Pfann  and  Thurston.  ’  These  transformed  ir  coefficients 
aro  given  in  terms  of  the  direction  cosines  of  the  arbi¬ 
trary  coordinate  system  with  respect  to  the  Si  crystal 
axes.  The  determination  of  ir  coefficients  appropriate  to 
an  arbitrary  coordinate  system  is  straightforward  but  is 
laborious. 


FIG.  1.  Specially  designed  Hall  bridge  pattern  permitting  an 
independent  measurement  of  electrical  properties  every  18  deg 
In  the  plane  of  the  film. 


ity  (Ap/po)  in  terms  of  piezoresistance  coefficients  rrn, 
ir12,  rrM  and  an  assumed  isotropic  stress  a,  in  the  Si  fiim 
due  to  differential  thermal  contraction  between  film  and 
substrate.  Since  he  was  dealing  with  spinei  substrates, 
it  was  assumed  that  the  thermal  stress  and  the  induced 
changes  in  resistivity  were  isotropic  in  the  piane  of  the 
film.  In  this  section  we  treat  the  anisotropy  in  stress 
and  piezoresistivity  and  present  formulae  which  are  .ap¬ 
plicable  to  the  anisotropic  case. 


We  first  consider  the  so-calied  longitudinal  piezoresis¬ 
tance  effect  for  the  case  in  which  there  is  only  one  com¬ 
ponent  of  current  («/,),  and  the  field  Ex  is  in  the  direction 


We  next  consider  the  determination  of  stresses  in  the  Si 
film  for  the  case  of  (001)  Si  on  (0lT2)  Ai203.  A  state  of 
compressive  stress  arises  in  the  Si  film  due  to  the  rei- 
ative  thermal  contraction  of  the  Si  and  the  A1203  sub¬ 
strate.  In  order  to  calculate  this  stress,  we  consider 
the  case  of  a  thin  fiim  on  a  relatively  (say  100—150 
times)  thicker  substrate  so  that  bending  of  the  composite 
can  be  neglected,  and  assume  that  the  strain  induced  in 
the  Si  is  proportional  to  the  difference  in  thermal  con¬ 
traction  between  the  two  materiais.  In  addition,  the 
thermai  contraction  is  assumed  to  take  piace  over  a 


Before  presenting  the  theoretical  results  obtained,  wo 
will  briefly  sketch  the  method.  A  more  detailed  theoret¬ 
ical  discussion  will  be  given  elsewhere.1  The  theoreti¬ 
cal  determination  of  the  effect  of  thermally  induced 
stresses  on  film  resistivity  and  mobility  can  be  divided 
into  three  parts:  (a)  calculation  of  the  piezoresistance 
effect  in  terms  of  stresses;  (b)  calculation  of  these 
stresses  in  terms  of  the  various  coefficients  of  thermai 
expansion;  and  (c)  successive  transformations  of  fiim 
and  substrate  coordinate  systems  to  achieve  the  proper 
relative  crystallographic  orientation  and  to  include  an¬ 
isotropy  in  the  plane  of  the  Si  fiim. 


The  resistivity  and  stress  are  related  through  the  equa¬ 
tion8  El=pijJl  +  rtijkl  Jt  Ttt ,  where  E  is  the  electric  fieid, 

/  the  current  density,  p  the  zero-stress  resistivity  ten¬ 
sor,  F  the  stress  tensor,  and  rr  the  piezoresistance  ten¬ 
sor.  For  cubic  crystals  such  as  Si,  pu  is  both  diagonal 
and  isotropic  and  equal  to  p0 5(;,  where  p0  is  the  zero- 
stress  resistivity.  This  tensor  equation  is  commonly 
contracted  to  a  single-subscript  notation  for  p  and  T  and 
to  a  double-subscript  notation  for  rr. 8  In  this  notation,  T 
can  be  written  as  a  6-component  column  vector  and  ir  as 
a  6x6  matrix.  Referred  to  the  Si  crystal  axes,  ir  has 
only  three  independent  coefficients:  ir„,  rrI2,  and  rrM; 
however,  ali  36  rr  coefficients  may  be  nonzero  when  re¬ 
ferred  to  a  Cartesian  coordinate  system  having  an  aibi- 
irary  orientation  with  respect  to  the  Si  crystal  axes. 


Fit..  2.  Typical  room-temperature  mobility  data  for  two  1001) 
Si/(0ll2)  A1203  samples.  The  solid  curves  are  least-squares 
theoretical  fits  to  the  experimental  points.  The  maximum 
mobility  Is  at  9  =  0”  and  is  along  the  [tool  Si  direction  and  the 
minimum  is  at  9  =  90”  along  the  [OlOl  Si  direction. 
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61  Degrees  I 


FIG.  3.  Low-temperature  (77  K)  mobility  data  for  one  (001)  St/ 
(0lT2)  AI2O3  sample  and  corresponding  least-squares  theoret¬ 
ical  eurve. 


the  substrate  and  required  for  calculation  of  the  piezo¬ 
resistance  effect  have  now  been  obtained. 

We  next  return  to  the  longitudinal  piezoresistance  effect 
and  to  Eq.  (1).  The  stresses  calculated  above  relate  to 
the  coordinate  system  and  parallel  Si/Al20,  relations 
given.  Using  the  results  of  Ref.  9,  the  tru  coefficients 
in  Eq.  (1 )  are  thon  transformed  to  the  same  coordinate 
systems  used  for  the  stress.  The  piezoresistance  (ApJ/ 
Po)  thus  obtained  represents  the  change  in  resistivity 
due  to  stress  for  a  current  Jt  and  field  £,  in  the  x  direc¬ 
tion  11  [100]  Si.  We  are,  however,  interested  in  the  pi¬ 
ezoresistance  as  a  function  of  angle  in  the  plane  of  the 
film.  A  second  transformation  on  both  the  tr’s  and  the 
stresses  is  then  performed  which  represents  a  rotation 
about  the  z  axis.  The  angle  8  of  this  rotation  is  mea¬ 
sured  from  the  x  axis  and  is  positive  toward  the  y  axis. 
The  longitudinal  piezoresistance  effect  corresponding  to 
current  flow  at  an  angle  8  from  the  [100]  Si  direction 
then  beco.-nes 


temperature  range  from  room  temperature  to  the  growth 
temperature. 10 

The  Cartesian  coordinate  system  employed  is  defined  as 
follows.  The  z  axis  is  normal  to  the  surface  of  the  Si 
film  and  thus  is  perpendicular  to  (001 )  Si  and  also  per¬ 
pendicular  to  (0lT2)  A1203.  The  x  and  y  axes  are  in  the 
plane  of  the  film  with  x  11  [lOO]  Si  and  II  [2110]  A1203,  and 
y  I'  [010]  Si  and  II  [01 11]  ai2o3. 

The  six  stresses  T, ,  T2,  T,,  T4,T,,Te  (=  T  ,T  ,T  ,T  T 

1>  Zl  31  <1  51  8  V  *II>  *yy>  *yt» 

rxy,  respectively)  completely  define  the  stress  in  the  Si 
film.  The  stresses  T3,  TA,  and  7;  are  zero  at  the  free 
surface  of  the  Si  film  and  if  we  further  assume  that  the 
film  is  thin  enough  that  all  stresses  are  umtorm  in  the 
z  direction,  then  T3=Tt  =  T5  are  zero  everywhere  in  the 
film.  r6  is  in  general  nonzero,  but  for  the  relative  ori¬ 
entation  of  film  and  substrate  under  consideration  is 
identically  zero.  Thus  only  7]  and  T2  are  nonzero.  De¬ 
noting  the  elastic  constants  of  Si  by  CtJ  and  the  stiffness 
coefficients  by  S{/,  we  find  the  following  implicit  rela¬ 
tions  for  the  stresses  7]  and  T2: 

HT,  +  T2)=  [(Cj,  +  C„C12  -  2Cj2)/2CM][l.  287(af  -  a, ) 

+  0.713(as-a2)]A7' 

=  [l/2(S,,  +S12)][l.  287(a,  -  a,) 

+  0.713(«3-a2)]A7\  (2) 

and 

?(T,  -  T2)=  (-0.  713)|(Cn  -C’I2)(a,  -  «2)A T 

=  -0.  713[S11/2(S11  -Sl2)(Su  +  2S12)](«,  -  «2)A7\ 

(3) 

Here  «s  is  the  (isotropic)  Si  thermal  expansion  coeffi¬ 
cient  and  a,  and  «2  are  thermal  expansion  coefficients 
for  the  A1203  substrate,  perpendicular  to  the  c  axis  of 
A1203  and  parallel  to  the  c  axis  of  A1203,  respectively. 
A7’is  the  temperature  difference  between  growth  and 
room  temperatures  and  in  Eqs.  (2)  and  (3)  is  understood 
to  be  a  positive  number.  In  the  limit  »,  =  «2,  7]  and  T2 
are  equal,  and  Eq.  (2)  reduces  to  that  given  by 
Schldtterer. 5  The  thermal  stresses  in  the  Si  film  due  to 


APi/ Po~  Ex/Jxp0  -  1 

=  5(7;  +  T2)(jtu  +  n,2)  +  %(T,  -  T2)(nu  -  n,2)cos28,  (4) 

where  7]  and  T2  are  given  in  Eqs.  (2)  and  (3).  To  avoid 
ambiguity  in  orientation,  we  repeat  that  the  angle  8  is 
measured  from  that  (100)  Si  direction  in  the  plane  of  the 
film  that  is  parallel  to  the  (2110)  A1203  direction  in  the 
plane  of  the  substrate.  Note  that  the  term  depending  on 
8,  which  determines  the  amount  of  orientational  aniso- 
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Sample 

^n*l 

4a 

a  a> 

rms 

4o 

error 

(%> 

Room- 

-temperature  data 

102 

446 

408 

427 

8.9 

1.6 

615 

18-1 

432 

397 

415 

8.6 

1.4 

598 

15 

472 

425 

449 

10.5 

2.5 

648 

3412-2 

513 

472 

493 

8.2 

1.9 

711 

3411-1 

516 

478 

497 

7.  0 

2.  1 

717 

18-2 

422 

387 

405 

8.6 

2.5 

584 

31 

451 

401 

426 

11  6 

2.2 

614 

13 

443 

394 

419 

11  7 

2.9 

604 

(aver¬ 

age) 

(462) 

(420) 

(441) 

(9,  5) 

(636) 

Liquid- 

-N,  temperature 

(77  K)  data 

18-1 

1014 

797 

906 

21 

5,  7 

1305 

31 

1063 

832 

948 

24 

7.0 

1367 

tropy  in  resistivity  and  mobility,  also  depends  directly 
on  nr,  —  aa.  That  is,  the  anisotropy  in  piezoresistance  in 
the  (00!)  Si  plane  depends  directly  on  the  anisotropy  in 
thermal  expansion  coefficients  of  the  A1,03  substrate. 
Using  o')  =  3.9xlO'6/°C,5  a,  =  8.  31  x  10-8/°C,  a2  =  9.03 
xlO"e/”C,u  A7’=1100"C,  and  the  known  elastic  con¬ 
stants  of  Si,13  Eqs.  (2)  and  (3),  yield  stresses  of 

j(7]+  r2)=-0.9206xl0,odyn/cm>,  (5) 

|(r,  -  T2)  =  +  0.  2852  x  109  dyn/cmJ. 

Substituting  values  of  piezoresistance  coefficients  for  n- 
type  Si  from  Smith, 13•,,  Eq.  (4)  becomes 

Ap,/po  =  0.  44192-0.  04449  cos2fl.  (6) 

Since  the  mobility  p  can  be  related  to  the  resistivity  p, 
by  P  =  RH/plt  where  Rw  =  the  Haii  constant,  the  theoret¬ 
ical  anisotropy  in  mobility  can  be  written 


4/(i0  =  U  +  (Aft/po)]*1  =  (1. 44192  -  0.  04449  cos20)-', 


(V 


where  p0  is  the  zero-stress  mobility.  From  Eq.  (7),  we 
note  that  the  mobility  will  be  a  maximum  along  the  [100] 
Si  direction  and  wiil  be  a  minimum  90°  away  along  the 
[010]  Si  direction. 


The  amount  of  anisotropy  can  be  described  conveniently 
by  a  parameter  A  which  we  define  as  A  =  (u  -  u  .  )/ 
4*,  where  pA  =  2(pmal  f  pml0)  is  an  approximate  average 
mobility  in  the  plane. 15  The  predicted  mobiiity  aniso¬ 
tropy  is  found  to  be  A  =  6.  2 %  and  the  average  mobiiity 
is  found  to  be  PA  =  0.  694  p0.  Thus  we  find  a  6.  2%  aniso¬ 
tropy  in  mobiiity  superimposed  on  a  substantiai  30% 
over -all  theoreticai  reduction  in  mobiiity  for  n-type  Si, 

The  theory  aiso  predicts  a  transverse  piezoresistance 
effect  corresponding  to  an  electric  fieid  E2  in  the  plane 
of  the  film  and  orthogonal  to  the  current  J,  direction. 
This  transverse  piezoresistance  effect  is  also  found  to 
depend  upon  orientation  and  is  given  by 

EzMpo  =  “5(^1  -  T3)(i tu  -  ,T12)sin2fl=  +  0. 04449  sin20  , 

(8) 
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and  is  zero  along  the  [100]  SI  direction  and  is  of  maxi¬ 
mum  magnitude  along  the  [010]  Si  direction.  The  sign 
convention  employed  here  is  that  E,  is  positive  if  J,  xEa 
is  out  of  the  plane  of  the  film  and  is  thus  in  the  +  z 
direction. 


The  transverse  piezoresistance  (E,/^)  associated  with 


a  field  E3  normal  to  the  Si  film  was  also  examined  and 


for  (001)  Si  on  (01l2)  Al2Os,  is  found  to  be  identically 
zero. 


Although  the  results  described  above  were  obtained  ex¬ 
plicitly  for  n-type  Si,  the  formalism  Is  also  valid  for 
p-type  Si  providing  that  appropriate  vaiues  for  the  n  co¬ 
efficients  are  used.  The  quantity  jrM  -ir12,  which  is  im¬ 
portant  for  the  anisotropic  effect,  is  found  to  be  approx¬ 
imately  a  factor  of  20  less  for  p-type  Si.  As  a  result, 
little  anisotropy  would  be  expected  in  the  Si  orientation 
under  discussion,  although  other  orientations  can  in 
general  show  pronounced  effects. 


EXPERIMENTAL 

The  films  used  for  this  study  were  grown  by  chemical 
vapor  deposition  techniques  in  two  different  reactor 
systems.  One  of  these  Is  a  vertical  system16  and  the 
other  Is  a  horizontal  system.  In  both  cases  the  films 
were  formed  from  the  thermal  decomposition  of  silane 
(SiHj,  with  H2  used  as  a  carrier  gas.  Intentional  doping 
to  concentrations  of  (1— 6)xl016  cm"3  was  achieved  with 
the  use  of  arsine  (AsH3).  Growth  temperatures  ranged 
from  955  *C  in  the  horizontal  system  to  1075 "C  in  the 
vertical  system,  and  film  thicknesses  varied  from  1  5 
to  1.8  p. 

After  growth  the  samples  were  annealed  at  1100  “C  in  02 
for  30  min  followed  by  a  N2  anneal  for  2  h  at  1100  "C  in 
order  to  stabilize  film  properties  and  electrically  neu¬ 
tralize  any  A1  impurities  in  the  film. 17  The  resulting 
oxide  film  was  removed  prior  to  making  electrical  mea¬ 
surements. 


Measurements  of  resistivity  and  Hall  effect  were  made 
on  each  film  using  a  speclaiiy  designed  Hall  bridge  pat¬ 
tern  etched  in  the  epitaxial  layers.  This  pattern,  shown 
in  Fig.  1 ,  consists  of  two  wheel-shaped  bridges.  The 
five  arms  of  a  singie  bridge  are  separated  by  72  deg 
from  each  other,  and  the  arms  of  one  bridge  are  rotated 
by  18  deg  with  respect  to  those  of  the  other  bridge.  This 
allows  an  independent  measurement  of  electrical  prop¬ 
erties  (Haii  coefficient  and  resistivity)  every  18  deg  in 
the  plane  of  the  film. 


Electrical  data  were  taken  on  eight  samples  of  {001}  Si/ 
{01T2}  A1203.  For  every  sample,  the  orientation  of  each 
leg  of  the  Hall  pattern  was  determined  with  respect  to 
one  of  the  Si  directions  of  the  form  (100)  in  the  plane  of 
the  film.  This  was  accomplished  by  obtaining  a  Laiie 
back-reflection  x-ray  pattern  for  each  substrate  and 
from  this  locating  a  (2lf0)  Al203  direction  in  the  piane 
of  the  substrate  that  is  parallel  to  a  (100)  Si  direction  in 
the  film  plane.  The  mobility  in  each  leg  can  then  be 
plotted  versus  the  angle  8  between  the  current  direction 
and  that  particular  (100)  Si  axis.  Typical  room-temper¬ 
ature  mobility  data  from  two  samples  are  shown  plotted 
In  Fig.  2.  The  mobility  can  be  seen  to  vary  with  angle 
and  appears  to  approach  a  maximum  at  6  =  0°  and  a  mln- 


b-4 


a 

j 


2308  A.J.  Hughes  and  A.C.  Thor  sen:  Anisotropy  in  properties  of  (001)  Si/(01l2)  Al,03 


imum  at  5*90°,  This  Is  consistent  with  the  behavior 
predicted  by  Eq.  (7)  and  suggests  that  the  theoretical 
form  of  the  mobility  could  be  used  to  fit  the  experimen¬ 
tal  results. 

The  theory  leading  to  Eq,  (7)  predicts  that  the  mobility 
anisotropy  should  be  of  the  form 

m/m0=  -6cos28)‘‘  ,  (g) 

where  a  and  b  are  constants.  Equation  (9)  cannot  be  em¬ 
ployed  directly  in  analyzing  the  experimental  results  in 
that  the  zero-stress  mobility  p0  is  not  known  and  would 
be  expected  to  be  a  function  of  growth  conditions  and  to 
vary  slightly  from  sample  to  sample  or  run  to  run.  Ac¬ 
cordingly,  we  have  fitted  the  experimental  results  by  the 
method  of  least  squares  to  a  theoretical  curve  of  the 
form 

M=  («' -/>'  cos20)-'.  (jO) 

The  numerical  curves  thus  determined  are  shown  by  the 
two  solid  curves  in  Fig.  2  for  two  samples.  The  aniso¬ 
tropy  parameter  A  is  independent  of  m  and  is  equal  to 
2b/a  =  2b'/a'. 

The  electrical  data  taken  on  the  eight  samples  of 
{001}  Si/{01 1  T"2}  Al2Oa  were  all  fitted  by  least  squares 
and  analyzed  in  terms  of  Eq.  (10).  The  maximum  mobil- 
minimum  mobility  pmln,  average  mobility  p 
and  anisotropy  parameter  A  were  calculated.  These  four 
quantities  are  tabulated  in  Table  I  for  the  samples  mea¬ 
sured.  it  can  be  seen  that  the  anisotropy  A  varies  from 
approximately  7.  6^.  to  11. 7%  with  an  average  value  of 
9.  . 

The  scatter  in  the  mobility  data,  as  evidenced  by  the  re¬ 
sults  shown  in  Fig.  2,  is  probably  due  primarily  to  in¬ 
homogeneities  in  film  properties  over  the  surface  of  the 
film.  For  example,  carrier  concentrations  measured  on 
the  separate  areas  of  the  Hall  pattern  on  a  given  film  are 
found  to  vary  an  average  of  ±  7?.  Even  though  these  dif¬ 
ferences  are  taken  into  account  m  calculating  the  mobil¬ 
ity  in  each  arm  of  the  bridge,  slight  errors  may  be  in¬ 
troduced  since  the  spatial  extent  over  which  the  resistiv¬ 
ity  is  measured  is  much  larger  than  that  over  which  the 
Hall  constant  (carrier  concentration)  is  measured  (see 
configuration  of  bridge  in  Fig.  1). 

The  rms  error  ((% )  between  the  experimental  points  and 
the  fitted  theoretical  curve  for  each  sample  is  shown  in 
Table  I.  The  error  ia  all  cases  is  sufficiently  small 
compared  with  the  anisotropy  to  indicate  that  the  fit,  and 
hence  the  correlation  between  theory  and  experiment,  is 
statistically  significant. 

Measurements  were  taken  at  77  K  on  two  selected  sam¬ 
ples  and  these  results  are  also  tabulated  in  Table  I. 
Representative  data  for  one  sample  are  shown  in  Fig.  3, 
The  parameters  deduced  from  the  curve-fitting  proce¬ 
dure  for  this  data  are  also  shown  in  Table  I.  The  aniso¬ 
tropy  A  is  found  to  increase  in  going  from  room  temper¬ 
ature  to  77  K  by  roughly  a  factor  of  3  for  those  samples 
measured.  Correlation  of  this  increase  in  anisotropy 
with  theory  would  require  information  on  the  variation  of 
the  Si  piezoresistance  coefficients  nu  and  7ria  with  tem¬ 
perature.  Data  for  tt,,  as  a  function  of  temperature  are 
available  in  the  literature18  and  indicate  that  x,,  in¬ 
creases  bv  a  factor  close  to  3  in  frim  fuuto  tem¬ 


perature  to  77  K.  Data  on  7rla  as  a  function  of  tempera¬ 
ture  do  not  appear  to  be  available;  however,  it  is  not 
unreasonable  to  assume  that  the  temperature  depen¬ 
dence  of  both  coefficients  is  the  same.  If  this  is  the 
case,  then  the  experimental  increase  la  anisotropy  at 
low  temperature  is  consistent  with  theory. 

The  last  column  in  Table  I  lists  values  of  v0,  the  zero- 
stress  mobility.  As  mentioned  earlier,  /i0  cannot  di¬ 
rectly  be  determined  from  experimental  data  but  must 
be  obtained  from  a  combination  of  theory  and  experi¬ 
ment.  The  relationship  between  zero-stress  mobility  p0 
and  the  average  mobility  may  be  written  as  15 

^0=  a  UAa  .  (II) 

We  then  assume  that  Eq.  (7)  holds  and  that  a  =  1. 44192. 
The  values  of  p0  thus  determined  are  given  in  the  last 
column  of  Table  I  and  range  from  a  high  of  717  to  a  low 
of  584,  with  an  average  of  636. 

The  rather  low  value  of  p0  compared  with  bulk  mobilities 
in  Si  (>  1000  cm2, At  sec  at  these  carrier  concentrations) 
should  be  noted  and  strongly  suggests  that  a  consider¬ 
able  mobility  rtdurUcr.  results  Inthtt*  films  from 
causes  other  than  thermal  stress. 

We  next  consider  the  transverse  piezoresistance  effect 
described  in  Eq.  (8).  We  have  attempted  to  measure  the 
transverse  piezoresistance  voltage  appearing  across 
appropriate  terminals  of  the  Hall  bridge  fV'W.'f/T  tb* 
small  size  of  the  voltage  makes  this  determination  ex¬ 
tremely  difficult.  The  voltage  is  measured  across  the 
same  terminals  as  the  Hall  voltage,  with  no  magnetic 
field.  As  a  consequence,  any  misalignment  in  the  ter¬ 
minals  would  result  in  an  "offset”  voltage  which  would 
tend  to  mask  the  transverse  voltage  and  introduce 
errors.  In  addition,  inhomogeneities  in  bridge-arm 
thickness  and  width  also  introduce  some  ambiguity. 

Notwithstanding  the  above  difficulties,  some  experimen¬ 
tal  transverse-effect  data  have  been  obtained.  In  obtain¬ 
ing  these  data,  the  sense  of  the  transverse  voltage  must 
be  properly  accounted  for  since  Eq.  (8)  predicts  a  neg¬ 
ative  voltage  over  some  range  of  angles.  The  sign  con¬ 
vention  employed  is  that  the  transverse  field  is  taken 
positive  if  the  current  direction  crossed  into  the  trans¬ 
verse  field  yields  a  vector  out  of  the  plane  of  the  film. 

Representative  data  on  one  sample  taken  at  300  K  are 
shown  in  Fig.  4.  A  least-squares  fit  of  these  data  to  a 
sin20  term  as  suggested  by  Eq.  (8)  is  shown  by  the  solid 
curve  in  Fig.  4.  The  resulting  coefficient  of  the  sir;?, 6 
term  is  approximately  equal  to  0.08,  which  is  larger 
than  predicted  from  Eq.  (8);  however,  the  fit  of  the  ex¬ 
perimental  points  is  not  good.  The  value  of  p0  used  in 
Eq.  (8)  is  the  value  of  determined  from  the  longitudi¬ 
nal  mobility  measurements  described  above. 

The  transverse  voltage  also  becomes  substantially 
larger  at  77  K,  and  low-temperature  data  are  also 
plotted  in  Fig.  4.  The  magnitude  of  the  increase  is  a 
factor  of  between  2  and  3,  similar  to  the  increase  found 
for  the  longitudinal  effect.  Note  that  the  least-squares 
fit  to  a  theoretical  sin20  term  is  much  better  at  77  than 
a1,  3WJ  K,  Jut  probably  to  the  larger  signal. 
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SUMMARY  AND  CONCLUSIONS 


Detailed  studies  of  the  electrical  properties  of  {001}  Si/ 
{OUty  AljOj  have  shown  the  existence  of  a  significant 
anisotropy  in  mobility.  This  anisotropy  can  be  accounted 
for  through  the  piezoresistance  effect  in  terms  of  a  sim¬ 
ple  modei  of  thermally  Induced  stress  taking  Into  ac¬ 
count  the  difference  In  thermal  expansion  coefficients  of 
Alp,  parallel  and  perpendicular  to  the  c  axis.  Both 
longitudinal  and  transverse  piezoresistance  effects  were 
considered  and  theoretical  formulae  developed  which  ac¬ 
count  for  the  anisotropic  effects. 


The  theory  for  the  longitudinal  piezoresistance  effect 
predicts  an  anisotropy  in  mobility  of  about  6.  2%  and  for 
the  transverse  piezoresistance  effect  predicts  a  sin2fi 
anisotropy  coefficient  of  about  0,044  for  the  normalized 
transverse  electric  field.  The  value  of  each  of  these 
quantities  is  strongly  dependent  upon  the  data  used  for 
the  thermal  expansion  of  A1203.  Recent  measurements13 
at  this  laboratory  using  a  differential  technique  have 
given  a  difference  in  thermal  expansion  coefficients 
"i  -«a  of  (-1-  08±0. 12)xlO‘8/°C  for  the  At203  sub¬ 
strate.  Using  this  value  for  o,  -  ct2  leads  to  a  predicted 
mobility  anisotropy  of  9.  3% 20  and  a  transverse  electric 
field  anisotropy  factor  of  0.067  which  are  in  good  agree¬ 
ment  with  the  experimental  results. 


A  few  concluding  points  suggested  by  our  investigation  of 
the  electrical  properties  of  {001}  Si/{0ir2}  Ala03  should 
be  mentioned.  The  first  is  simply  that  the  anisotropy  in 
mobility  reported  here  must  be  taken  into  account  in  the 
evaluation  and/or  optimization  of  film-growth  pro¬ 
cesses.  Previous  practice  has  apparently  been  to  ignore 
the  orientation  of  current  flow  in  the  plane  of  the  film, 
and  thus  mobility  data  for  various  {001}  Si  films  could 
possess  a  built-in  *10%  scatter  due  to  this  anisotropy. 
Such  scatter  renders  the  task  of  definitively  evaluating 
changes  in  film-growth  processes  somewhat  difficult. 


Second,  for  some  Si  film  applications  in  which  the  mo¬ 
bility  is  important  it  would  be  desirable  to  orient  the 
current  flow  (along  the  (100)  SI  direction  II  <2110>  A1203 
direction)  to  obtain  the  maximum  mobility. 


The  third  and  final  point  concerns  the  interpretation  of 
the  origin  of  the  mechanisms  determining  mobility  in 
these  heteroepitaxial  films.  The  excellent  agreement 
between  theory  and  experiment  for  the  mobility  aniso¬ 
tropy  suggests  that  the  anisotropy  can  be  substantially 
accounted  for  in  terms  of  thermal  stresses  induced  by 


the  Al,Oj  substrate. 


In  principle,  residual  growth  stresses  in  the  Si  due  to 
lattice  mismatch  and  effects  of  dislocation  distributions 
in  the  films  could  also  lead  to  anisotropies  In  carrier 
mobility.  Our  results,  however,  Indicate  that  the  aniso¬ 
tropy  can  be  adequately  explained  without  recourse  to 
these  effects.  On  the  other  hand,  the  rather  low  average 
value  (*636  cm2/V sec)  deduced  for  the  zero-stress  mo¬ 
bility  p0  from  analysis  of  the  theoretical  and  experimen¬ 
tal  data  indicates  that  thermally  induced  stress  is  not 
the  dominant  mechanism  in  lowering  the  over-all  aver¬ 
age  film  mobility  from  bulk  Si  values.  Thus,  defect 
structures,  such  as  dislocations,  would  appear  to  tlav 
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an  Important  role  in  determining  heteroepitaxiai  Si  film 
mobilities. 


At  any  rate,  these  studies  of  mobility  anisotropy  have 
yielded  more  detailed  information  than  has  previously 
been  available  from  studies  of  mobility,  This  also  sug¬ 
gests  that  additonal  detailed  Investigations,  which  ex¬ 
amine  the  variation  of  anisotropy  with  film  thickness, 
growth  romflilcns,  arid  temperature,  may  well  piovt 
valuable  in  better  understanding  and  improving  Si/ALO, 
films. 
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Studies  of  the  resistivity  and  Hail  effect  in  n-tvoe  (221)  si/mrn  ai  n  rn  u 

“  <■*  «'“•'»»  "»“»*>  i*  • 

4(1?  rhVar  °ns  n  mobility  between  hieh‘  and  low-mobility  directions  are  approximately’ 

S’,  ,C  fT'  in  n,obility  ls  exl>laincd  on  the  basis  of  pie^oresistLce  eHeeU  in 

1  ut> sU- aie^ S S Lh  ‘ndUCCd  by  the  dlffcrcnUal  thermal  contraction  between  the  Si  and  the  Al203 


An  implicit  assumption  that  is  usually  made  in  the  elec¬ 
trical  characterization  of  silicon  heteroepitaxial  films 
on  sapphire  (A1203)  substrates  is  that  the  electrical  prop¬ 
erties  are  isotropic  in  the  plane  of  the  film.  Although 
this  assumption  is  approximately  correct  in  the  (100)  Si 
growth  plane,  1  we  report  results  in  this  letter  showing 
that  the  electrical  properties  are  extremely  anisotropic 
in  the  (221)  Si  plane.  It  is  shown  that  the  anisotropy  can 
be  explained  on  the  basis  of  the  piezoresistance  in  Si 
associated  with  stresses  induced  by  the  differential 
thermal  contraction  between  the  Si  and  the  A1203 
substrate. 


rate  of  approximately  2  pin/min.  After  growth  the 
samples  were  annealed  at  1100  C  in  02  for  30  min,  fol¬ 
lowed  by  a  N2  anneal  for  2  h  at  1100°C  in  order  to  stabi¬ 
lize  film  properties  and  electrically  neutralize  any  Al 
impurities  present  in  the  film.4  After  removal  of  the 
resulting  oxide  film,  electrical  measurements  were 
carried  out  on  a  specially  designed  double -Hall -bridge 
pattern  etched  in  the  Si  film  by  the  use  of  photolitho¬ 
graphic  techniques.  The  bridge  allows  an  independent 
measurement  of  Hall  mobility  (p)  every  18r  in  the  plane 
of  the  film. 


The  (221)  orientation  of  Si  is  of  particular  interest  since 
the  piezoresistance  effect  leads  to  the  prediction  that 
there  will  be  directions  in  the  plane  of  an  w-type  film  for 
which  the  carrier  mobility  exceeds  that  uf  balk  Si,  and 
other  directions  for  which  the  mobility  is  less  than  the 
bulk  value.  As  a  consequence,  it  is  possible  to  separate 
the  c  v.i.Httrtijt,  jnJtiifcly  dfcglauauun  in  me  si/ai2U3 
system,  due  to  differential  thermal  expansion  stresses 
from  degradation  due  to  other  causes.  The  results  in¬ 
dicate  that  differences  in  thermal  expansion  are  not  in 
themselves  sufficient  to  explain  the  reduced  mobility  in 
Si/A^Oj  films. 


Plots  of  Hall  mobility  (p)  vs  angle  (6)  measured  from 
the  [110]  Si  direction  are  shown  for  representative  sam¬ 
ples  in  Fig.  1.  For  each  sample,  a  least-squares  fit  of 
the  data  to  the  %.*Uun  (p)'1  <n  <mus2v  was  obtained 
(see  below).  An  anisotropy  factor  A  can  be  defined  ac¬ 
cording  to  A  --  100(p  -  P.J/P,  =  2006/fl,  where 

*****  "fa,.,  Ibt  maid MUIII  a„J  „,i,un,u,ii  , hoiJiliues 

determined  from  the  least-squares  fit,  respectively, 


j  —i - ito^cv-uvciy, 

and  pA  =  2 (Pmu  f  Pmln).  Table  I  shows  the  values  of  these 


The  data  reported  in  the  present  paper  are  from  mea¬ 
surements  taken  on  Si  grown  by  chemical  vapor  deposi¬ 
tion  on  a  A1203  orientation  ~3°  from  the  (1122)  plane, 
although  nearly  identical  effects  are  observed  on  the 
(1122)  plane  or  on  other  planes  near  the  (1122)  plane. 
Previous  orientation  studies2  have  shown  that  (221) 

Si  ii  ( 1122)  A1203  and  that  [IlO]  Si  H [lloo]  A1203  in  the 
plane  of  the  film.  The  reactor  system  and  methods  of 
growth  are  similar  to  those  described  previously, 3  and 
utilize  the  pyrolysis  of  silane  (SiH4). 


—  -  - - -  IIICOC 

parameters  for  a  number  of  samples.  The  observed 
anisotropy  is  seen  to  be  appreciable,  as  evidenced  by  an 
average  room -temperature  anisotropy  factor  of  approxi¬ 
mately  39%. 


For  this  experiment,  As -doped  films,  having  thick¬ 
nesses  from  1.  9  to  2.  5  pm  and  net  donor  concentrations 
of  (2  -7)x  1016  cm'3,  were  grown  at  1100  °C  at  a  growth 


In  order  to  explain  this  behavior,  we  have  calculated  the 
reduction  in  mobility  due  to  the  piezoresistanc»  iffi.rt 
resulting  from  the  difference  in  thermal  expansion  co¬ 
efficients  of  Si  and  A1203.  For  the  sake  of  completeness, 
we  have  also  taken  into  account  the  difference  in  thermal 
expansion  coefficients  in  A1203  parallel  and  perpendicu¬ 
lar  to  the  [0001J  axis.  The  details  of  the  calculations 
will  be  given  in  a  more  complete  account  of  this  work5- 
however,  using  the  data  of  Smith6  for  the  piezoresistance 
coefficients  (*„  =  -  102x  10'12,  54x  10'12,  „(4=-14 

l10  -  in  uni‘s  of  cmVdyn),  values  of  elastic  constants 

if:. .  —  1R  Rv  mil  ~  1 1  ~ 
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(Cn=16.5xlOn,  Cjj  =  6.  4x  10u,  C.,=  7.  93x10“  i 
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units  of  dyn/cm2)  from  Hearmon, 7  the  thermal  ex’pan- 


-T:U<U'  J  ^ra.nct.rs  ol  afcvfclai  (2a1,  Si/Ul2z)  a¥J3  films. 
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FIG.  1.  Variation  of  Hall  mobility  with  angle  from  the  [ll0|  Si 
direction  for  two  (221)  Sl/(1122)  Al203  films. 


sion  coefficients  asl  =  3.  9x  lO'V'C, 6  (lie)  =  9.  03 

x  10  0/  C,  «ai,o3  (1c)=8.31x  1O"0/ C, 9  and  a  temperature 
difference  of  1100  °C,  over  which  the  thermal  contraction 
takes  place,  we  find  that  (n/na)‘l=  l.  06375  -  0. 21891 
Xcos2 0,  where  fi0  is  the  mobility  in  the  absence  of 
stress  and  8  is  defined  as  above.  The  theory  predicts  an 
anisotropy  factor  of  41%,  in  good  agreement  with  the 
experimental  values  given  in  Table  I. 

An  interesting  feature  of  the  piezoresistance  in  this 
plane  is  that  the  mobility  is  predicted  to  be  larger  than 
the  unstressed  mobility  for  0°  <6  <37°,  and  smaller  than 
the  unstressed  mobility  for  37°  <  6i90°  (if  we  consider 
only  the  first  quadrant).  The  theory  predicts  at  37°  the 
effect  of  thermal  stress  on  the  mobility  is  zero,  so  that 
p(37°)  =  fi0,  the  unstressed  mobility.  The  difference 
between  the  bulk  mobility  and  the  value  of  mobility  mea  ¬ 
sured  at  this  angle  thus  represents  the  mobility  degrada¬ 
tion  due  to  the  causes  other  than  thermally  induced 
stress.  The  fact  that  the  value  of  p0=600  cm2/Vsec  is 
much  smaller  than  the  bulk  mobility10  strongly  suggests 
that  thermally  induced  stress  is  not  the  dominant 
mechanism  in  lowering  the  mobility. 

The  piezoresistance  theory  also  predicts  an  electric 
field  orthogonal  to  the  current  direction  Independent  of 
any  externally  applied  magnetic  field.  This  transverse 
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field  (Et)  is  given  theoretically  by  ET/jp0  =  0.21227 
xsin20,  where  j  is  the  current  density  and  p0  is  the  re¬ 
sistivity  in  the  absence  of  stress.  Accordingly,  we  have 
measured  this  transverse  field  (across  the  same  termi¬ 
nals  used  to  measure  the  Hall  voltage),  and  show  the 
quantity  ET/jp0  plotted  for  two  samples  in  Fig.  2.  The 
quantity  p0  is  found  from  p0=RH/p0,  where  RH  is  the 
Hall  constant  and  p0,  as  before,  is  the  mobility  at 
8=37  .  A  least-squares  fit  of  the  data  to  the  predicted 
equation  yields  an  average  value  for  the  coefficient  of 
sin28  of  0.  2,  in  good  agreement  with  theory. 

The  temperature  dependence  of  the  anisotropy  has  not 
been  examined  in  detail,  although  measurements  on 
selected  samples  at  77  °K  show  a  much  larger  anisotro¬ 
py,  as  evidenced  by  the  data  in  Table  I.  Such  an  in¬ 
crease  is  consistent  with  what  might  be  expected  for  the 
variation  of  piezoresistance  with  temperature.  Although 
detailed  measurements  of  all  of  the  piezoresistance  co¬ 
efficients  have  not  been  carried  out,  jtu  appears  to  be  an 
order  of  3  times  greater  at  77  °K  than  at  300  °K. 11  If  it 
is  assumed  that  all  of  the  coefficients  vary  by  the  same 
amount,  it  can  be  shown  that  the  anisotropy  factor  A 
would  be  expected  to  be  an  order  of  2. 7  times  greater 
at  77  °K.  The  data  on  two  samples  show  an  average  in¬ 
crease  of  a  factor  of  2.4,  consistent  with  that  predicted. 


0  (DEGREES) 

FIG.  2.  Variation  of  transverse  field  ET  with  angle  from  the 
lllo)  Si  direetion. 
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The  enhanced  mobility  predicted  for  some  current  di¬ 
rections  suggests  that  the  (221)  plane  of  Si/Al203  may 
yield  superior  electron  mobilities.  In  comparison  with 
(100)  Si  and  (111)  Si  grown  at  this  laboratory,  the  (221) 
plane  has  indeed  yielded  higher  mobilities.  It  would  be 
expected  that  some  types  of  devices  would  show 
improved  performance  if  oriented  in  such  a  way  as  to 
take  advantage  of  the  improved  mobility.  In  MOS  devices 
the  stress-induced  anisotropy  would  have  to  be  con¬ 
sidered  in  conjunction  with  the  anisotropies  previously 
found  in  n-  and  ft  -channel  devices,  12  which  are  a  conse¬ 
quence  of  surface  quantization  effects  in  the  channel 
resulting  from  the  high  electric  field  at  the  surface. 

The  mobility  anisotropy  discussed  in  this  letter  is  be¬ 
lieved  not  to  be  associated  with  surface  quantization 
since  there  is  no  applied  electric  field  normal  to  the 
surface  and  the  oxide  film  formed  during  the  aftergrowth 
annealing  sequence  has  been  removed  prior  to  the  mea¬ 
surement  of  film  properties.  In  addition,  the  high  donor 
concentrations  in  the  films  tend  to  mask  contributions  to 
the  electrical  conduction  from  surface  effects.  Experi¬ 
mentally,  we  have  seen  no  significant  effect  of  ambient 
or  surface  chemical  treatments  on  the  conductivity. 
However,  to  quantitatively  assess  the  magnitude  of  jxis- 
sible  surface  effects  would  require  the  measurement  of 
electrical  properties  in  an  MOS -device  configuration 
under  flat -band  conditions. 

The  good  agreement  found  between  the  experimentally 
determined  mobility  anisotropy  and  the  calculated 
piezoresistance  lends  confidence  to  the  assumption  that 
the  anisotropy  is  related  primarily  to  thermally  induced 
stress.  If  this  is  indeed  the  case,  further  experiments 
which  examine  the  variation  of  anisotropy  with  film 
thickness,  growth  conditions,  or  other  variables  may 
well  prove  valuable  in  determining  the  relation  of  stress 
to  these  variables. 
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ABSTRACT 

The  effect  of  stress  on  the  electrical  properties  for  a  number 
of  orientations  of  Si/Ai^O^  has  been  theoretically  determined. 

The  stress  model  developed  and  employed  is  based  upon  the  differen- 
L i n  1  in  thermal  expansion  between  Si  film  and  Ai^O^  substrate. 

The  anisotropy  in  substrate  thermal  expansion  coefficients  and  in 
Si  film  stresses  is  rigorously  included  and  treated  within  the 
framework  of  a  thin-film  t h i ck - sub s t rat e  approach.  The  phenomeno¬ 
logical  piezoresistance  formalism  is  used  to  determine  the  effect 
of  stress  on  resistivity  which  is  then  in  turn  related  to  the 
carrier  mobility.  In  calculating  the  mobility  change  due  to 
thermal  stress,  it  is  assumed  that  bulk  piezoresistance  coeffi¬ 
cients  can  be  employed.  Mobility  calculations  were  performed  for 
the  common  Si  orientations  -  (001),  (111)  and  (221)  -  and  for  a 
number  of  orientations  not  yet  grown  and  investigated  experimen¬ 
tally.  For  all  Si/A£203  orientations  arid  modes  of  epitaxy  treated, 
an  anisotropy  in  mobility  is  predicted  with  two  directions  of 
maximum  mobility  and  two  directions  of  minimum  mobility  in  the 
plane  of  the  Si  film.  For  n-type  Si  the  mobility  is  either 
reduced  or  enhanced  depending  upon  mode  of  epitaxy  and  orientation. 
For  p-type  Si  the  mobility  is  always  enhanced  and  for  some  orien¬ 
tations  by  a  significant  factor.  Experimental  data  are  reviewed 
for  n-type  (001)Si  and  ( 2 2 1 ) S i  and  found  to  agree  with  the  theo¬ 
retically  predicted  mobility  anisotropy,  lending  confidence  to  the 
use  of  the  thermal  stress-piezoresistance  model. 

*  Supported  in  part  by  ARPA  under  Order  1585  ,  monitored  by  USAMICOM, 
Redstone  Arsenal,  AL,  under  contract  DAAH01- 70-C- 131 1 
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INTRODUCTION 


The  principal  objective  of  the  present  investigation  has 
been  the  determination  of  the  effect  of  thermal  expansion 
stress  on  the  carrier  mobility  in  n-  and  p-type  Si/A£203  films. 

The  theoretical  analysis  of  mobility  and  mobility  anisotropies 
and  the  responsible  physical  mechanisms  began  as  a  response  to  the 
recent  experimental  discovery  at  this  laboratory  that  a  mobility 
anisotropy  was  present  in  (OOl)Si^,  (221)Si^  and  (lll)Si^3) 
heteroepitaxial  films  grown  on  A£203  and  would  probably  occur  in 
all  orientations  of  Si/AS^O^,  It  was  recognized  that  different 
physical  mechanisms  which  could  produce  some  form  of  mobility 
anisotropy  would,  in  general,  be  expected  to  lead  to  different 
magnitudes  of  anisotropy  and  to  different  orientations  for  the 
mobility  maxima  and  minima  in  the  plane  of  the  film.  Thus, 
detailed  studies  of  the  mobility  anisotropy  could,  in  principle, 
provide  a  more  powerful  means  than  has  previously  been  available 
for  determining  the  role  of  various  phenomena  in  establishing  the 
carrier  mobility  in  Si/A£203  films. 

Si/A£203  epitaxy  has  been  actively  pursued  for  a  number  of 
years  now  and  is  >f  considerable  practical  and  commercial  signifi¬ 
cance  for  certain  classes  of  devices.  However,  a  fundamental 
understanding  and  adequate  theoretical  models  capable  of  delineating 
quantitatively  the  various  physical  phenomena,  and  their  possible 
role  in  producing  a  mobility  anisotropy,  are  almost  completely 
lacking.  Since  stresses  in  Si/A£203  are  known  experimentally  to  be 
on  the  order  of  10  -10^  dyn/cm2  ancj  since  si  is  a  strongly 

p i e z o r e s i s t i ve  material,  it  was  felt  that  a  piezoresistance- effect 
anisotropy  model  would  be  a  useful  theoretical  tool  and  should 
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relate  substantially  to  experimental  mobility  anisotropies.  It 
is  our  view  that  the  room  temperature  stress  in  Si/AHn  fiIms  is 
predominantly  due  to  thermal  expansion  mismatoh  stress  rather  than 
to  lattice-constant  mismatch  or  other  growth  stresses.  This  view 
coupled  with  the  fact  that  thermal  stresses  appear  much  more 
amenable  to  theoretical  treatment  than  do  growth  stresses,  dictated 
that  an  investigation  in  terms  of  a  model  combining  thermal  expan¬ 
sion  stresses  and  piezoresistance  effects  would  be  n  useful  step  in 

better  understanding  and  improving  Si/At^  films.  This  is  the  basi 
approach  of  the  paper. 
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II.  MATHEMATICAL  DESCRIPTION  OF  HETEROEPITAXIAL  ORIENTATION  RELATIONS 
In  the  present  paper,  the  treatment  of  the  piezoresistance 
effect  for  Si/AH^  heteroepitaxial  films  will  require  that  various 
physical  property  tensors,  e.g.,  piezoresistance  tt,  elastic  constant 
C,  stress  T,  and  strain  e,  be  transformed  for  rotated  coordinate 
systems  in  the  Si  and  in  the  A*^  which  correctly  embody  the  relative 
orientations  of  the  film  and  substrate  for  the  different  growth  modes 
of  Si/A*203  heteroepitaxy.  A  rigorous  mathematical  formalism  for 
this  has  been  developed  which  facilitates  treating  an  arbitrary  Si 

fllm  orientation  and  which  is  applicable  to  all  the  major 

growth  modes  of  Si/A£203  heteroepitaxy. 

The  experimental  description  of  a  mode  of  epitaxy  is  accomplished 
by  employing  x-ray  diffraction  (using  back-reflection  Laue  and  full- 
circle  goniometer  techniques)  to  determine  the  relative  crystallo¬ 
graphic  relationships  between  film  and  substrate.  The  relative 
crystallographic  relationship  between  a  Si  film  and  the  A^O,  sub¬ 
strate  can  be  completely  and  uniquely  specified  by  giving  two  non- 
collinear  directions  in  the  substrate  which  are  parallel,  respectively, 
to  two  film  directions.  This  information  for  the  four  major 

growth  modes  for  Si/A£203  is  presented  in  Table  I  and  is  taken  from 
the  work  of  Manasevit,  Holder  and  MoudyC6).  The  relationships 
shown  for  C 2 2 1 ) S i / / C 1 1 2 2 ) A £ 2 0 3  in  Table  I  were  not  given  in  Ref.  6 
and  have  been  recently  obtained  by  Moudy(7^.  Note  that  Table  I 
presents  both  parallel  plane  and  parallel  direction  information. 

For  Si,  the  direction  [a,b,c]  is  normal  to  the  plane  (a,b,c)  having 
the  same  indices.  For  A£203,  the  direction  [h,k,i,£]  is  not  normal 
to  the  plane  (h,k,i,£)  unless  £  is  zero. 
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The  information  in  Table  I  allows  at  least  two  sets  of  two 
parallel  directions  to  be  obtained  for  each  of  the  modes  of  epitaxy. 
From  these  given  directions,  two  sets  of  three  parallel  directions 
can  be  constructed  which  are  mutually  orthogonal.  The  resulting 
orthogonal  direction  vectors  will  be  denoted  by  D*f,  D>f,  D>f  for 
the  Si  film  (f)  and  these  are  parallel,  respectively,  to  D^s,  d>s, 
D3S  in  the  AZ203  substrate  (s).  D>f,  D^f,  D>f  are  mutually  ortho¬ 
gonal  and  form  a  right-handed  set.  The  same  is  true  of  the  set 
DjS,  D2s,  D*s  in  the  substrate.  Thus  far,  the  D*f,  D*f,  D>f  and 
°1  ’  °2  »  °3 S  are  abstract  vectors.  In  order  to  proceed  further 
toward  a  mathematical  theory  of  epitaxial  orientation  relations, 
specific  coordinate  systems  in  both  Si  and  AA^  must  be  defined 
and  employed. 

Thus,  let  a  Cartesian  set  of  axes  (unit  vectors)  be  defined 

in  the  film  and  denoted  by  (Xf°,  Yf°,  Zf°) .  In  Si,  these  are 

referenced  to  the  cubic  crystal  axes.  Let  a  similar  Cartesian  set 

be  defined  in  the  substrate  and  be  denoted  by  (X  0  Y  0  7  °i 

J  s  * 

The  Cartesian  set  of  axes  (x/,  H  ° ,  z,.0)  is  referenced  to  the  At203 

hexagonal  axes  a*.  a*,  aj ,  aj ,  with  X  0  along  the  a*  axes;  Y  ° 

>  s  1  s 

90  deg  from  the  a1  axis  and  toward  the  a^  axis;  and  Zs°  along  the 
C  (or  a*)  axis.  In  general,  (X  °,  Y°  Z,°)  and  (X  °,  Y  °,  Z  °) 
differ  in  relative  orientation.  The  relative  orientation  between 
these  two  coordinate  systems  embodies  the  experimental  parallel 
relation  information  of  Table  I  and  determines  the  relative  crys¬ 
tallographic  orientation  of  Si  film  and  AA20  substrates.  The 
important  determining  constraint  is  that  D*f  in  the  Si  film  be 
parallel  to  the  D^s  in  the  substrate  for  i  =  1,2,3. 
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Let  (SF)  denote  the  matrix  transformation  between  (X  0  Y  0  z 

v  S  * 

and  (Xf  ,  Yf  ,  Zf  ).  Consider  the  coordinate  axes  to  be  spatially 
coincident;  (SF)  will,  therefore,  simply  be  a  rigid  body  rotation 
matrix.  The  matrix  (SF) ,  therefore,  is  orthogonal  and  (SF)  =  CSF) “  1 
We  then  have 


(SF) 


(1) 


vectors  D^f  and  D^s  for  i  =  1,2,3 


Now  let  the  parallel  direction 
be  expressed  in  Cartesian  components  relative  to  their  respective 
coordinate  systems.  We  then  form  the  matrices  (Imn)  and  (LMN) : 

'l\ 

and  (LMN)  = 


( lmn) 


m. 


m„ 


M, 


N1  ^ 


M, 


M. 


J 


(2) 


Where  11»  V  n !  =  Dlx*  ° 1 y ’  Dll>  h 


V  Nx 


s 

lx' 


s  n  s 
ly  lz 


-  D.  .  D.3,  D.s,  etc.  Thus, 


dmn)  is  comprised  of  row  vectors  relating  to  the  Si  film  and 

(LMN)  is  comprised  of  row  vectors  relating  to  the  A£203  substrate. 

In  forming  the  matrices  (lmn)  and  (LMN),  the  vectors  D>  have  been 
assumed  to  be  orthogonal  and  normalized  to  unit  magnitude. 

The  matrices  (lmn)  and  (LMN)  then  define  new  primed  coordinate 


systems  (Xf  ,  Yf,  Z  ')  and  (X 


i  » 

Yc  ,  Z  ) 

C  9  C  ' 
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The  primed  vectors  (Xf  ,  Yf  ,  Zf  )  must  be  parallel  and  thus  equal 
to  the  primed  vectors  (Xg '  ,  Yg '  ,  ij).  Using  Eqn.  (1),  we  have 


(4) 


which  yields  the  following  relation  between  (lmn)  and  (LMN)  matrices, 

(lmn)  =•  (LMN)  (SF)  and  (LMN)  =  (lmn)  (SF)  , 

which  must  be  obeyed  by  an  (lmn)  matrix  and  an  (LMN)  matrix  in  order 
to  be  consistent  with  the  parallel  relations  associated  with  a  given 
mode  of  epitaxy.  If  (lmn)  and  (LMN)  are  known,  then  (SF)  can  be 
determined  from 

(SF)  =  (LMN)  (lmn)  and  (sT)  -  (lmn)  (LMN). 


Next,  suppose  we  operate  on  (X  ' ,  Y  ' ,  Z  ' )  and  (X  '  Y  '  Z  ' 1 

s  s  s'  *■  f  »  f  »  L  f  } 

in  Eqn.  (3)  with  the  same  trans fomat ion  natrix  (a)  yielding 


)  ^e  the  transformation  which  now  connects  (X^°,Y^°  Z  °) 
and  (Xg  ,  Ys°>  Zs°),  we  have 

(SF)  -  (as)  (LMN)  (lmn)  (af)  *  (ct^)  (SF)  (af )  . 

In  general,  (SF)  will  be  different  from  (SF)  originally  adopted, 
but  on  general  principles  must  correspond  to  a  physically  equivalent 
relative  orientation  of  film  and  substrate  and,  thus,  to  the  same 
mode  of  epitaxy. 

By  considering  (af)  to  post-multiply  (lmn)  rather  than 

(Xf  ’  Yf°’  Zf0)  and  t0  similarly  post-multiply  (LMN),  we  have 

the  result  from  Eqn. (10)  tha  (lmn)  (af)  is  symmetry-equivalent  to 
(lmn)  and  (LMN)  (ag)  is  symmetry-equivalent  to  (LMN).  A  prescription 
has  thus  been  given  for  obtaining  all  equivalent  sets  of  (lmn),  (LMN) 
and  (SF)  matrices  from  an  original  arbitrarily  chosen  set. 

In  general,  the  existence  of  symmetry-equivalent  parallel 
directions  and  relative  orientations  of  film  and  substrate  presents 
no  problems.  However,  in  particular  cases  an  internal  consistency 
problem  of  practical  significance  can  arise.  Suppose,  for  example, 
that  parallel  relations  for  a  given  mode  of  epitaxy  have  been 
established.  This  amounts  to  having  arbitrarily  labeled  a  few 

principal  x-ray  reflections  which  then  'ixes  the  coordinate  systems 
and  directions. 

Next,  suppose  that  a  substrate  yielding  this  mode  of  epitaxy 
is  cut  a  few  degrees  "off-orientation."  X-ray  studies  providing 
the  required  characterization  of  the  substrate  and  film  surfaces 
can  then  lead  to  afferent  but  equivalent  pole  diagrams.  However,  only 
one  of  these  pole  diagrams  is  mathematically  consistent  with  the 
previously  determined  parallel  relations.  The  inconsistency  here 


is,  in  effect,  the  simultaneous  use  of  different  coordinate 
systems. 

The  reason  for  emphasizing  equivalent  f i lm / sub s t r a t e 
orientations  is  that  information  obtained  from  cry s t allogr apher s 
characterizing  a  mode  of  epitaxy  and  the  f i lm / sub s t ra t e  surfaces 
can  lead  to  confusion.  From  the  point  of  view  of  crystallography, 
a  direction  and  its  equivalent  are  equally  acceptable.  However, 
equivalence  is  not  adequate  for  mathematical  consistency. 

One  final  type  of  operation  will  now  be  discussed  which  will 
be  used  in  subsequent  sections  of  the  paper  to  treat  piezoresistance 
for  Si  films  grown  on  substrates  of  arbitrary  surface  orien¬ 

tation.  We  shall  identify  the  third  row  of  the  ( LMN )  and  (lmn) 
matrices  with  directions  normal  to  the  substrate  and  film  surfaces, 
respectively.  Thus,  if  a  substrate  is  cut  "off-orientation,”  a 
modified  set  of  matrices  is  required.  Suppose  the  direction  vector 
characterizing  the  original  substrate  surface  is  taken  from  the 
third  row  of  the  (LMN)  matrix  and  written  as  a  column  vector.  Then 
a  matrix  operation  (y)  acting  on  this  vector  yields  a  new  vector 
characterizing  the  o f f -o r i en t at  ion  substrate  surface. 


Thus , 


Since  the  operations  (y)  and  (B)  in  Eqns.(13)  and  (14) 
preserve  the  original  identical  mode  of  epitaxy,  it  must  be 


possible  to  represent  both  of  th 
as  in  Eqns .  (7)  and  (8)  ; 


em  in  terms  of  an  operation  (qi) 


(a)  (LMN)  =  (LMN)  (y)  ;  a  =  (LMN)  (y)  (LMN) 

(cx)  (imn)  =  (lmn)  (T) ;  «  =  (lmn)  (?)  (!?„), 

Using  Eqn.  (16),  we  find  that  Eqns.  (17)  are  consistent  and  an 

operation  (a)  can  be  defined  which  is  operationally  identical  to 

(Y)  and  (3)  as  indicated. 

The  parallel  relations  assumed  and  the  transformation  matrices 

(SF),  which  have  been  developed  for  the  four  Si/AJ^Oj  modes  of  major 

interest,  are  presented  in  Table  II.  One  mode  of  Si  epitaxy  on 

spinel  (MgA£204)  is  also  treated  in  Table  II. 

In  Table  II,  the  Zf  and  Zg  parallel  directions  are  given 

in  terms  of  normals  to  film  and  substrate  planes.  X  '  Y  '  and 
tt  f  f 

and  Xs  ,  Ys  are  given  explicitly  in  terms  of  direction  vectors. 

This  distinction  is  convenient  in  subsequent  sections  which  will 
focus  attention  on  film  and  substrate  surface  planes.  The  parallel 
directions  (planes)  for  the  AA^  substrate  are  given  both  in  a 
hexagonal  four-index  notation  and  in  Cartesian  coordinates. 

The  orientation  relations  (lll)Si// (1124) Af,^  and 
(221)Si// (1122)A£203  were  both  determined  experimentally  and  in 
Table  II  are  identified  as  belonging  to  the  identical  mode  of 
epitaxy  (Mode  II).  The  basis  for  this  is  the  following.  The 
(lll)Si//(H24)A£203  parallel  directions  and  transformation  matrix 
(SF)  for  Mode  II  were  assumed.  The  (lll)Si  surface  plane  was  then 
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(17‘ 


transformed  to  a  (221)Si  surface  plane.  The  associated  (LMN) 
matrix  for  (22 1 )  Si/ / (1 122)  Afl^O 3  was  then  calculated,  employing  the 
formulae  given  in  this  section,  and  compared  with  that  obtained 
directly  from  the  experimental  data.  The  two  (LMN)  matrices  were 

f 

found  to  be  substantially  identical.  For  example,  Zg  for 
(22 1 )  Si/ / (1 1 22)  is  0.(1122)  while  the  transformed  Zs  derived 

from  (lll)Si// (1124)A£203  is  j.(l,l,2,  2.03412).  The  angle  between 
these  two  planes  is  equal  to  0.31506  deg  and  is  smaller  than  the 
experimental  error  in  either  the  (lll)Si  or  ( 2 2 1 ) S i  orientation 
determinations.  We  shall,  based  on  a  low-index  argument,  assume 
that  the  (lll)Si  orientation  relation  is  the  more  fundamental  and 
take  this  Si/A^O^  case  as  the  "natural"  defining  growth  for  Mode  II 
epitaxy.  (221)Si  is  then  in  the  same  mode  of  epitaxy  but  is 
"rotated"  following  the  mathematical  formalism  given  in  this 


section. 


III.  STRESS  AND  STRAIN  IN  Si/AJ^Oj  EPITAXIAL  FILMS 

The  principal  objective  of  the  present  paper  is  the  theoret¬ 
ical  determination  of  stress-induced  anisotropies  in  carrier 
mobility  for  Si/Aj^Oj  films.  In  this  section  we  shall  discuss 
briefly  the  origin  of  stress  and  strain  in  Si/AS^Oj  films  and  then 
describe  the  stress  model  we  have  employed  to  calculate  thermal 
expansion  mismatch  stresses  for  use  in  the  piezoresistance-effect 
studies . 

The  origin  of  stress  and  strain  in  various  epitaxial  films 
has  been  discussed  by  Hoffman^®^  and  by  Chopra^^.  These  authors 
both  divide  the  stresses  in  epitaxial  films  into  two  categories: 

(1)  intrinsic  stresses  and  (2)  thermal  stresses. 

For  the  case  of  Si/AS.^  epitaxial  films,  our  view  is  that 
the  probable  major  origin  of  intrinsic  stress  is  lattice  misfit  or 
mismatch  between  film  and  substrate.  However,  it  is  not  presently 
possible  to  actually  calculate  the  lattice-mismatch  intrinsic  stress. 
In  general,  this  intrinsic  stress  is  probably  of  less  importance 
in  establishing  the  total  room  temperature  stress  level  in  the  films 
than  is  the  thermal  stress  to  be  discussed  subsequently.  Intuitively 
this  seems  a  reasonable  assumption,  in  that  the  films  are  grown  at 
elevated  temperatures  and  the  lattice  mismatch  or  growth  stress  can 
be  relieved  by  the  generation  of  dislocations.  High  dislocation 
densities  in  Si  films  have  been  observed,  consistent  with  this 
assumption . 

Let  us  next  consider  the  second  category  of  stresses  in 
epitaxial  films,  thermal  stresses.  Thermal  stress  in  Si/AJ^Oj  is 
due  to  the  combined  constraints  of  strong  film-substrate  bonding 
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and  dissimilar  coefficients  of  thermal  expansion  and  contraction 
for  film  and  substrate. 

One  formula  often  quoted*-9-*  is 

Tth  =  (aF'as^ATEF*  (18) 

Where  Tth  is  the  thermal  stress,  af  an  average  film  expansion 
coefficient,  and  ag  an  average  substrate  expansion  coefficient. 

AT  is  the  temperature  differential  betwee:.  growth  and  room  temperature 
and  Ep  is  an  appropriate  Young's  modulus  of  the  film. 

Equation  (18)  has  been  rather  widely  used  as  a  means  of 
estimating  film  stresses  due  to  thermal  expansion  mismatch. 

However,  there  are  several  qualifications  and  remarks  which  pertain 
to  this  equation. 

First,  the  equation  is  valid  only  for  the  thin-film  thick- 
substrate  case.  Jefkins*-10*  has  presented  formulae  for  the  stress 
m  Si  films  on  ASL 20 3  substrates  based  on  an  isotropic  simple  elastic 
model.  These  formulae  indicate  that  for  films  of  interest 
[  ~  2Um-  thi  ck  Si,  1 0  -  2  Omi  1  -  th  i  ck  A^Oj]  the  stress  relief  due  to 
bending  is  slight  and  the  isotropic  (average)  stress  calculated 
using  Eqn .(18)  agrees  to  within  about  1\%  .  This  small  difference 
would  justify  the  use  of  the  thin-film  thick-substrate  approximation 
for  calculation  of  the  average  stress  in  a  thin  film. 

The  second  point  to  be  noted  about  Eqn.  (18)  concerns  the 
question  of  anisotropy.  Anisotropy  is  generally  defined  as  the 
dependence  of  a  given  physical  property  upon  crystallographic 
orientation.  Thus  Eqn.  (18)  is  anisotropic  in  that  Young's  modulus 
Ef  depends  upon  the  Si  film  orientation  and  is  different  for,  say, 
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1  J  and  (1U,Si  fU"S-  ‘»  «•  context  of  thl.  paper 

anisotropy  .U,  be  taken  to  refer  pt,„llf  t0  the  depende„ce  j  . 

8lVBn  PhySlCal  Pr°Pert>'  “«“>«  orientation  in  the  piane  „f 

;  glVe"  fU"/SUbStrate  «■  this  sense  the  model  of 

,n’  (18)  lS  1SOtr°PiC  and  yields  information  on  anisotropic 
stresses  m  epitaxial  films. 

The  question  of  anisotropy  has  apparently  not  been  treated 

m  the  literature.  Schlotterer fl  pmmi  . 

employed  a  formula  similar  to 

Eqn-  (18)  in  treating  stress  in  Si  f i i 

Th  111  Sl  fllms  0n  M^2°4  ^bstrates. 

The  fum  orientation  was  properly  tak^n  •  . 

P  P  X  taken  into  account  in  calculating 

aPPr°Priate  Y°U"8'"  V  However,  the  appiied  thermal 

;  ;;  and  -«  — P-  as  is  appropriate 

for  (001)  and  (lll)Si  on  MgA£  0  T» 

t.  ,  *  2°4-  The  case  of  AJ1203  with  a  built-in 

rma  stress  anrsotropy  was  not  treated. 

*V,  -bstrates,  in  general,  yie.d  a  thermal  strain  which 

“  1"trlnslcall>'  anisotropic  because  of  an  anisotropy  i„  thermal 
expansion. 

mod  ,  w  .  The  thln-fil"  thick-substrate 

del  which  we  have  employed  , .  . 

explicitly  treats  substrate 

induced  anisotropy. 

our  model  is  based  upon  the  relative  thermal  expansion  between 

.  ‘  ‘  ^  AV*  SUbStPate  "d  ^  strain  in  particular  directions 
«  lc  this  expansion  induces,  following  the  thermal  expansion 

r  1  and  M*#3  is  described  by  second-rank  symmetric  tensors 

referenced  to  Cartesian  coordinate  systems. 

Ut  be  tbe  eXPanSi°"  «“»“*««  for  si  and  let  and  „ 

»e  the  thermal  expansion  coefficients  respectively  perpendicular  L 

Parallel  to  the  (0001)  axis  (or  c-axis)  for  M  „  ..  , 

ror  ax'2°3*  The  relative 

hermal  expansions  between  film  a„d  substrate  are  then 
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for  a  temperature  difference  AT.  Here,  Cartesian  axes  in  Si  are 
along  the  cubic  axes.  For  AA^  the  x-axis  is  parallel  to  aj,  the 
y-axis  is  at  90  deg  from  the  a^  axis  and  toward  the  a^  axis. 

The  z-axis  is  parallel  to  the  AJL^  c  (or  aj)-axis. 

Equation  (19)  is  written  in  terms  of  the  so-called  tensor 
thermal  expansions  a... .  For  convenience  we  now  switch  to  the  so- 
called  "engineering"  thermal  expansi ons C 12 }  Ei  given  by  e  = 

“ll*  E2  =  a2  2  *  £3  =  a33*  e4  =  2a23*  £5  n  2oii3»  and  e6  =  20^. 

In  this  contracted  single  index  notation,  Eqn.  (19)  yields 

=  (WAT»  e2=(as-ai^T  =  ei*e3=(0ts-°12:iAT»  e4  =  0*  e5  =  0»  and  e6  =  0 . 
For  application  to  the  Si/A£203  system  we  shall  require  the 

thermal  expansions  for  an  arbitrary  film/substrate  orientation. 

»  I  I 

Let  (X , Y , Z )  represent  a  new  set  of  Cartesian  axes  with  X' and  Y'  in 
the  plane  of  the  film  (substrate)  surface  and  z'  normal  to  the  film 
(substrate)  surface  plane.  The  coordinate  transformation  effecting 
this  change  of  axes  would  then  be  obtained  from  Eqn.  (3)  of  Section 
The  Si  thermal  expansion  matrix  is  invariant  to  an  arbitrary  (lmn) 
coordinate  transformation.  Therefore  the  relative  thermal  expansion 
matrix  on  the  right  hand  side  of  Eqn.  (19)  may  be  transformed 
directly  using  only  the  (LMN)  matrix  appropriate  to  the  substrate. 

Thus,  let  the  coordinate  transformation  matrix  for  the 
substrate  be 
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The  engineering  thermal  expansions  in  the  arbitrary  surface  plane 


(x'  and  y'  in  the  plane,  Z*  normal  to  the  plane)  are  then  given  by 


(1-N^Ej  ♦  N j2e 2 


(i-n22)Ei  .  n22E3 


(1-N3Vi  *  Nj% 


"  2N2N3  (el'e3-* 


"2N1N3  (ej-ej) 


-2NiN2(ere3) 


We  emphasize  that  the  above  are  thermal  expansions  and  have  not  yet 


been  identified  as  strains  per  se 


the  interface  and  constrained  in  the  X '  and  Y '  coordinate  directions 


in  the  plane  of  the  interface.  Thus,  the  thermal  expansions 
e  1 ,  e'2  and  given  in  Eqn.  (21)  are  to  be  regarded  as  applied 
strains.  There  are  no  other  components  of  applied  thermal  strain 


since  the  film/substrate  materials  are  free  to  expand  or  contract 


in  these  other  directions.  However,  the  e|,  e'2  and  applied 

thermal  strains  do  act  to  produce  other  components  of  strain 


through  Poisson-type  effects. 


The  stress  model 


in  terms  of  the  applied  thermal  strains  e1,  e 2  and  e^.  The  model 


The  film  and  substrate  are  assumed  to  be  firmly  bonded  across 


Will  be  in  terms  of  the  engineering  stresses  as  given  by  Nye^3), 

where  .1^,  Tyiy  =Tyy ,  Ty^-T^,  TyT^T^,  TyT^-T^ 

and  T6*T12=Txy.  The  matrix  ef  elastic  constants  also  is  assumed 
in  the  form  given  by  Nye. 

The  essence  of  the  model  is  simply  that  the  applied  thermal 
strains  cy  y  and  are  assumed  applied  to  the  Si  film.  These 
strains  interact  through  the  elastic  properties  of  the  Si  to  produce 
stresses  in  the  Si.  In  this  model  the  only  role  of  the  substrate 
is  to  produce  the  initial  applied  strains,  and  the  elastic  properties 
Of  the  substrate  need  not  be  considered. 

The  stresses  in  the  Si  film  T* ,  are  then  given  in  terms  of 
the  Si  elastic  constant  matrix  elements  c[.  and  the  strains  e!  by 


C11  C 1 2 


C13  C 1 4  C 1 5 


where  Ej,  e2>  are  the  applied  thermal  strains  and  e'4> 

are  Poisson  reaction  strains  to  be  determined  by  boundary  conditions. 

The  Si  film  surface  is  a  free  surface  and  must  have  no  stresses 

surface.  Thus  i  3,  7^,  Tg  are  all  zero  at  the  fra©  surface  of  the 
film.  Slnce  the  film  is  thin,  we  assume  the  stresses  are  approxi_ 

mately  constant  in  the  z  direction.  Continuity  then  allows  us  to 
set  these  three  stresses  to  zero  throughout  the  Si  film.  This  is 
the  second  basic  assumption  of  the  model. 


The  Si  elastic  constants  C  „  ( i  ■  1 ,  6;  j  =  1 , 6)  for  the  arbi¬ 


trary  Cartesian  coordinate  system  can  he  obtained  from 

those  for  the  original  crystallographic  system  along  the  Si  cubic 


axes  by  suitable  transformation.  Wortman  and  Evans(14)  have  worked 


out  the  tedious  transformation  based  upon  the  fourth  rank  elastic 

tensor  properties  and  we  have  employed  their  results  for  the  Si 


matrix  elements  C'.  .  . 

ij 


Assuming  the  to  be  known  for  the  desired  coordinate 


system,  and  the  applied  thermal  strains  e  e'  and  e'  to  be  known 
from  Eqn.  (21),  the  reaction  strains  e',  e^,  e',  not  given  by 
Eqn.  (21)  can  be  determined  by  matrix  inversion.  The  determination 
of  the  thermal  expansion  stresses  Tj,  T'  and  T'  required  for  the 
carrier  mobility-piezoresistance  study  can  then  be  carried  out. 

This  summarizes  the  stress  model  which  is  basic  to  the  present 
work.  Within  the  framework  of  what  is  essentially  a  thin-film 
thick-substrate  model,  the  anisotropy  in  thermal  expansion  coeffi¬ 
cients  and  the  anisotropy  in  the  Si  stresses  have  been  rigorously 
included  and  treated. 

The  stress  model  developed  is  a  very  simple  model. 

Our  present  opinion  is  that  for  the  Si/A£203  films  of  interest 
(1.5-2ym  thick  Si,  1 0  -  20mi  1  -  th  i  c  k  A^O-^  any  limitations  in 
application  of  the  stress  formula  will  be  due  in  large  measure  to 
features  which  lie  outside  the  framework  of  the  elastic  theory, 
rather  than  in  the  level  of  rigor  employed  in  working  this  particular 
problem. 
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IV.  THE  PIEZORESISTANCE  EFFECT  IN  Si 

The  basic  phenomenological  formalism  for  the  linear  piezo¬ 
resistance  effect  in  a  thin  "bulk-  Si  film  has  been  established 
for  some  time  now.  However,  there  has  previously  been  no  develop¬ 
ment  and  application  to  heteroepitaxial  substrate-induced  mobility 
anisotropies.  In  this  section  the  formal  equations  which  we  have 
employed  in  calculations  of  mobility  anisotropy  are  discussed. 
Most  of  the  emphasis  is  upon  a  linear  model  but  non-linear 


corrections  of  interest  for  the  high  stress  levels  in  Si/AH203  films 


are  also  presented  and  employed. 

The  linear  piezoresistance  effect  has  been  discussed  by  various 
authors  *  .  These  authors  begin  by  assuming  that  the  electric 

field  components,  E.",  are  functions  of  the  current  density  components, 
Jj  and  the  stress  components,  T^ .  E'  is  then  expanded  in  a  series 
about  the  state  of  zero  stress  and  zero  current.  The  differential 
increments  in  the  series  expansion  are  then  replaced  by  the  variables 
themselves  yielding  the  basic  tensor  equation 


b  •  ~  p  ' .  J  '  +  it  ,  j  "  t  ' 

1  ij  1  ljkfc  j  kJl 


(23) 


For  cubic  crystals  the  resistivity  tensor  p.".  is  diagonal  for  zero 


stress  and  has  only  the  components  where  Pq  is  the  zero- 


stress  resistivity  and  6..  is  the  Kronecker  delta.  For  these 


materials  it  is  customary  to  factor  out  the  zero-stress  resistivity 


from  the  tt/'  .  Then 


Ei*  =  P 


(6  .  .  +  -iiiLLx  '  )  j  ' 

(  1J  P0  k£J  Jj  * 
where  j  k  £.  /  P  0  are  the  Piezo  resistance  coefficients  as  defined 


(24) 
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and  employed  by  most  authors  in  quoting  experimental  values  for 
the  coefficients.  In  applications,  Eqn.  (24)  is  usually  shortened 
by  using  the  single-subscript  notation  for  the  stress  components 
and  a  doub le- sub s cr ipt  notation  for  the  piezoresistance  coefficients 
[11^1*  22-*-2,  33-*3,  23-M,  13+5,  1 2+6  ]  ^  ^  ^  ^  .  The  double-subscript 
17 ' s  then  can  be  displayed  as  ?.  6x6  matrix,  which  for  Si  and  for  a 
Cartesian  coordinate  systen  directed  along  the  cubic  crystal  axes 
is  of  the  form 


(25 


The  three  constants  tt  1  ^ ,  tt  12,  tt44  then  completely  characterize  the 

bulk  piezoresistance  matrix  for  Si. 

Equations  (23)  and  (24)  refer  to  an  arbitrarily  oriented 

Cartesian  coordinate  system, and  the  associated  matrix  of  coefficients 

will  not  be  of  the  simple  form  in  Eqn.  (25).  Pfann  and  Thurston ( 16 5 
C 17  1 

and  Thurston  have  presented  tables  listing  all  -n'  coefficients 

for  rotated  (transformed)  axes.  In  terms  of  the  primed  (rotated) 
coordinate  system,  no  elements  of  the  tt'  matrix  are  required  to 
be  zero  and  the  matrix  is,  in  general,  asymmetric  in  that  some 
17  '  s  differ  from  the  tt  i  '  s  by  a  factor  of  two. 

For  an  arbitrarily  oriented  coordinate  system  (X'Y'Zl  or 
(XijXzjXj)  in  which  X'  (X^)  and  Y'  (XC,)  are  in  the  plane  of  the 
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Si/A£203  film  and  Z'(X')  is  normal  to  the  Si/A£203  film,  the 
Eqn.  (24)  then  becomes 


E^0 

= 

E^o 

= 

V\Vl'**6V2%\6T6>  ■ 

E?Po 

= 

El^o 

= 

E?*o 

E3/P0 

= 

The  it..  include  the  resistivity  Pq  and  in  some  cases  a  factor  of 
two  in  the  definition.  In  Eqns .  (26)  we  have  assumed  that  the 
stresses  Tj,  T2  and  T&  are  the  only  non-zero  stress  components 
in  the  Si  as  discussed  in  Section  III.  In  the  first  three  of 
Eqns.  (26)  it  was  assumed  that  the  current  flow  was  solely  in 
the  1  direction  (only  j'  non- zero),  and  in  the  last  three  of 
Eqns.  (26)  the  current  is  solely  in  the  3  direction  (only  j'  non¬ 
zero).  These  two  cases  are  adequate  to  cover  all  experimental 
device  configurations  of  interest. 

The  first  and  last  of  Eqns.  (26)  can  be  written  in  a  form 

directly  related  to  the  carrier  mobility.  These  two  equations 
deal  with  so-called  longitudinal  piezoresistance  effects  in  that 
the  field  and  current  components  are  both  in  the  same  direction. 

We  let  Pj  and  p3  be  the  resistivities  in  the  1  and  3  coordinate 
directions,  respectively.  Then,  in  terms  of  the  zero-stress 

resistivity  p  , 
o 
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P1  _  po+Apl 


Ap: 


j:  p 

lpo 


’  p, 


7TnTi+  ffi2T2%  "iVJ  « 


P3  po+Ap3 


APJ 


j:p 
3  Ko 


^iV  ^2T5  +  Tr36T6  • 


(27; 


The  Hall  mobility  y  in  a  given  direction  can  be  related  to  the 
resistivity  p  by  y  =  RH/p,  where  R^=  the  Hall  constant.  Thus 


Vuo 


l+ApJ/Po 


;  vJ/vo 


l+Ap'/p 

3  o 


(28) 


where  APj/pQ  and  Ap^/po  are  as  defined  in  Eqn.  (27).  Our  primary 
interest  in  the  Hall  mobility  will  be  in  mobilities  for  current 
flows  in  the  plane  of  the  Si  film.  In  subsequent  discussion, 
a  mobility  y  without  a  subscript  will  be  taken  to  relate  to  y^. 
Eqns.  (28)  are  fundamental  and,  in  the  linear  piezo¬ 

resistance  model,  all  connections  between  mobilities  and  resis¬ 
tivities  will  assume  y  =  R^p  and  Eqns.  (28). 

The  above  discussion  relates  to  a  linear  piezoresistance 
model  in  that  the  change  in  resistivity  is  proportional  to  the 
stress.  However*  since  stress  levels  in  Si/Ai^Oj  films  are  high 
and  relate  to  strains  much  larger  than  those  usually  encountered 
in  "bulk"  Si  piezoresistance  applications,  it  is  of  interest  to 

consider  non-linear  corrections  to  the  linear  model. 

A  non-linear  theory  of  piezoresistance  does  not  presently 
exist.  The  rigorous  development  of  such  a  theory  appears  to  be 
a  formidable  task  and  will  not  be  attempted  here.  We  shall 
simply  consider  corrections  to  the  longitudinal  piezoresistance 
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Pi  =  ^i/Jj  and  =  E^/J'  which,  while  not  rigorously  correct, 
are  probably  an  improvement  on  the  linear  model  for  those  cases 
in  which  the  predicted  change  in  resistivity  due  to  stress  is 
substantial . 


Schlotterer 
stress  using 


(11) 


has  related  the  change  in  resistivity  to  the 


—  =  k  T+k  T2  , 
P  12 


(18) 


following  Mason,  Forst  and  Torni  1  lo  1  J .  The  latter  authors 
employed  a  Si  many-valley  model  and  considered  the  effects  of  a 


uniaxial  stress.  No  direct  use  was  made  of  the  .  coefficients. 


Schlotterer  then  assumed,  based  upon  the  work  of  Sanchez  and 
Wright^19-*,  that  for  practical  use  the  relation  k  *0.5  k  2 


was 


2  “•*'  "1 

valid.  For  present  purposes,  we  shall  simply  assume  this  relation 
is  approximately  correct  and  will  not  attempt  to  justify  it. 

The  procedure  followed  by  Schlotterer  would,  in  our  case  of 
anisotropic  stress,  yield  mobilities 


Vo 


and 


P,/U 

o  0 


where  the  implied  sum  on  j  includes  j=l,2  and  6  for  Si/A£203  films. 


The  primes  on  tt  and  T  signify  that  it  and  T  are  to  be  calculated 
for  the  arbitrary  coordinate  system  under  consideration. 
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(29) 


(30) 


It  might  be  argued  that  if  Eon  f2<n  ,c 

4  •  is  approximately 

for  the  stress  and  strain  levels  of  Ref.  (19),  then  the 
larger  stresses  in  Si/A^o,  fll„  may  requlre  t„.  Incluslon  Qf 

hlgher  »e  have  derived  a  correction 

“hlCh  "  ln  th°  SPirU  °f  ^  but  which  includes  higher 

order  terms. 

We  consider  E,„.  (23)  for  j-i.  Differentiating  with  respect 

t.T  „  l,  A.  _  ? 


t0  Tk£  we  obtain 


7t  : 


6(El'/Jir)  -  Pi  (  )  («Tk't, 


1  1 


Schlotterpr  ^  ^  r-* 

presents  a  graph  of  a.,  coefficients  for  a  range 

O  f  r’  n  r*  ■v'  i  - - 


ij  - iur  a  range 

Of  carrier  concentrations.  The  experimental  Ps  are  approximately 
constant  for  moderate  changes  in  carrier  concentration.  This 
suggests  to  us  that  the  terms  C^ut/p'l)  in  Epn.  (31,  are 

approximately  constant  for  moderate  changes  in  resistivity  pf.. 

We  shall  assume  this  to  be  the  case  Nnf  “ 

e  case.  Noting  that  (E' /J ')~p' m  =  p' 

Eqn.  (31)  integrates  to  yield  **  * 


e  ♦C'ijTj) 


111  tHe  C°ntraCted  subs cr i pt  notation.  Here  p'(Tr  =  0)=p  is  the 

zero-stress  (scalar)  resistivity.  The  longitudinal  mobilities 
for  this  model  are  then 


Vl/Vo 


1  V, 


1 


e*f”3'jTp 


E,"S'  ‘28)’  (3°>  anii  (33)  thus  represent  mobilities  based  on 

(1st  order),  second,  and  exponential  order  stress  models. 
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(31) 


(32) 


(33) 


respectively,  and  are  all  mutually  consistent  to  within  their 
respective  orders. 


Numerical  data  far  comparison  of  mobilities  in  the  linear 
and  non-linear  models  are  given  in  Table  III.  The  linear 
mobility  is  taken  as  the  independent  variable  and  the  non-linear 
mobilities  are  then  calculated  from  Eqns.  (28),  (30)  and  (33). 


In  this  paper,  the  linear  model  will  be  used  in  all  analytic 
expressions  and  formulae  presented  for  the  various  orientations 
of  p-  and  n-type  Si/A^o,  films  studied.  In  most  cases  where  the 
resulting  linear  and  non-linear  numerical  mobilities  differ  sig¬ 
nificantly,  both  the  linear  and  exponential  order  mobilities  will 
be  given. 

A  more  mathematically  rigorous  non-li„ear  model  than  in 
fqns.  (31)  and  (32)  has  also  been  investigated.  This  model  is,  in 
principle,  applicable  to  the  transverse  as  well  as  the  longitudina 
piezo  resistances.  The  result  is  a  complicated  matrix  analog  of 
the  scalar  Eqn.  (32)  and  numerical  data  cannot  be  obtained  in 
closed  form.  It  thus  appears  to  be  without  profit  to  further 
pursue  this  approach  at  this  time. 

In  concluding  the  discussion  of  non-linear  p i e tores i s tan ce 
effects,  the  position  taken  will  be  re-emphasized.  Neither  of 
the  correction  terms  employed  in  Table  III  is  considered  to  be  rigo 

ously  correct.  Nonetheless,  in  those  cases  where  linear  and  non¬ 
linear  mobilities  differ  markedly,  it  is  felt  that  the  non-linear 
mobility  (in  particular  the  exponential  order  mobility)  is  much 
more  accurate  than  the  linear  mobility  and  must  be  employed. 
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A  prime  was  used  in  all  piezoresistance  formulae  presented 
thus  far.  The  prime  signifies  that  the  it's  and  T's  are  referenced 
to  an  arbitrary  fixed  coordinate  system.  We  now  wish  to  become 
specific  to  the  case  of  Si/AS^Oj  films.  For  a  given  Si/A£203 
film  growth,  we  are  interested  in  the  anisotropy  of  piezoresistance 
and  mobility  as  a  function  of  orientation  or  direction  in  the  plane 
of  the  film. 

Thus,  let  x'y'z'  represent  reference  axes  with  X"  and  Y'  in 
the  plane  of  the  film  and  Z'  normal  to  the  plane.  The  coordinate 
transformation  relating  this  set  of  reference  axes  to  the  Si 
crystal  axes  will  be  denoted  by  (£mn)  as  in  Eqns .  (3)  or  (8). 

Now  let  (a)  be  a  matrix  representing  a  rotation  about  the  Z 
direction.  Then 


(a) 


( 


co  s  0 


s  inO 


s  in0 
c  o  s  0 


^  0  0 


where  the  sense  of  the  rotation  is  that  a  positive  angle  0  rotates 
the  axis  toward  the  Y'  axis.  Following  Eqn.  (8),  the  (£mn)" 
type  matrix  explicitly  involving  the  angle  0  and  referenced  to 
the  x'y'  z'  axes  i s 


(£mn)  =  (a  )  ( £mn)  . 
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IS  rotation  could  also  be  applied  to  the  Ai.,03  substrate 
coordinate  system, but  for  present  purposes  this  is  unnecessary. 

Let  us  return  to  Eqns .  (26).  Now  the  prime  in  the  equations 
will  be  interpreted  to  refer  to  a  coordinate  system  at  an  angle  0 
from  the  reference  coordinate  system  with  X  and  V  in  the  Si  plane 
and  Z  normal  to  the  Si  film  plane.  The  n  coefficients  are  then 
obtained  using  (tun)' in  Eqn.  (35)  and  the  table  in  Ref.  (15). 

The  n  coefficients  will  then  be  explicit  functions  of  0  the 

angle  of  rotation.  The  stresses  TJ,  T'  and  TJ  are  similarly 
interpreted  and  are  given  by 


T-  _,VT2,  /i-T2 

1  ^  2  '  +  (  2  ^  cos20  +  T^sin20, 


•  VT2 


TrT2. 


r'  r  1  ,12 

2  2  J  ‘  ( — 2 - 3  cos20  -  T6sin20, 

T1"T2 

T6  =  — 2 )sin26  +  T  cos20, 


vhere  now  and  are  the  stresses  in  the  film  at  an  angle  0 

from  the  reference  and  T,,^,^  are  the  reference  stresses  at  0  =  0 


APPLICATION  TO  (001),  (221),  and  (lll)Si  FILM  GROWTH 
m  this  section  the  calculation  of  mobility  anisotropy  for 
Si  film  growths  of  three  different  orientations,  (001),  (221), 
and  (in),  will  be  discussed.  Measured  experimental  anisotropies 
and  calculated  theoretical  results  for  (001)31^  and  for  (221, Si'2) 
have  been  previously  reported  by  the  author.  The  purpose  of  the 
present  section  if  to  present  a  more  detailed  theoretical  description 
and  provide  information  not  reported  in  the  earlier  accounts. 


A.  (001)si/(01l2)At203  (Mode  I  Epitaxy) 

The  case  of  (00 1)  Si// (01 12)  At^  is  an  example  of  Mode  I 

epitaxy  as  defined  in  Section  II  and  in  Table  II.  For  calculation 

of  the  pie  tores i s  tance  and  mobility  anisotropy  for  this  orientation 

we  shall  adopt  reference  coordinate  systems  given  as  (x'y'  z') 
re,'  f  f'  f 


and  (Xs,Ys,Z^),  respectively,  in  Table  II.  The  (LMN)  matrix 
relating  this  coordinate  system  in  the  A»2(>3  substrate  to  the 
original  Cartesian  crystal  axes  is  then  given  by  (Eqn.  (3)) 


c 

0 

0  > 

/* 

0  o  ' 

(LMN)  = 

l0 

a 

1 

/  n 

l 

cn 

0/3  i 

=  0 

+.535703144  -.8444063841 

\° 

l 

n/3 

V 

+  .8444063841  +.535703144  j 

where  r)  is  a 

normali zation 

constant  = 

(1/3+  (a/c)  2)"z  and  a  (c/a) 

ratio  of  2.730159 

has 

been 

assumed  for 

^2^3*  ^he  (^mn)  matrix 

relating  to  our  reference  axes  in  Si  and  the  (tan)'  matrix  incor¬ 
porating  the  rotation  about  the  z'  axis  are  both  trivial  and  are 
given  below. 
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(3  7) 


,■ 


(£mn) 


(Amn) 


Employing  the  (Ann)'  matrix  in  Eqn.  (38),  the  i r  transformation 
table  in  Ref.  (15),  and  Eqns.  (.>6)  we  find  the  piezoresistance 
equations  (Eqns.  (26))  become 


(38) 


E  APi 

(—  'D  = 

o  1  Ko 


-  ( 


VT2. 


e2  TrT2 

(p  J-)  =  ■(  2 

Mo  1  z 


T  -T 

(irll  +  1T12)  +  2 }  (^11-^12)00520 


+  T67T44  sin2G 


)  Orn-Tr12)  s  :L  n  2  6  +  Tir  Cos20, 


(39) 


(40) 


-1)  = 


ap; 


PoJ3  P, 


-  ,W.. 

(  2  )2tt12. 


(41) 


For  (001)81,  (Z3/P0J\)  -  tfL/P0J'5)  -  (Eyp^)  -  0.  Eqns.  (39) 

and  (41)  reduce  to  Schlotterer ' s ( 11 5  Eqns.  (5a)  and  (5c)  in  the 

limit  where  Tj=T  and  T  =0.  In  the  above  equations  ir,  ir  it 

If  12’  44 

are  the  basic  piezoresistance  constants  shown  in  Eqn.  (25)  and 
referenced  to  the  Si  crystal  axes. 

We  next  solve  for  the  stresses  Tj,T2  and  T&  from  Eqns.  (22). 

Here  the  primes  on  the  C..'s  and  the  e^s  relate  to  0  =  0  and  hence 
to  the  (£mn)  matrix  in  Eqn.  (38)  and  to  the  (LMN)  matrix  in  Eqn.  (37), 
respectively.  In  accordance  with  our  notation,  the  primes  on  the 
stresses  are  to  be  suppressed.  We  obtain 
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rj-lvtirsUt 


Kfir&n&iS&Jti-  ^JL'MXrSAIZiS*. 


T1 

C 1 1 E 1  + 

C12e2 

+  C12£3 

T2 

=  C 1  2  £  1  + 

C 1 1 e  2 

+  C12£3 

T3 

=  C12eJ  + 

C12e2 

+  Clle3 

T4 

c44e4» 

T5 

=  C  4  4  e  5  ’ 

T6 

C44e6* 

(42) 


where  e^,  e2,  and  are  applied  thermal  strains  and  e'  e'  and  e' 

OH  5 

are  Poisson  reaction  strains. 

In  accordance  with  our  stress  model,  T_=T  =T  =0.  We  then 

3  4  5 

find  that  for  this  orientation  of  film  and  substrate,  the  reaction 
strains  ej  and  ej  are  zero.  Solving  for  Tj,  1 2,  T6,  we  obtain 
the  implicit  relations 


Tl+T2  _  C11  +C11C12‘2C122  29  90 

2  TCTj - -  It2-Nj  -N2  )e,  .  (N^tNjhej], 


T1-T2  (circi2>  .  2  2 

2 - 2 -  'V)  (Ej-Ej), 


(43) 


6  -  -2C44  NjN2  (Ej-e,), 


where  Nj  and  N2  are  defined  in  Eqns.  (3)  and  (20)  and  given 
numerically  in  Eqn.  (37).  The  crystallographic  axis  strains 
ei  and  e3  have  been  given  in  Section  III  and  are 

£1  =  (VVAT> 

(44) 

E3  =  (as-a2^AT> 
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where  ag  is  the  thermal  expansion  coefficient  for  Si  and  c^,  a2 
are  the  thermal  expansion  coefficients  for  A^2°3  in  directions 
perpendicular  and  parallel,  respectively,  to  the  c-axis. 

The  stresses  in  Eqn.  (43)  can  now  be  evaluated.  For  elastic 
constant  data  for  Si  we  use  that  given  by  Hearmon^^,  por  ^e  §i 
thermal  expansion,  an  average  value  of  ag=3 . 9xl0_6/deg  C^11)  was 
employed.  For  our  primary  data  on  A^Oj  thermal  expansion  we  use 
that  of  Austin  with  (^  =  8  .  31xl0-6/deg  C  and  a2»9 . 03xl0_6/dog  C 

Thus  ( a  i  ~  a2 )  =  _  0  •  72x  10  /deg  C  for  Ai^O^.  Recent  measurements 
at  this  laboratory  using  a  differential  technique  have  given  a 
difference  in  thermal  expansion  coefficients  a^-o^  of  (-  1 . 0  8  +  0  .  12 )  xlO -fVdqg  c 
for  the  Ai^Og  substrate.  Calculations  have  been  performed  using  both 
sets  of  expansion  data  for  AJ^O^  as  a  means  of  determining  the 
sensitivity  of  the  calculated  results  to  the  expansion  data^^. 

The  results  are  given  below  for  AT»1100C  deg. 


Case  1;  (oij-o^)  =  -0 . 72xl0-6/deg  C^21) 
2-(T^  +  T2)  =  -0.9206x10^  dyn/cm^ 

I(T1"T2)  =  +0-2852x109  dyn/cm2 
T6  "  °-° 

Case  2;  (Oj-c^)  =  -  1 . 08x  10 -6/de  g  c(22) 

y(T^+T2)  =  -0.9105x10*0  dyn/cm^ 

2"(Tj-T0)  =  +0.4278xl09  dyn/cm^ 

T6  -  O’0 
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(45) 


‘■•mini .. 


Substitution  of  these  values  for  the 
piezoresistance  coefficients  for  Si  from 


stresses  and  values  of 
Smith (24. 25)  yields  the 


results  shown  in  Table  IV. 

-et  us  first  consider  the  p-type  Si  results  of  Table  IV 

For  this  (00 1 ) Si  orientation  the  p-type  piezoresistance  effects 

(Api/po,  ^2/ P o J 1  *  and  Ap3/p0)  are  all  quite  small  and  would  yield 
variations  in  mobility  much  smaller  than  the  experimental  scatter 

in  data  normally  present  in  such  measurements.  For  all  practical 


purposes  then,  piezoresistance  in  p-type  (001)Si  films  can  apparently 
be  ignored.  On  the  other  hand,  we  shall  find  that  p-type  piezo¬ 
resistance  effects  for  some  other  Si  orientations  are  substantially 
larger  than  the  n-type  piezoresistance  effects. 

Next  consider  the  n-type  Si  longitudinal  piezoresistance  data 
Api/P0  of  Table  IV.  Since  the  mobility  y  can  be  related  to  the 

resistivity  p1  by  M-Rn/Pj,  where  f*H  =  the  Hall  constant,  the  theo¬ 
retical  mobility  can  be  written 


Case  1 

M/M0  =  [l  +  Ap|/po]-1  =  [ 1 . 441924- . 044488  cos20]_1, 


M/yQ  =  [ l+Ap j/Pq] " 1  =  [1.437032-. 066732  cos20["1, 

where  yQ  is  the  zero-stress  mobility.  From  Eqn.  (46)  we  note  that 
the  mobility  will  be  a  maximum  along  the  [100]Si  direction  and 
will  be  a  minimum  90  deg  away  along  the  [010]Si  direction. 

The  amount  of  anisotropy  --  the  parameter  which  can  be  most 
conveniently  and  definitively  extracted  from  experimental  data  -- 
can  be  described  conveniently  by  a  parameter  A  which  we  define  as 
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A  =  ^max'^min^V  Wh6re  PA  =  l/ 2  ^max^min5  is  an  aPProximate 
average  mobility  in  the  plane ^26^.  In  terms  of  these  parameters 

we  find  A  =  6.17%  and  yA  =  0.694yQ  for  Case  1  and  A  =  9.29%  and 
PA  =0.697yo  for  Case  2  values  of  thermal  expansion  data.  V.’e  thus 
find  an  anisotropy  in  mobility  of  from  6.2%  to  9.3%  superimposed 
on  a  substantial  ~30%  overall  theoretical  reduction  in  mobility 
for  n-type  (OOl)Si. 

The  theoretically  predicted  transverse  piezoresistance  effect, 
corresponding  to  an  electric  field  in  the  plane  of  the  film  and 
orthogonal  to  the  current  J'  direction,  is  also  given  in  Table  IV. 
This  effect  is  characterized  by  the  magnitude  of  the  coefficient 
of  the  sin20  term  which  is  +0.044488  and  +0.066732,  respectively, 
for  Case  1  and  Case  2. 

The  longitudinal  piezoresistance  effect  Ap^/po  and  the 
transverse  effect  E^^o^l  have  both  been  measured  experimentally 
and  reported  in  Ref.  (1).  The  values  of  the  anisotropy  parameter  A 
determined  experimentally  ranged  from  7.6%  to  11.7%,  with  an 
average  value  of  9.5%.  The  measured  transverse-effect  data,  due  to 
the  smallness  of  the  voltage  being  measured,  exhibited  more  scatter 
than  did  the  longitudinal  effect.  A  least-squares  fit  to  the  data 
yielded  a  sin26  coefficient  of  0.08.  The  theoretically  predicted 
longitudinal  and  transverse  piezoresistance  effects  for  Case  2 
thermal  expansion  data  are  thus  in  reasonable  agreement  with 
experiment.  Thus,  the  conclusion  in  Ref.  (1)  was  that  the  mobility 
anisotropy  and  transverse  effect  could  be  substantially  accounted 
for  in  terms  of  substrate-induced  stresses  without  recourse  to 
other  phenomena.  For  a  detailed  discussion  the  reader  is  referred 


B.  (221)/Si/ (1122)A£203  (Mode  II  Epitaxy) 

The  case  of  (22 1)  Si/ (1 122)  AS^O^  is  an  "off-orientation" 
example  of  Mode  II  epitaxy  as  discussed  in  Section  II.  In  calcu¬ 
lating  the  piezoresistance  effects  for  this  film  orientation,  we 
shall  adopt  reference  coordinate  systems  given  as  (X£,Y£,Z~)  and 
(XI>Y',Z') ,  respectively,  in  Table  II  under  the  (221) Si/ (11 2 2) Afl,_0 , 
heading . 

The  parallel  relations  in  Table  II  lead  to  the  following  (LMN) 
matrix  relative  to  our  reference  axes: 


(LMN) 


// 3 

2 

-1/2 

0 

(a/c) 

/T  (a/c) 

-l 

2n 

2n 

0 

/3 

(a/c 

V  2n 

2n 

n 

N  =  0.0 


N2.-0.93899, 


where  n  is  a  normalization  constant  =  (l+(a/c)  )2  and  (c/a) =2 . 730159 
was  assumed. 

The  (imn)'  matrix  incorporating  the  rotation  about  the  Z£  axis 
is  given  by 


(S,mn) 


+3cos9+sin0 

-3sin0+cos0 

2/2 


-3cos0+sin0  -4sin0  \ 

+  3sin0  +  cos0  - 4 c o s 0  /  , 

2/2  /I  / 


where  0  is  measured  from  the  [110] Si  direction  toward  the  [ 1 1 4 ] S i 
direction.  Using  this  (S,mn)'  matrix  and  the  transformation  properties 
of  the  it's,  the  ( 2  2  1 )  Si/ / (1 1 2  2 )  AS,  ^0^  piezoresistance  equations 
become 


ii] 


We  next  solve  for  the  stresses  T  T  and  T  f  c 

1  *  ‘2  and  *6  from  Eqns.  (22) 

and  for  this  f ilm/substrate  orientation  obtain 


Cllel  *  C12e2  *  C13e3  4  S-tV 


C'l2<  4  C22e2  4  C'3e'  ♦  C'4e', 


T3  =  C13el  4  C2V2'  *  C'3e'  .  C34e', 

T4  =  Cl>l'  4  C2>2'  *  C34 e3  *  C4'4E4', 

T5  =  C55e5  4  C56 e6 * 

T6  ■  C56  e5  4  C66  e6 ' 

"here  the  pri.es  on  the  C...s  refer  to  ..rotated-  elastic  coeffi- 

ClentS  USi"8  °f  E^-  («)  »ith  the  angle  6-0  and  the  trans¬ 

lation  table  of  Ref.  (14).  This  system  ,f  equations  J§  fuf£i_ 

cient ly  complicated  to  render  an  analytic  solution  of  Eqns  (501 

ln  termS  °f  the  aPPHed  V  and  e'  impractical.  „e 

have  therefore  simply  solved  the  equations  numerically,  using  the 
elastic  constants  for  Si. 

“sing  Eqns.  (21)  and  the  (LMK)  matrix  in  Eqn.  (47,  we  obtain 

H  "  f“s-“l)AT.  E!'(“,-»|)«  4  0.8817095  (Uj-a^iT,  and  E-.0i 
Again  setting  lyr^T^o,  and  solving  for  Tj  and  T?  we  find  the 
following  implicit  relations: 

T14T2  , 

2  -  t22.338504  (as-a1JiT  *  9 . 226 1 599 (a, -a,) M] no1 >dyn/ Cm2 


Tl-T2 


[1.4106208  (os-Oj)iT  -  4.64325004  (0^0  ,)  AT]  xlo 1 ’dyn/cm 


We  note  in  Eqn.  (51)  that  the  thermal  stresses  would  be 
anisotropic  even  in  the  limit  a^o^,  in  which  the  thermal  expansion 
is  isotropic;  this  is  in 'contrast  to  the  (OOl)Si  case  in  which  the 
stress  anisotropy  depends  solely  and  directly  on  the  thermal  expansion 
anisotropy . 

The  stress  equations  evaluated  for  Case  1  and  Case  2,  our 
two  sets  of  thermal  expansion  data,  and  for  T=1100C  deg  yield 

Case  1 

T  (Tl  +  T2)  =  -  1 . 1 567 1 2x10 1 °dyn/ cm2 , 

\  f^i"^2^  =  “0 • 3 16 5467xl09dyn/ cm2 , 

T6  =  °-° 

(52) 

Case  2 

j  (T1  +  T2)  =  -1 .  149017xl010dyn/cm2 , 

\  (T1-T2)  =  -0. 1047437xl09dyn/cm2, 

T6  "  °-°- 

The  average  stress  (T^+T2)/2  is  approximately  the  same  in  the 
two  cases.  However,  (T  .-T^)  /2  decreases  substantially  for  the 
larger  (ot^-c^)  case  and  reflects  competing  origins  for  the  stress 
anisotropy,  as  is  apparent  from  Eqns.  (51). 

Substitution  of  the  above  stresses  yields  the  data  for 
A p  1  / PQ »  E2^PoJ1  and  E3^PoJ1  shown  in  Table  V.  The  remaining 
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piezoresistance  effects,  Z'2/poJ3  and  Ap3/po‘  are  also  ^on-zero  for 
this  (221)Si  orientation  and  can,  if  desired,  be  evaluated  using 
the  information  given. 

One  general  characteristic  of  the  data  of  Table  V  compared 
to  the  data  for  (OOl)Si/ C 0 1 1 2 ) A 5, 2 0 3  should  be  noted.  For  (OOl)Si, 
the  anisotropy  in  mobility  and  the  transverse  effect  E ^/P  J'  both 
depend  directly  on  the  anisotropy  in  thermal  expansion  and 

are  thus  relatively  different  for  Case  1  and  Case  2  expansion  data. 
However,  for  ( 22 1 ) S i  an  important  contribution  to  the  piezoresistance 
effects  is  independent  of  the  anisotropy  in  thermal  expansion. 

The  data  in  Table  V  thus  show  little  relative  change  between  Case  1 
and  Case  2  expansion  data. 

We  next  consider  the  transverse  effect  E'/p  j'  data  in 

3  o  1 

Table  V.  This  effect  depends  principally  upon  the  term  (tt  ,-tt  .it  ■) 

111244'* 

which  is  numerically  similar  for  n-  and  p-Si.  Thus  all  E:/p  J' 

3  o  1 

data  in  Table  V  are  substantially  identical.  This  effect  involves 
the  generation  of  a  field  E normal  to  the  Si  film  plane,  induced 
by  the  current  in  the  plane  of  the  film.  Whether  this  field  is 
of  practical  significance  is  not  clear  and  no  experimental  studies 
have  been  performed.  In  principle,  if  the  zero-stress  resistivity 
and  mobility  are  functions  of  z-cocrdinate  position  (depth)  into 
the  film  --  as  is  always  the  case  to  some  degree  in  thin  hetero- 
epitaxial  films  then  the  field  generated  could  be  non-uniform 
in  the  z-direction  and  act  to  produce  a  modified  charge  carrier 
distribution  in  the  film. 

We  next  consider  the  transverse  effect  E'/p  j'  Again,  all 

£  ox 

data  in  Table  V  are  similar  for  the  reasons  given  above.  For 


n-type  Si  the  coefficient  of  the  sin20  term  is  +0,21227  for  Case  1 
and  +0.2314387  for  Case  2.  These  results  are  in  good  agreement 
with  those  of  Ref.  (2)  in  which  a  least-squares  fit  to  experimental 
data  yields  a  coefficient  of  +0.2. 


We  next  consider  the  longitudinal  piezoresistance  Ap'/p 

1  o 

and  the  associated  mobility.  Using  u  =  RH/p,  where  RH  is  the  Hall 
constant,  the  mobility  y/yQ-l (1+Ap J/Pq) .  The  characterizations  of 
the  mobility  m  terms  of  dmax,  ymin,  MA  and  A  are  given  in  Table  V. 
Mobility  data  derived  from  the  exponential  order  non-linear  model 
are  also  shown  in  Table  V  and  denoted  by  ue  ,  ue .  ue  and  Ae 
This  notation  will  be  used  also  in  subsequent  discussion  and  employs 
the  superscript  e  to  denote  the  use  of  the  exponential  order  non¬ 
linear  mobility  formulae  of  Eqn.  (33). 

The  general  conclusion  of  the  theory  is  that  for  p-type  (221) Si 
a  substantial  mobility  enhancement  should  result  from  thermal 
stresses. and  mobility  anisotropies  approaching  50%  should  be  obtained. 
However,  no  experimental  data  is  available  and  comparison  of  theory 
with  experiment  cannot  be  made. 

For  n-type  ( 2 2 1 ) S i  a  mobility  maximum  occurs  along  the  [ 1 T 0 ] S i 
direction  and  a  mobility  minimum  occurs  90  deg  away  along  the  [ 1 1 4 ] S i 
direction.  The  theory  predicts  that  at  an  angle  of  -36.6  deg  from 
the  [  110] Si  direction  the  effect  of  thermal  stress  on  the  mobility 
is  zero,  so  that  p(36.6  deg)=yQ,  the  unstressed  mobility.  The 
difference  between  the  bulk  mobility  and  the  value  of  mobility 
measured  at  this  angle  thus  represents  the  mobility  degradation  due 
to  causes  othex1  than  thermally  induced  stress. 
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Experimental  data  for  n-type  (221)Si  mobility  measurements 
are  given  in  Ref.  (2)  and  indicate  an  average  value  of  y  the 
unstressed  mobility,  of  600  cm2/V-sec  and  an  average  anisotropy  of 
A=39%.  The  fact  that  yQ  is  so  much  smaller  than  the  bulk  mobility^27-* 
strongly  suggests  that  thermally  induced  stress  is  not  the  dominant 
mechanism  in  lowering  the  mobility  in  epitaxial  films.  The  average 
experimental  anisotropy  of  39%  is  in  good  agreement  with  that 
predicted  theoretically  and  lends  confidence  to  the  theoretical 
model  employed. 

The  enhanced  mobility  predicted  for  some  current  direction 
suggests  that  the  (221)  plane  of  n-type  Si/Aii^Og  may  yield  superior 
electron  mobilities.  In  comparison  with  (001)Si  and  (lll)Si  grown 
at  this  laboratory,  the  (221)  plane  has  indeed  yielded  higher 
mob ilities^2^ . 

C.  (lllJSi/A^O^Mode  II,  III,  IV  Epitaxy) 

There  are  at  least  three  substrate  orientations  and  different 
modes  of  epitaxy  yielding  (11 1 ) Si  growth.  In  terms  of  our  model 
and  piezoresistance  calculations,  these  different  modes  differ 
only  in  the  substrate- induced  thermal  strains  and  can  therefore  be 
treated  all  at  the  same  time. 


Following  tie  procedure  illustrated  in  the  (noi)Si  and  ( 2  2 1 )  S i 
applications,  the  (£mn)"  matrix  for  (lll)Si  is 


(£mn)=  — 

/6 


{fb  cosG  +sin8 
cos 0  -/3  sine 

V/2 


- /3  cos0  +sin0 
cos  0  +  /3  s in0 
fl 


-  2  s  i  n  0 

-  2  c  o  s  0 


/2  / 


(5  3) 
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where  the  angle  9  is  measured  from  the  [  110] Si  direction  toward 
the  [112]Si  direction.  The  piezoresistance  equations,  for  (lll)Si 
and  explicitly  involving  the  angle  0,  become 


•  VT2,  C2’llt4”l2*"44)  ,VT2,  flfir"l2*2"44) 

o  J  5  +  l - 1  J - - 


APi/P0  -  - —  3--  ZJ 


cos  2  9 


Oin-w12»2»44) 

+  T  - = -  s  m2  9 , 


(VV 
e2/poj{  =  t — h1- 


TT.  .  -7T  _  +  2  7T 

sin20  +T  cos20]  (— — - — - — )  , 


T  -T 
1  2 


/2 


E3/PoJ1  =  t  C - 2 - 5  Sin6  +T6cos0]  —  (1Til-ni2-TT44)  * 


(54) 


e;/poj3  =  E3/p0j;  * 


T  -  T 

E2/poJ5  =  cos6  -T6sin0]^|  C11-'’12-"44)  . 


Ap'/p0  =  ( 


Tl  +  T2  /2 

2  (TT11  +  2TT12-1T44) 


(28) 


The  stress  equations  (Eqn.  (22))  are  then  solved  next,  subject 
again  to  the  boundary  condition  T3=T4=T5=0.  The  form  of  the 
elastic  constant  matrix  is  similar  to  that  for  the  ( 2 2 1 ) S i  case 
but  numerically  simpler.  The  reaction  strains  e'  and  e'  are 
non-zero.  e  ^  is  non-zero  unless  c'  is  zero,  as  will  be  seen  to 
be  the  case  for  some  substrate  orientations. 

The  results  for  the  stresses  T^,  and  T  ,  given  in  terms 
of  N  ^  and  N  ^  relevant  to  the  particular  substrate  orientation,  are 
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Data  for  (lll)Si  Mode  II,  III,  IV  epitaxy  are  given  in 
Tabla  VI,  The  values  of  the  substrate  transformation  direction 
cosines  and  are  shown,  and  the  resulting  stresses  calculated 
for  Case  1  and  Case  2  expansion  data  are  also  given. 

In  Table  VI  only  the  longitudinal  piezoresistance  effect 
Apj/pQ  has  been  considered  and  the  mobility  derived  from  this  is 
shown.  Because  the  anisotropy  formulae  ail  involve  functions  of 
cos20and  sin20  there  will  be  four  mobility  extrema,  with  two  maxima 
and  two  minima.  Adjacent  maxima  and  minima  are  separated  by  90  deg. 
For  n-type  Si  the  mobility  maxima  for  Modes  II  and  IV  occur  at  0=0 
which  is  along  the  [1100]  and  [  1  2  10 ]  A ^  directions,  respectively. 
For  Mode  III,  a  maximum  mobility  occurs  at  0  =  38.242  deg,  which  is 
also  along  the  [llOOjAX^O^  direction. 

The  [ 1 1 0  0 ] ,  [1210]  and  [110  0]  directions  for  Modes  II,  III, 

IV,  respectively,  all  are  perpendicular  to  the  C-axis  or  [0001] 
direction.  These  directions  thus  pertain  to  directions  yielding 
the  minimum  thermal  expansion  lying  within  the  plane  of  the 

film.  Clearly  the  mobilities  in  these  directions  should  represent 


minima  locations  in  the  plane  are  interchanged  for  p-type  Si  as 
compared  with  n-type  Si. 

The  effects  of  stress  on  mobility  for  n-type  (lll)Si  are 
small,  as  indicated  in  Table  VI.  The  average  mobilities  are  close 
to  the  unstressed  mobility,  and  the  anisotropies  range  from  a  low 
of  2.46%  (Case  1,  Mode  IV)  to  a  high  of  9.63%  (Case  2,  Mode  III). 
The  Case  2,  Mode  1 1 1,  anisotropy  is,  however,  high  enough  that  it 
should  be  considered  in  employing  the  mobility  as  a  means  of 
evaluating  the  effects  of  change  in  growth  parameters. 

The  effects  of  thermal  stress  on  p-type  (lll)Si  are  much 
more  significant  than  for  n-type  (lll)Si.  A  substantial  overall 
mobility  enhancement  approaching  a  factor  of  two  is  predicted. 

Based  on  the  (OOl)Si  and  (221) Si  data,  Case  2  thermal  expansion 
data  would  be  most  applicable.  For  this  case,  the  anisotropy 
factors  A  range  from  6%  up  to  15%  and  thus  could  be  of  experimental 
significance  in  measurements  of  thin  film  mobilities. 
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VI.  APPLICATION  TO  AN  OFF-PLANE  (lll)Si  ORIENTATION 

In  the  previous  section  the  piezoresistance  and  mobility 

anisotropy  were  calculated  for  Mode  III  epitaxy  with  a  (1  l3o)  AJl  203 

substrate  orientation.  An  examination  of  Fig.  4  of  Ref.  (5), 

showing  a  stereographic  projection  illustrating  Si/A£203  epitaxial 

relations,  indicates  that  a  (ll20)A£203  substrate  orientation  is 

either  on  or  close  to  the  boundary  between  Mode  II  and  Mode  III 

epitaxy.  Therefore,  in  order  to  insure  good  quality  Si  film  growths 
of  only  one  orientation 

the  A£203  substrate  is  often  intentionally  cut  a  few  degrees  off 
the  (1120)  orientation^29-^  .  In  this  section  we  shall  apply  the 
theory  developed  in  this  paper  to  a  substrate  orientation  which 
is  *5  deg  off  the  (1120)  plane.  This  application  illustrates  some 
of  the  formalism  developed  in  Section  II  and  also  provides  infor¬ 
mation  on  changes  in  mobility  to  be  expected  for  slight  off-plane 
Si/A£203  film  growths. 

The  question  of  equivalent  orientations  and  internal  mathe¬ 
matical  consistency  was  briefly  discussed  in  Section  II.  This 
problem  arises  naturally  in  the  characterization  of  some  off-plane 
-  ( 1 1 20 ) A ^2 0 3  substrates  used  in  this  laboratory  to  obtain  good 
quality  555  (lll)Si  film  growths. 

The  determinations  cf  the  substrate  surface  orientation  and 
of  the  film/substrate  parallel  relations  both  involve  the  identi¬ 
fication  and  labeling  of  certain  principal  x-ray  reflections  from 
the  Si  film  (e.g.,  {311})  and  from  the  A^2°3  substrate  (e.g.,  {1120}) 
The  ambiguity  in  labeling  these  r e f 1 e c t i ons - - 1 h ree  choices  each  for 
Si  and  A£203--leads  to  different  but  equivalent  coordinate  systems 
and  characterizations  of  substrate  surface  and  parallel  relations. 
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The  question  of  mathematical  consistency  comes  in  when  sub¬ 
strate  surface  orientation  characterization  and  Je t erminat ion  of  the 
parallel  relations  are  not  made  with  an  identical  la  eling  of  x-ray 
reflections.  In  this  case,  care  must  be  taken  to  insure  that 
the  pole  diagrams  selected  to  characterize  the  film  and  substrate 
surfaces  are  consistent  with  the  parallel  relations  assumed  to 
hold  for  the  f i lm/ s ub s t r a t e  composite. 

Pole  diagrams  for  one  substrate  o r i ent at i on* o f  interest 


are  shown  in  Fig.  1.  Figures  la,  lb,  lc  are  equivalent  in 
terms  of  describing  the  film  and  substrate  surfaces,  but  as  shown 
relate  to  different  parallel  relations.  The  orientation 
shown  as  Fig.  la  applies  to  the  parallel  relations  previously 
adopted  in  Section  II  and  will  be  employed  in  the  subsequent 
calculations.  The  substrate  orientation  indicated  in  Fig.  m  was 
cut  +5  deg  from  the  (1120)  plane  toward  the  (0110)  plane  and  along 
the  zone  indicated.  The  axis  of  rotation  for  the  cut  is  the  [0001] 
direction.  For  this  mode  of  epitaxy,  the  (lll)Si  plane  is  ||  to 
the  (1120)A£203  plane.  However,  the  relevant  crystallographic 
zones  in  AJ^Oj  and  Si  are  not  exactly  parallel  so  that  the  Si 
axis  of  rotation  is  only  approximately  the  [Oil]  direction.  The 
pole  diagram  shown  for  Si  is  thus  only  approximate  and  will  not 

be  used  to  generate  the  (£mn)  matrix  for  Si  required  to  calculate 
piezoresistance  effects. 


The  rotation  in  A*^  of  Fig.  la  is  about  the  hexagonal 
c-axis  which  is  the  Cartesian  z-axis.  If  (0)  is  the  rotation 
matrix  acting  upon  a  column  vector,  and  which,  for  0=+3O  deg, 
rotates  the  Cartesian  vector  related  to  [1120]  into  the  Cartesian 
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related  to  [OlrO],  then  the  (LMNf'.atrix  for  this  off-plane 
substrate  is  given  by  (LMN)'=  (LHN)  (6),  where  (B)  is  the  transpose 
of  CB)  and  ( LMN)  relates  to  the  "on-plane"  (1120)  case.  Since 
the  rotation  is  about  the  z-axis,  the  direction  cosines  Nj  and  N., 

Which  define  the  thermal  strains  are  constant  and  independent  of2 
the  angle  of  rotation 

From  Section  II,  the  («,„)' matrix  associated  with  the  off- 

P 13116  Si  fil”  is  8ive"  ^  (*■")'•  (LMN)'(SF)  .  CLMN)C?5(SF),  Where 
(SF)  is  taken  from  Table  II.  For  (limn)' we  obtain  (e  =  *5  deg) 


fll/J  °  0  \  Al.0S2475  -0.9128301  *0. 0580268 

0  1//?  0  I  ♦1.0716329  *1.078041  -  1.920788 


\°  0  17  /  \+0-8742111  041  88  +  1.0724898/ 

where  the  normalization  factors  used  wore  chosen  to  facilitate 

comparison  with  the  on-plane  simple  (1 1 1 )  St  /  (1  ljo)  Al^  case. 

From  the  third  row  of  E,„ .  (56)  we  see  that  the  Si  film  surface 

lies  approximately  along  the  zone  between  the  (111)  and  (011)Si 
planes. 

The  above  discussion  relates  to  an  off-plane  substrate 
orientation  which  is  *5  deg  away  fro.  the  (1120)  plane.  Calcu¬ 
lations  have  also  been  performed  for  a  -S  deg  rotation  yielding 
plane  substrate  surface  which  lies  between  the  (1120)  and 
(1010)  planes.  Although  perhaps  of  less  experimental  interest, 
this  latter  -5  deg  rotation  data  coupled  with  the  on-plane  data 
and  the  *5  deg  rotation  data  provide  insight  into  the  sensitivity 
of  piezoresist. nee  and  mobility  to  small  changes  in  substrate 
orientation  for  this  mode  of  (lll)Si  epitaxy. 


The  formulas  for  the  longitudinal  piezoresistance  for 
0=  +5  deg  are  rather  lengthy  and  will  not  be  given.  We  note  only 
that  there  are  contributions  to  the  anisotropy  in  piezoresistance 
which  do  not  vanish  in  the  limit,  of  stress  isotropy.  This 

result  is  similar  to  that  found  for  the  (221)Si  films,  as  one  would 
expect. 

The  stresses  T.,T0  and  T,  are  numerically  different  for 

I  4  O 

these  various  orientations,  even  though  the  direction  cosines 
N  ^  and  N,  determining  the  thermal  strains  are  independent  of 
rotation  angle,  since  the  Si  elastic  constant  matrix  is  itself 
different  for  the  different  orientations.  Numerical  results  for 
the  (11203+5  deg,  (1120),  and  (1120)-5  deg  orientations  are  shown 
in  Table  VI.  The  average  stress  (T^+T^)  is  nearly  constant. 

However  (T-T23/2  and  T&  vary  appreciably. 

Numerical  data  for  the  mobilities  derived  from  the  longi¬ 
tudinal  pi e zo re s i s t ance  Ap^/MQ  are  given  in  Table  VII.  Values 
ft v  u  /u  ,  u  .  /y  ,6  --the  angle  from  the  reference  axis  at 

max  o’  ">m  0  W)nax 

which  the  mobility  is  a  maximum--and  the  anisotropy  factor  A  are 

given  for  both  linear  and  exponential  models.  Since  the  formulae 

evaluated  involve  cos29  and  sin28  terms,  adjacent  maxima  and  minima 

are  separated  by  90  deg  and  equivalent  maxima  and  minima  are 

separated  by  180  deg.  Exact  comparisons  of  the  angles  0  are 

Mmax 

somewhat  (less  than  +5  deg)  ambiguous  since  the  reference  axis  is 
different  in  each  case. 

The  results  for  both  n  and  p-type  Si  are  different  in  the 
various  cases  and  reflect  competing  phenomena  in  establishing  the 
total  anisotropy.  As  expected,  the  effects  of  stress  on  n-type 
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Si  Hall  mobilities  are  much  less  than  for  p-type  Si.  The  maximum 
mobility  hmax  varies  by  +7.5%  and  the  anisotropy  Ae  changes  by  *15% 
for  p-type  Si  and  a  +5  deg  change  in  substrate  orientation.  Thus 
the  mobility,  based  on  thermal  expansion  stresses,  would  be  expected 
to  be  a  sensitive  function  of  substrate  orientation. 

The  dependence  of  p-type  Si  mobilities  on  orientation  suggests 
that  proper  selection  of  substrate  orientation  can  lead  to  enhanced 
mobilities.  Of  course,  physical  phenomena  other  than  thermal  stress, 

may  also  play  an  impor¬ 
tant  role  and  have  to  be  considered  in  selection  of  optimum  sub¬ 
strate  and  film  orientations. 
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VII.  APPLICATION  TO  THE  GENERAL  (XXI)  Si  ORIENTATION 

The  applications  of  the  thermal  s t r es s - p i e z ore s i s t an c e  model 
have  thus  far  been  limited  primarily  to  film  and  substrate 
orientations  which  are  of  current  experimental  interest.  In  this 
section  we  shall  apply  the  model  in  a  predictive  sense  to  general 
Si  film  orientations  of  the  form  (xxl)  to  determine  the  effect  of 
thermal  stress  on  the  mobility  for  other  Si  orientations  which 

have  not  been  researched  experimentally. 

* 

Workers  in  the  field  of  Si  epitaxy  have  long  recognized  that 
Si  film  quality  depends  upon  a  number  of  physical 

phenomena.  The  use  of  thermal  stress  model  mobility  data  to 
determine  Si  orientations  of  possible  experimental  interest  thus 
focases  attention  on  only  one  part  of  the  epitaxy  problem.  None¬ 
theless,  the  effect  of  thermal  expansion  stress  on  Si  film 

,  .  .  .  .  the 

mobilities  is  substantial  and  the  use  of/thermal  stress  model  to 

predict  Si  film  orientations  of  potential  experimental  interest 
is  an  important  and  useful  new  theoretical  tool. 

Extensive  computer  calculations  for  the  (xxl)Si  film  orien¬ 
tation  have  been  performed.  By  varying  x  from  zero  to  infinity, 

several  Si  film  orientations  of  the  form  (xxl)  lying  along  the 

( 2  2  1 ) Si 

zone  (s )  defined  by  (OOl)Si,  (Ul)Si,/and  (110)Si  were  treated  and 
numerical  data  obtained.  Since  Si  film  quality  is  usually  highest 
for  (001)  and  (111)  types  of  growth,  this  general  treatment  is 
expected  to  encompass  the  orientations  of  major  interest. 

In  the  theoretical  treatment,  the  four  modes  of  Si/A£203 
epitaxy  and  the  mode  of  Si/MgAS^O^  epitaxy  defined  in  Section  II 
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are  all  included.  The  procedure  is  to  assume 

for  purposes  of  the  calculation  that  each  given  mode  exists 
for  all  Si  orientations  lying  between  (OOl)Si  and  (110)Si.  In 
reality  this  is  not  true  and  boundaries  between  the  various  modes 
must  be  considered.  In  general,  in  proceeding  away  from  the 
"natural"  plane  for  a  given  mode--say  (001) Si// (0112) AAO  (Mode  I)-- 
the  quality  of  the  film  deteriorates  and  the  Si  film  may  become 
mu  1 1 imoded  or  defected.  An  adequate  number  of  growth  parameter 
studies  have, however,  not  been  directed  specifically  toward  this 
question.  Thus,  the  exact  relationship  between  film  quality  and 
the  number  of  degrees  off  of  the  orientation  of  a  pure  mode  is  not 
known  in  detail,  although  some  general  information  has  been  estab¬ 
lished  by  earlier  workers(6).  The  (22 1 ) Si/ / ( 1 1 2 2) A£2 0  is  an 
example  of  an  orientation  that  is  =16  deg  off  the  "natural" 

(lll)Si// (1124)A£203  plane  and  yet  yields  good  quality  films^2\ 

The  reference  axes  for  all  of  the  (xxl)Si  plane  calculations 
are  assumed  to  be  X'-[lIo],  Y'-[l,l,2x]  and  Z|-[xxl].  The  (imn) 
matrix  for  this  choice  of  axes  is  then 


(£mn) 


/2  /2 


where  b=(l  +  2x2)2.  The  substrate  surface  orientations  associated 
with  the  (xxl)Si  plane  can  be  obtained  from  the  row  vector  (L3.M3.N3) 
in  the  (LMN)  matrix  by  using  Eqn.  (6)  and  the  (SF)  transformation 
matrices  given  in  Section  II.  The  ability  to  mathematically  treat 
general  Si/Ai^O-j  orientations  of  the  above  and  other  arbitrary 
forms  is,  of  course,  the  main  reason  a  general  mathematical 


description  of  Si  epitaxy  was  sought,  developed,  and  presented  in 
Section  II. 

The  mobility  has  been  determined  from  evaluation  of  the 
longitudinal  piezoresistance  aP'/Pq  for  four  modes  of  Si/A*20 
epitaxy  and  one  mode  of  Si/MgAi^  epitaxy.  Calculations  were 
performed  for  twelve  values  of  x  in  Cxxl)Si  ranging  from  (OOl)Si) 
through  ( 1 10) Si ,  and  both  n  and  p-type  Si  were  treated.  For  the 

Sl/,A2,2°3  modes>  Case  2  thermal  expansion  data  [ Ct1  ~  a2  ~  ~  1  •  0  8x  1 0  "  6  /  deg  C) 
was  employed.  For  the  case  of  MgA£204  an  average  [isotropic) 
expansion  coefficient  of  8.1xl0_6/deg  C  was  employed C 10) .  A 

temperature  difference  of  AT-1100C  deg  was  assumed  for  both  cases. 

The  numerical  data  obtained  are  given  in  Tables  I X -  X I  1 1 . 

In  the  tables,  the  first  column  lists  the  value  of  x  in  (xxl)Si. 

The  second  column  shows  the  number  of  degrees  the  subject  plane  is 
from  the  COO  1 ) Si  plane  and  toward  the  [110)Si  plane.  The  remaining 
columns  list  the  maximum  mobility  ratio  u  /u  ,  the  minimum 

Ml  a.  A  V 

mobility  ratio  lJmin/d0,  and  the  anisotropy  parameter  A.  For  p-type 
Si,  mobilities  calculated  from  the  linear  and  exponential  models 
differ  markedly.  Three  additional  columns  have  therefore  been 
added  to  Tables  IX-XIII  for  U®ax/uo,  y®in/yQ  and  Ae  data  based  upon 
the  exponential  non-linear  correction  term.  However,  in  applications 
and  discussion,  only  the  non-linear  data  will  be  employed. 

Before  discussing  the  numerical  data,  it  may  be  useful  to 
again  comment  briefly  on  the  various  origins  of  the  anisotropy  in 
mobility.  Basically  it  is  the  anisotropy  [in  the  plane  of  the  Si 
film)  of  the  piezoresistance  coefficients,  of  the  elastic  constants 
and  of  the  s ub s trate- indu ced  thermal  strain  that  leads 
tropies  in  carrier  mobilities. 
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Kith  respect  to  the  Si  fil.,  the  carrier  mobility  anisotropy 
originating  in  the  Si  piezoresistance  coefficients  and  in  the  elastic 
constants  can  be  considered  as  intrinsic  to  the  Si  and  the  mobility 
anisotropy  originating  in  substrate  thermal  expansion  anisotropy  as 
extrinsic  to  the  Si.  In  general,  the  calculated  mobility  anisotropy 
will  represent  an  admixture  of  intrinsic  and  extrinsic  effects  and 
no  useful  subdivision  is  possible. 


For  Si/MgAt20,  (Mode  V)  the  substrate  thermal  expansion  is 
isotropic  and  thus  extrinsic  mobility  anisotropy  is  zero.  The  Si 
piezoresistance  coefficients  and  elastic  constants  are  isotropic 
for  the  cases  of  (OOl)Si  and  (lll)Si.  Thus  the  intrinsic  mobility 
anisotropy  is  also  zero  for  these  two  orientations. 


For  the  Si/Af203  modes  of  epitaxy,  the  calculated  anisotropy 
for  (00 1 ) Si  and  (lll)si  is  extrinsic  and  smaller  than  the  intrinsic 
anisotropy  for  most  other  orientations.  Thus  the  anisotropy 


parameter  A  as  a  function  of 
local  minima  for  the  (001) Si 


orientation  would  l  expected  to  exhibit 
and  ( 1 1 1 ) Si  orientations.  The  data  of 


Tables  I X  -  X 1 1 1  illustrate  this  point  with  local 
at  or  near  these  two  orientations. 


minima  in  A  occurring 


The  thermal  stress  mobility  data  in  Tables  IX-XIII  are  extensive 
enough  to  permit  a  systematic  investigation  of  theoretical  mobility 
as  a  function  of  orientation  in  the  defined  zone.  We  first  discuss 

those  characteristics  which  are  qualitatively  similar  for  n-  and 
p-type  Si  films. 

The  maximum  mobility  datn  n  ,,e  , 

y  aata-  “max^o  °r  Umax/lV  sho"  that  both 

and  p-type  mob  lilies  are  minimum  for  (001)Si  and  increase  as 
one  proceeds  across. the  defined  zone  toward  (110)Si.  A  local 
minimum  in  the  maximum  mobility  is  evident  near  (lll)Si.  The 
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theoretical  maximum  mobilities  for  the  (001),  (lU)>  and  (221)si 
orientations  thus  can  be  ranked  with  (00  1 ,  <0^  ( 11 ,  ,  <„  (22  ,  , 

FU”S  8r°“"  in  thiS  laboratory  appear,  on  the  average,  to  confirm 
this  prediction*-1,2^. 

The  minimum  mobilities,  or  y^/^,  ^  ^ 

tatively  similar  for  n-  a,  d  p-type  films  and  exhibit  minima  for 

(001)  and  (110)Si  and  maxima  near  (lll)Si.  The  fact  that  the 

maxima  in  U  .  /u  or  ue  /u  ,  .  . 

mm  o  Mmin/Uo  anc*  t^1e  minima  in  y  /y  or  ye  / 

max  0  max/Mo 

mobilities  are  near  n i  n  _  c 

(lll)Si  is,  of  course,  due  to  the  minima  in 

the  anisotropy  parameters  A  or  A*  which  occur  near  (lll)Si. 

Although  the  n-  and  p-type  mobility  data  do  exhibit  quali¬ 
tative  similarities,  there  are  marked  quantitative  differences. 

Lev  us  consider  the  n-type  Si  mobility  first.  The  maximum  mobility 
“max  ran*es  from  a»30%  reduction  in  mobility  near  (001)Si  up  to  a 
=  25-1  enhancement  in  mobility  near  (llo)Si,  as  compared  with  the 
zero-stress  value  UQ .  The  minimum  mobility  is  always  less 

than  uo  by  amounts  ranging  from  *30-40%  near  (001)  and  (llO)Si 
down  to  perhaps  5%  near  (lll)Si. 

For  p-type  Si  films,  the  effect  of  stress  is  always  to 

increase  the  mobility,  and  both  and  are  greater  than  „ 

across  the  entire  tone.  The  stress  enhancement  of  the  maximum 

-obilitv  is  minimal  for  (001)Si  and  increases  rapidly  in  going 

across  the  tone.  Values  of  the  maximum  mobility  ratio  ve  /u 
„ .  ,  max  o 

in  the  range  1.06-1.09  for  (001)Si,  1.80-1.96  for  (111, Si, 

2-10-2.27  for  (221, Si  and  2.32-2.44  for  (llo)Si.  Thus,  as  compared 

to  (001)Si,  substantially  larger  theoretical  mobilities  can  be 

achieved  by  choice  of  f i 1 m/ sub s t ra t e  orientations  lying  within 
the  zone . . 
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Iho  am  so, ropy  parameter  Ae  for  p-,ype  Si  films  is  als0  a 
s  t  rong  f unc t  i  on  of  orientation  and  is  in  the  range  0-S%  for  (001)Si, 
lot  (llltSi,  4  S  -  5  5 1  for  (22])si,  and  74-  80%  for  (llfl)Si. 

rrCdiC,Cd  “"i%«ropy  along  most  of  the  tone  is  large  enough 
that  orientation  direction  within  the  plane  of  the  film  should  be 
taken  into  account  in  experimental  mobility  measurements  to 
evaluate  the  effect  of  changes  in  film  growth  parameters.  The 
genera]  characteristics  of  both  and  p-,ype  Si  are  illustrated 

Figs-  (2)  and  (3)  where  th  shaded  regions  indicate  the  range 
Ot  mobilities  associated  with  Modes  1-iv. 

The  theoretical  mobility  data  presented  in  Tables  1X-XII1 

^  fi8S-  2  °"d  3  thC  "suits  of  a  "mode ! "  and  the  assumptions 
of  this  model  must  he  borne  in  mind  in  employing  this  data  to 
estimate  experimental  mobilities. 

■ho  first  point  concerns  the  assumptions  underlying  the 
(OOlJ-CMOjSi  "continuous. tone"  calculation.  Mathematically  we 
have  assumed  that  all  modes  of  epitaxy  hold  across  the  entire  tone 
Realistically,  h  owe  ve  r  —  s  i  g„  i  f  i  cance  should  he  attached 
-'ability  data  for  points  in  the  tone  which  are  far  removed  from 
-na,  regions  in  which  the  modal  relations  have  been  shown  experi- 
mentally  to  hold  and  good  quality  films  obtained. 

The  basic  assumptions  of  the  thermal  s t re s s - p i e t o r es i s t ance 
model  already  have  been  discussed.  For  the  case  of  „.type  Si, 
comparisons  between  theory  and  experiment  for  the  mobility  anisotropy 
transverse  fteid  effect  have  been  made  for  (001)  and  (221)Si 
orientations.  Substantial  agreement  was  obtained,  thus  lending 
confidence  to  the  thermal  stress  mode, <‘.2).  This  suggests  that 
the  continuous-zone  theoretical  predictions  for  „-type  Si  should 


be  useful  for  zonal  points  in  the  vicinity  of  the  (001)  and  (221) 
orientations  for  the  appropriate  modes  of  epitaxy. 

For  p-type  Si,  no  mobility  anisotropy  measurements  have  as 
yet  been  performed  and  compare<  with  the  stress  model  calculations. 
Nonetheless,  we  believe  the  calculated  effects  of  stress  on  p-type 
mobilities  are  relevant  and  significant  and  that  the  zonal  calcu¬ 
lations  are  useful  in  a  predictive  sense  to  indicate  new  orientations 
worthy  of  experimental  study. 
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VIII.  DISCUSSION  AND  CONCLUSIONS 

Ihe  principal  objective  of  the  present  investigation  has  been 
the  theoretical  determination  of  the  effects  of  th er ma 1  - rxpan s i on  - 
induced  stress  on  the  carrier  mobility  in  n-  and  p-type  Si/Ai  0 
films.  2  3 

A  mathe¬ 
matical  description  of  he t e roep i t ax i a  1  relations  was  developed  in 

Section  II,  and  while  original,  is  mathematically  straightforward. 
Ihe  main  utility  of  the  formalism  is  in  mathematically  treating 
"off-orientation"  Si/AA  0  films. 

A  model  for  calculating  t he r ma 1 - e xpans i on- i nduce d  anisotropic 
stresses  was  developed  in  Section  III.  Within  the  framework  of  a 
thin-film-thick-substrate  approach,  the  model  rigorously  includes 
and  treats  the  anisotropy  in  thermal  expansion  coefficients  and 
the  ui. isotropy  in  the  Si  stresses. 

A  few  points  concerning  the  thermal  stresses  calculated  for 

the  Si/A£203  orientations  treated  in  the  paper  should  be  mentioned. 

The  average  stress,  (Tj+T2)/2,  ranged  from  ~-0.83xl010  dyn/cm2  up 
10  2 

to  --1.20x10  dyn/cm  ,  depending  upon  mode  of  epitaxy  and  orien¬ 
tation  within  the  mode.  These  values  are  consistent  with  available 
experimental  estimates*-4*5^. 

The  second  point  concerns  the  existence  of  shear  stresses  in 
the  S i  films.  Apparently  it  is  generally  assumed  in  the  literature 
and  by  workers  in  the  field  that  there  are  no  shear  stresses  in 
thin  S i / A Z 2 0 3  filmst30).  On  the  other  hand,  a  non-zero  shear 
stress  Tg  in  the  plane  of  the  film  arises  naturally  in  our  aniso¬ 
tropic  stress  model  and  will  be  discussed  briefly. 
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In  the  thin-film-thick-substrate  model,  the  Si  fil.  street 

is  characterized  by  T, ,  T2  and  T6  evaluated  for  so,,  choice  of 

reference  axes  with  x  and  y  in  the  plane  of  the  film.  It  turns 

out  that  the  existence  of  the  shear  stmcc  t 

near  stress  T&  iS  a  consequence  of 

the  stress  anisotropy  and  will  he  present  unless  T^T  The 

argument  is  as  follows  u  n  . 

'  Ut  the  stress  ^rms  (TrT2)/2  and  T6  be 

defined  tor  the  reference  axes  Than 

rence  axes.  Then  rotate  the  coordinate  system 

by  45  deg.  Ftc,  E,„.  (36),  and  T'=  -  (T , -Tj )  / 2  and  the 

two  stress  terms  have  interchanged.  Thus,  nyxp/2  and  T,  are 

physically  equivalent  and  hence  it  is  only  for  the  isotropic  case, 

Tj=T2,  that  there  is  physically  no  shear  stress. 

For  the  calculations  performed,  the  values  of  (T|-T2)/2  and  T6 

ranged  from  zero  up  to  M0%  of  the  average  stress.  Values  equa,  t! 

5-7*.  of  the  average  stress  were  more  typical  for  most  orientations 

and  modes  of  epitaxy.  The  relative  magnitudes  of  both  (Tj-T.J/Z  and 

'‘6  a"  Str°"g  fUnCtlo"s  of  the  mode  and  of  the  film  orientations 
and  vary  appreciably. 

The  basic  piezoresistance  effect  was  discussed  in  Section  IV 
and  formulae  employed  in  the  paper  were  presented.  The  major 
emphasis  was  upon  the  linear  piezoresistance  model.  However,  a 
no„-li„ear  correction  proceedure  was  also  developed  which  attempts 
to  account  for  higher  order  stress  effects.  The  resulting  non¬ 


linear  model  was  employed  extensively  for  p.type  si 

data  for  which  the  linear  model  is  believed  to  overestimate  the 
effects  of  stress  on  the  resistivity. 

For  both  the  linear  and  non-linear  models,  it  was  further 
assumed  that  the  thin  Si  film  is  thick  enough  that  bulk  piezo¬ 
resistance  coefficients  and  bulk  coefficient  transformation  proper- 
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t.es  arc  applicable.  Thus  an,  device  situation  in  which  surface 
conduction  is  important  could  lie  outside  the  framework  of  the  model 
per  se,  and  would  have  to  he  examined  on  an  individual  basis. 

Ihe  Si  pietoresistance  coefficients  it.,  are  functions  of  the 
carrier  concentration.  Assuming  that  this  relationship  is  approxi- 
mately  the  same  in  Si/A^O,  and  in  Si  films  of  the  same  carrier 
cone en t ra t i on,  coe f f i c i ent  data  such  as  that  summarized  in  Ref.  [11) 
indicate  the  numerical  calculations  of  the  present  paper  are 
applicable  for  concentrations  up  to  .  1  0 1 7/cnT  5 .  For  higher  concen¬ 
trations  different  and,  with  the  exception  of  generally  smaller 

coefficients  would  he  applicable. 

Application  of  the  therma 1 -stress-pie zoresistance  model  to 
the  discrete  (001),  (221)  and  (lll)Si  orientations  was  made  in 

Section  V.  In  general, the  longitudinal  effect  ApJ/po  leads  to  a 
mobility  variation  of  the  form  ( U/ Uo )  =  [ a - b c o s 2 6 ] ‘ 1 ,  while  the 


transverse  effect  leads  to  a  field  variation  of  the  form  El'0  J>bsin20 

t  o  1 

The  zero-stress  mobility,  u  and  rpsistivftv  „  „  , 

Mo’  ani  resistivity,  pq,  are  not  know,'  and 

cannot  be  determined  directly  from  experimental  data.  They  are  thus 
essentially  derived  parameters  and  would  be  expected  to  be  a  function 
of  growth  conditions  and  to  vary  slightly  from  sample  to  sample. 

The  anisotropy  parameter  A  and  the  angular  orientations  in 
the  Si  plane  yielding  mobility  maxima  and  minima  are  independent 
of  and  thus  permit  direct  comparisons  between  theory  and  experiment. 
These  parameters,  rather  than  the  mobility,  per  se,  should  therefore 
be  emphasized.  Similarly,  the  orientations  of  maxima  and  minima  in 
the  transverse  field  are  basic  in  directly  correlating  the  trans¬ 
verse  field  data. 
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comparison  between  theory  and  experiment  for  rh  . 
and  the  location  of  the  ^  an,s°‘ropy  A 

C«USif2>  gave  good  ag  Varl<>Ui  ^  (001>SiU)  and 

mobility  an i sof  ^  argUGS  that  the 

effects  or  ph  ^  ^  ^  eXPUined  WitH°Ut  recourse  to  other 
P  ysical  phenomena.  However  thP 
cannot  entirely  role  agreement  obtained 

'  °Ut  °ther  Phenomena.  thls  vein 

occurrence  of  cosl,  and  sinlh  terms  in  the  • 

is  not  specific  t  Piez°resistance  effects 

peci fxc  to  any  one  physical  „h« 

Physical  phenomenon  but  instead  ie 

3  symmetry  requirement.  Slmply 

Measurements  of  mobility  anl.nt 

y  anis°tropy  and  transverse  fola  u 
also  been  ca-ri»H  .  averse  F  eld  have 

ca.rted  out  in  these  laboratories  for  a 
Clll)Si/C1120)A*  n  and  funs  /,,  -  "umber- of  n-type 

2  “  1  i  1 1 )  Si/ (lo  14  )  A£  0  f  1 1 T, 

anisotropy  factor  A  2  3  •  he  mobility 

PX  factor  A.  averaged  over  a  number  of  sa.ples  uas 

and  mIM  for  (III)11/(||1#  „ 

*"‘50tr0"  fa"°”  •»  substantially  greater  th  a  *  *' 

fr°"  the  thermal  stress  mode,  and  indicate  in  "  ‘  P"diC,'d 

'**  cases,  that  other  '  ""‘h’5'  *°  **  '°01)S‘ 

a  K“0,iuBiena  are  rennimj  • 

Application  of  the  thermal  stress  pie 
off-plane  .(,11, Si  re s s -pr e to  re s i s t ance  model  to 

onentations  was  made  in  Section  VI  Cal  i  • 

-re  performed  for  substrate  orientations  lying  „  „  '  °"S 

(11201  ni  •  y  ng  ±5  deg  from  the 

J  plane  in  a  general 

general  (non - symme t ry )  direction  tw 

results  indicate  that  cha  ’  **  numerical 

ord  "  SUbStrate  -ientation  of  this 

rder  can  produce  noticeable  changes  in  • 
changes  in  the  ■  •  anisotropy  and  substantial 

'  thS  l0Catl°"  °f  maxima  and  minima 

SeCti°n  VI1  aPPiicatio„  of  the  model  was  mad 
C.  .  e 1  was  made  to  genpmi 

31  “mentations  of  the  form  (xxl)  T(  .  a 

lAAij.  The  model  wa«  « _  i  . 

predictive  sense  to  H  .  •  employed  in- a 

ense  to  determine  the  effect  of 

”ect  of  thermal  stress  on  the 
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mo H 1  I ! ' y  for  other  Si  orientations  which  have  not  yet  been 
researched  experimentally.  In  these  calculations,  twelve  Si 
t  r i ent a t ion s  lying  along  the  zone(s)  defined  by  (OOl)Si,  (llllsi, 

(22 1 ) Si  and  (110)Si  were  treated  and  numerical  data  obtained  for 
the  four  common  modes  ot  Si/A*^  epitaxy  and  for  one  mode  of 
Si/N'gAi^O^  epitaxy. 

The  mobility  data  are  rather  similar  numerically  for  the 
various  modes  and  indicate  only  a  weak  dependence  upon  the 
particular  mode  of  epitaxy.  The  n-  and  p-,ype  maximum  mobility 

m  a  X  ^ 1  o  01  J  m  a  x  ^ 1J  o  ’  are  qualitatively  similar  in  that  the 
moh i  is  a  minimum  for  (OOl)Si  and  increases  as  one  proceeds 

across  th  defined  zone  toward  (HO)Si.  Be,h  and  p.type  maxil|um 

mobilities  exhibit  a  slight  local  minimum  in  the  vicinity  of  tU1,Si. 

'-r  n-,ype  Si,  the  effect  of  thermal  stress  on  the  maximum 

mobility  ranges  from  a=)30%  reduction  in  mobility  for  (ooi)Si  up  to 

a  ~  2  enhancement  in  mobility  near  fnnKi  t- 

'  near  UlOJSi.  Tne  ratio  of  (110)  Si 

l)Si  mobilities  ^max/UQ  is  ^ 1 . 7 -  1 . 9  which  represents  a  sub- 

stantial  and  significant  increase  in  theoretical  mobilities. 

Hcasuicmcnts  '  in  these  laboratories  tend  to  support  this  general 

characteristic  in  that  n-type  (221 3  Si  films  have  yielded  higher 

mobilities  than  ( 1 1 1 ) Si / ( 1 1 20) A t 2„,  films, which  in  turn  have 

yielded  higher  mobilities  than  (001)Si  films. 

'or  p - 1 ype  Si.  the  effect  of  thermal  stress  is  always  to 

produce  an  enhancement  in  both  maximum  and  minimum  mobilities 
e  e  * 

Umax/Mo  and  Umin/Mo*  The  maximL’m  mobility  /  \x  is  -1.07  for 

(OOl)Si  and  increases,  as  one  proceeds  across  the  defined  zone,  to 

=  2.-1  for  CHOJSi.  These  results  indicate  that  the  thermal  stress 

effect  on  mobility  has  a  strong  orientation  dependence  and 


that  substantial  larger  mobilities  may  be  achieved  by  choice  of 
film/substrate  orientation. 

Collectively,  the  theoretical  thermal  stress  mobility  data 
for  n-  and  p-type  Si/Af^  suggest  that  i,  would  be  worthwhile  to 
experimentally  explore  Mode  I  epitaxial  film  growths  in  the  range 

of  -0-25  deg  from  (001)Si  and  Modes  II,  II,  and  IV  epitaxial  film 
growths  in  the  range  +  2S  deg  from  (lll)Si. 

Summarizing,  the  effect  of  sub s t rat e - induced  thermal  expansion 
ress  on  mobility  in  Si/AljOj  films  has  been  theoretically  determined 
Within  the  framework  of  a  t h in- f i 1 m- 1 h i ck - subs t rat e  piezoresistance 
model  which  incorporates  the  anisotropy  in  substrate  thermal 
expansion  and  in  Si  film  stresses. 

It  is  well  known  that  Si/Al^  film  growth  quality 
dependsupon  a  number  of  factors  including  intrinsic  stress,  dislo¬ 
cation  generation  and  impurity  or  heavy  natal  complexes,  in  addition 
to  the  thermal  stress  effect  treated  in  this  paper.  All  of  these 
phenomena  could  produce  effects  on  mobility  which  are  dependent 
upon  crystallographic  orientation.  Although  a  formidable  task, 
additional  theoretical  and  experimental  work  directed  toward  under- 
standing  these  phenomena  should  be  attempted. 

The  present  study  of  thermal  stress  effects  represents  the 
first  systematic  investigation  of  mobility  as  a  function  of 
orientation.  We  believe  the  model  is  a  useful  theoretical  tool, 
and,  in  conjunction  with  experimental  investigations  and  additional 
selected  orientation  growth  studies,  can  lead  to  a  better  under¬ 
standing  and  to  means  of  optimizing  Si/AtjOj  films  for  some 
appl i cat  ions . 
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figure  captions 


Fig  . 


Fig.  2 


Fig.  3 
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equivalent  pole  diagrams  characterizing  an  off- 
orientation  *<111)S1  film  growth.  The  parallel  relatlons 

associated  with  these  three  pole  diagrams  are: 

A>  !1iOISl//[2201]A2„Oj  and  ( 1 11 ) s 1/, ( llio) A tj o3 , 

B)  [I01|Sl//t2201|A2203  and  ( 11 1 ,  s  1/ / (1 120)  AH^ ; 

O  U10]S1//[2021JA2203  and  (111) SI// (1 210) AIjOj . 

Although  all  three  diagrams  are  equivalent,  only  that 
°f  Fig.  1A  Is  consistent  mathematically  with  the 
Parallel  relations  assumed  In  Section  U  and  has  been 
employed  in  the  calculations. 


ill 


ratio 


f0r  n"type  Si  is  plotted 


versus  orientation  in  the 

tne  (xxl)Si  zone.  The  shaded 

regiois  indicate  the  r  n  1  1  o  -i 

the  collective  range  in  mobility  for 

the  four  modes  of  Si/Ai  n 

/A*2°3  epitaxy.  (Data  for  the 

individual  modes  Is  given  In  the  tables.) 

The  mobility  ratio  ( H * / U 0 )  based  upon  the  exponential 
non-linear  correction  factor  for  p-type  SI  Is  plotted 

versus  orientation  In  the  (xxl)sl  rone.  The  shaded 

regions  indicate  the 

the  collective  range  in  mobility  for 

the  four  modes  of  Si/A£  n 

/A*2°3  epitaxy.  (Data  for  the 

individual  modes  is  given  in  the  tables.) 
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Table  I.  General  Silicon-Sapphire  Orientation  Relationships 


Mode 


Orientation  Planes  Orientation  Direction 
Silicon//Sapphire  Si licon//Sapphire 


I 

II 

II 

III 

IV 


(001)// (0112) 

(111)  7(1124) 

(no)//  ( 1 1 2 o ) 

i  -  2 1 )  /  /  ( 1 1 2  2 ) 

(111)  //(il2o) 

(112) //(i!o5) 

(111)// (1014) 


[ 100] // [ 21 10] 

[ llo] // [ lloo] 

[ llo] // [ llco] 

[ 110] // [2201] 

[ 110]  // [ i2io] 

[111) // [2021  ] 
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TABLE  II.  Parallel  Directions  and  Transformation 


£  I  i  [ 1 1° ] ;  X'\  | [ 1 lOO ] 4  and  |  |  [ /T  ,  -l,0]c  For  (SF) ,  use  ( 1 1 1 ) Si |  |  ( 1 1 24) Ai 

^ |  |  [114J;  Y'||[1123]4  and  |  | [ 1 ,  S% ,  -2(c/a)Jc  and  see  text  for  explanation 

f  J.(221 )  ;  Z'-t(H22)  and  |  |  [  1 ,  /1 ,  2(a/c)l 


TV(OJII) | JTS(III)  III  3P°w 


Table  III.  Comparison  of  Linear  and  Non-Linear  Mobilities 


(m/l0) 

Linear 


0.70 
0.80 
0.90 
1 . 00 
1  .  10 
1 .  20 
1  .  30 
1  .  40 
1  .  50 
1.75 
2. 00 
2.50 
3.00 
4.00 
5.00 
6.00 
8.00 


00 


(y/M0) 

Second  Order 


(ii/n0) 

Exponential  Order 


0.658 
0 . 780 
0.895 
1 . 000 
1 . 095 
1.180 
1 .  256 
1.324 
1 .385 
1.508 
1 .600 
1  .  724 
1  .  800 
1  .882 
1.923 
1 .946 
1.969 

2.000 


0.651 
0.779 
0 . 895 
1 . 000 
1 . 095 
1.181 
1  .  259 
1.331 
1  .  396 
1 .535 

1  .649 
1  .  822 
1.948 
2.117 
2 .225 
2.301 

2  .  399 

2.718 


i 


TABLE  IV 

(001)Si//(0ll2)Al20 


Case  1 

N-type  Si 

Ap{/po 

=  +.441924  -.044488 

cos  2  0 

E2^PoJ1 

=  +.044488  sin29 

Ap3/po 

=  9  S’  4  3  2  9 

Case  1 

P-type  Si 

ap;/p0 

=  -.060424  +.002186 

cos2u 

E2/poJI 

=  -.002186  sin2  9 

iPj/Po  ’ 

=  +.010440 

Case  2 

N-type  Si 

S/Oo  : 

=  +.437032  -.066732 

cos29 

E2/poJ1 

=  +.066732  sin29 

Ap3/po  - 

-  -.983323 

Case  2 

P-type  Si 

Apl/po  =  "-059755  +.00328  cos26 
E2/,poJl  =  --00328  sin2 0 
Ap3/pc  =  + • 010325 


TABU!  V 


( 2  2 1 ) S i |  |  (1122)A£70? 

C <i s_e  1  n -  type  Si 

^ 3 / P 0 J i  =  -.15  3  2308  sinO 

^2Apo‘*l  =  +.2122694  sin20 

Apj/c^  =  +.06  1  8867  -.2  122694  cn^O 

u—  •  *•  1770V  “min  “0. 784831, o;  p'aj  =  1.1623„  ; 


max 


mi  n 


iA  =0. 98092  p  ;  A  =  39 .98' 


1 0 . 9  6  1  2  5  p 


=  41 


Case  2  n-typc  Si 

1-  ^  /  *  0 -1  i  -.1640888  sinO 

E',Apo^l  =  +.2314  3  87  sin26 

Ari/Po  =  +  -0(,588(j4  -.23  1  2  447  cos  2  0 

C  .  =  1 •  1 9  8 1  u  ;  p  .  =0.7709 3 p  ;  pC 

max  o  min  Mo ’  Mmax 


1.1798U  ;  pe. 

o’  Mmin 


UA  =0,<-)8452  o  >  A  =  45-39?«  I  PA  =  0 . 9  6  1  3  7  p  o ;  Ae 

Case  1  p-typc  Si 

1  3  /  DoJ  1  =  '•  1406087  sinO 

E2/po'J  i  =  "  •  2  4  2  5  34  2  sin20 

Ap!/po  r  * .5286124  -.2425342  cos20 

■  4-3696V  -.in  ■  •  2.H.22V  . 

UA  -  :.884-juo  ;  A  •  102. 9",  ;  v*  =  1 . 7407uo  ;  A1’  .  4: 


-ise  2  p-typc  Si 
E3/PoJ1  =  ~ • 12285  sinO 
E2Apo^l  =  + • 2281746  sin26 
Apl/Oo  =  '• 8 210478  -.228174  co?2( 


ma 


v  -  3.9876po;  p^  =  1.4142po;  %ax  =  2.1153^, 


aA  =  2.7009pq  ;  A  =  95.28% 


JA  =  l-7278po 


e 

P  • 
min 


=  44 


0 . 7602  1  p 
.83% 


0  .  74  2 9 4  p 
4  5.4  4% 


1  .  3  3  1  2  p 

o 

5  7% 


1  .  3403po 
.  85% 
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TABLE  VI 


(lll)Si/A420 


Mode  II 

Mode  III 

Mode  IV 

N1 

0.0 

♦0.618993 

0 . 0 

N  2 

-0 . 80673 

-0.785396 

-0.618993 

Case  1 

(T1+T2)/2 

-1.1671xl010 

-  1 .  19  87xl010 

-  1 . 14  2S  xlO10 

CTj-T2)/2 

♦0. 34343xl09 

♦0. 12333xl09 

♦0. 20220xl09 

T  , 

0 

0.0 

♦0.51312X109 

0 . 0 

Case  2 

CVT2)/2 

-  1  .  1513xl010 

-  1 . 1987xl010 

-1. 1150X1010 

(TrT2)/2 

>0 . 5151Sxl09 

♦0 . 18499xl09 

♦0. 30330X109 

T6 

0.0 

♦0. 76968xl09 

0.0 

n-type  Si  Case  1 

xr 

3 

P 

X 
• — 
xr 
o 

1.013 

1 .025 

1 . 005 

lmin^o 

0.972 

0.961 

0.9804 

A  (%) 

4.18 

6.42 

2.46 

n-type  Si  Case  2 

u  /u 

max  o 

1.024 

1.042 

1.011 

u  .  /u 

min  o 

0.962 

0.946 

0.974 

A  (%) 

6.27 

9 .63 

3.69 

p-type  Si  Case  1 

^max^o’  ^max^o 

2.604 

1.851 

2.851  1.914 

1  2.442  1.805 

^min^o'  ^min^o 

2.227 

1.735 

2.219  1.732 

|  2.234  1.737 

A  (%)  ;  Ae  (%) 

15.60 

6.50 

24.91  9.98 

|  8.92  3.81 

p- type  Si  Case  2 

1 

u  /u  ;  Ue  /y 

mi x  o’  max' Mo 

2.662 

1.867 

3.069  1.962 

1  2.416  1.797 

um,_/yn;  y®.  /y„ 

min  o  min  o 

2.  114 

1.694 

2.103  1.689 
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ELECTRON  MICROSCOPE  IN  SITU 
NUCLEATION  AND  GROWTH  STUDIES 


Introduction 


•  ...  A  seardi  of  the  technical  literature  involving  the  electron  microscopic 

in  situ  observation  of  the  nucleation  and  growth  of  thin  films  has  been  performed 
The  resultant  bibliography,  it  is  hoped,  includes  all  previous  research 
interest  to  someone  undertaking  such  a  task.  Pertinent  literature  in  all  maior 

electron  microscopy  conferences  listed  in  section  G  has  been  included,  as  wel? 

as  n^nerous  articles  from  the  general  scientific  literature.  Insofar  as  is 
possible  for  someone  in  the  physical  sciences  on  a  limited  time  budget,  the 
lological  literature  was  searched  for  information  on  hydration  or  environ- 

orni-mJro5ren?Snf  ^  artjcleS  encountered  concerning  auxiliary  pumping  systems 
rpn-^q  0VfTe  1  ?!  the  c°"ta™natlon  rate  were  included,  primarily  in  section 
article  not^in^dV"^  ^  aPPreciate  ' .iformation  concerning  any  pertinent 

To  date,  no  one  has  yet  examined  i_n  situ  the  nucleation  and  growth  of  a  thin 

Tilin  from  a  gaseous,  molecular  phase  onto  an  electron  transparent  substrate 
aP’/he  excenent  works  of  Pashley  et  al  ,  and  Poppa  (references  C1-C8 

deposition  specifl^d1]^  t0  Physical  vapor  deposition,  not  chemical  vapor 

eposition.  Tie  experiments  listed  in  section  E  concern  either  (a)  oxidation  of 

a  t Mn  film,  (b)  growth  of  solid  whiskers  from  pyrolysis  of  a  gaseous  phase 

(cl  decomposi  tion  of  a  solid  phase  into  a  different  solid  phase  plus  a  gaseous 
phase.  None  concern  the  growth  of  a  thin  solid  film  from  the  pyrolysis  oU 
g  eous  phase,  an  in  situ  experiment  which  the  author  is  undertaking.  The 
rt  erences  citea  have  been  valuaLle  chiefly  for  their  collective  wealth  of  ideas 
and  experimental  techniques. 
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D7.  G.  A.  Bassett,  Condensation  and  Evaporation  of  Solids,  Ed.  E.  Rutner, 

P.  Goldfinger,  and  J.  P.  Hirth  (Gordon  and  Breach,  New  York,  1962), 
p  599.  "Continuous  Electron  Microscope  Study  of  Vacuum  Evaporated 
Metal  Films,: 

D8.  D.  W.  Pashley  and  M.  J.  Stowell,  p  GG  1  in  reference  G3  (5th  Int'l, 
Philadelphia,  1962).  "The  Growth  of  Evaporated  Silver  Layers  Inside 
an  Electron  Microscope." 

D9.  G.  A.  Bassett,  p  270  in  reference  G7  (2nd  European,  Delft,  1960). 

"Continuous  Observation  of  the  Growth  of  Vacuum  Evaporated  Metal  Films." 

DIO.  H.  Poppa,  Zeit.  Naturforsch.  19a,  835  (1964).  "In-Situ  Studies  of 
Epitaxial  Thin  Film  Growth."  " 

Dll.  H.  Poppa,  J.  Appl .  Phys.  38,  3883  (1967).  (See  Also  National  Aeronautics 
and  Space  Agency  Technical  Note  D-4506,  1968).  "Heterogeneous  Nucleation 
of  Bi  and  Ag  on  Amorphous  Substrates  (In-Situ  Electron  Microscopy)." 

D12.  A.  E.  Curzon  and  K.  Kimoto,  p  507  in  reference  G2  (6th  Int'l,  Kyoto,  1966) 
"The  Condensation  of  Liquid  Lead  on  Molybdenite  in  an  Electron  Microscope. 

D13.  B.  Coopersmith,  A.  E.  Curzon,  K.  Kimoto,  and  N.  D.  Lisgarten,  Basic 
Problems  in  Thin  Film  Physics,  Ed.  R.  Niedermayer  and  H.  Mayer 
(Vandenhoeck  and  Ruprecht,  Gottengen,  1966),  p  83.  "The  Condensation 
of  Lead  Inside  an  Electron  Microscope." 
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D14 , 

D15. 

D16. 

D17. 

D18 . 


D19 , 


D20. 

021  . 

D22 . 

D23 . 

024. 

D25. 

D26 . 

027. 


pda  S;.Sennet^«  and  G-  Scott,  Canadian  J.  Research 

— ’  r  U-SO).  Continuous  Observations  with  the  electron  Microscope 
on  the  Formation  of  Evaporated  Films  of  Silver,  Gold,  and  Tin." 

A.  barna,  P.  B.  Barna,  and  J.  F.  Pocza,  p  445  in  reference  G1  (7th  lnt'l, 
Grenoble,  1970).  Simultaneous  Investigation  of  Structure  and  Electrical 
Properties  of  Vacuum  Deposited  Thin  Films  by  In-Situ  Electron  Microscopy." 

T.  E.  Hutchinson,  p  370  in  reference  G12  (26th  EMSA ,  New  -leans,  1908) 

In  svtu  E  ectron  Microscope  Investigations  of  the  Growtn  jnd  Structure 
ot  inin  r i Ims . 

H  Poppa  p  378  in  reference  G12  (26th  EMSA,  New  Orleans,  1968).  "Thin 
Film  Nucleation  and  Growth  Under  Controlled  Conditions." 

fh  Ff  PrCZ5’  A'  Barna*  and  P-  Ba™a>  Basic  Problems  in  Thin  Film  Physics, 
lQfif ‘n  cqenn^er\an?-H‘  Hayer  (VancJenlioeck  aTTT %7ecBt, “Gottingen, 

I960,  p  153.  Nucleation  and  Growth  Processes  in  Vacuum  Deposited 
Germanium  Films . 

T.  E.  Hutchinson  and  0.  Sherman,  U.  Minnesota  Annual  Technical  Prooress 
Report  III-  COO-179C,  3  Jan.  1971;  IV  -  C00-1790  -  4  Jan.  1972, 

In  Situ  E. ectron  Microscope  Investigation  of  the  Nucleation  and  Growth 
of  Sputtered  Thin  Films." 

D.  M  Sherman  and  T .  E.  Hutchinson,  p  516  in  reference  G10  (28th  EMSA, 
Houston,  1970).  Nucleation  and  Early  Growth  of  Sputtered  Thin  Films." 

?q75;  M°°rhe?d*  H:  P°PPa*  P  516  in  reference  G15  (30th  EMSA,  Los  Angeles, 
1972).  In  Situ  Deposition  of  Gold  on  MgO."  y 

H.  Sato  and  S  Shinozaki,  J.  Vac.  Sci .  and  Technol .  8,  159  (1971  ). 
Morphology  of  Nuclei  and  Epitaxial  Behavior  of  Au  and  Ag  on  MgO." 

J?  nncQ\S’  Shinozaki  and  L.  J.  Cicotte,  J.  Vac.  Sci.  and  Technol.  6, 

62  (1969).  Direct  Observation  of  Epitaxy  on  MgO." 

GnnH°n^0*Aund  H-  Sat0’  p  505  in  reference  G2  (6th  lnt'l, 
Kyoto,  1966).  Direct  Observation  of  the  Epitaxy  of  Ag  on  MgO  Substrates." 

S.  Shinozaki  and  H.  Sato,  J.  Appl.  Phys.  43,  701  (1972)  "Cleaving  of  MgO 
inside  an  Electron  Microscope." 

r0bay^h1\n;Jak^fna9il  and  G-  HonJ°*  P  439  in  reference  G1 
(7th  Int  1,  Grenoble,  1970).  In  Situ  Observation  of  Epitaxial  Growth  of 
Metals  and  Non-metals. 

j’  ['  Poc“’  A.  Barna,  and  P.  B.  Barna,  J.  Vac.  Sci.  Technol.  6,  472 
vl969).  Formation  Processes  of  Vacuum  Deposited  Indium  Films~and 

Thermodynamical  Properties  of  Submicroscopi c  Particles  Observed  by  In  Situ 
Electron  Microscopy. 
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t.  (1-19)  Related  Topics:  Gas  Inlet  Devices 


El.  D.  W.  Pashley  and  A.  E.  B.  Presland,  Phil.  Hag.  7,  1407  (1962).  "The 

Movement  of  Dislocations  During  the  Observation  of  Metal  Films  Inside  an 
Electron  Microscope." 

E2.  H.  G.  Heide,  J.  Cell  Biology  1_3,  147  (1962).  "Electron  Microscopic 
Observation  of  Specimens  Ur.der  Controlled  Gas  Pressure." 

E3.  H.  G.  Heide,  p  87  in  reference  G4  (4th  Int'l,  Berlin,  1958).  "Die 
Objektverschmutzung  und  ihre  Verhutung." 

E4 .  I.  G.  Stoyanova  and  G.  A.  Mi khai lovski ,  Biofizika  4,  483  (1959). 

E5.  J.  L.  Williams,  K.  Heathcote,  and  E.  J.  Greer,  p  100  in  reference  Gil 

(27th  EMSA,  St.  Paul,  1969).  "A  New  1000  kV  Electron  Microscope." 

E6.  I.  G.  Stojanowa,  p  82  in  reference  G4  (4th  Int'l,  Berlin,  1958).  "Ein 
Kamner  fur  die  Untersuchung  von  Objekten  mit  Gasumgebung  (A  Chamber  for 
the  Examination  of  Objects  Surrounded  by  Gas)." 

E7 .  H.  G.  Heide,  Naturwissenschaf ten  47,  313  (1960).  "Electron  Microscopy 
of  Objects  Under  Atmospheric  Pressure  or  Pressures  Which  Prevent  Drying 
of  Objects ." 

E8.  V.  R.  Matricardi,  G.  G.  Hausner,  and  D.  F.  Parsons,  p  542  in  reference 
G10  (28th  EMSA,  Houston,  1970).  "Hydration  Chamber  for  Jeolco  200kv 
Microscope." 

E9 .  R.  C.  Moretz,  G.  D.  Hausner,  and  D.  F.  Parsons,  p  544  in  reference  G10 
(28th  EMSA,  Houston,  1970).  "Studies  on  Water  in  the  Hydration  Chamber 
of  a  Modified  Electron  Microscope." 

E10.  R.  C.  Moretz,  G.  G.  Hausner,  and  D.  F.  Parsons,  p  44  in  reference  G14 
(29th  EMSA,  Boston,  1971).  "New  Design  for  a  Differentially  Pumped 
Hydration  Chamber  for  a  Siemens  1A." 

Ell.  D.  F.  Parsons  and  R.  C.  Moretz,  p  497  in  reference  G1  (7th  Int'l, 

Grenoble,  1970).  "Microscopy  and  Diffraction  of  Water  in  the  Electron 
Microscope." 

El 2 .  D.  L.  All inson,  p  169  in  reference  G1  vol  I  (7th  Int'l,  Grenoble,  1970). 
"Environmental  Cell  for  Use  in  a  High  Voltage  Electron  Microscope." 

E13.  A.  Fukami,  K.  Adachi,  and  M.  Katoh,  p  263  in  reference  G2  (6th  Int'l, 

Kyoto,  1966).  "On  a  Study  of  New  Micro  Plastic  Grid  and  its  Applications. 

E14.  A.  Fukami,  T.  Etoh,  N.  Ishihara,  M.  Katoh,  and  K.  Fujiwara  p  171  in 
reference  G1  (6th  Int'l,  Kyoto,  1966).  "Pressurized  Specimen  Chamber 
for  Electron  Microscope." 


Preceding  page  blank 
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G.Dup°uy,  F.  Perrier,  and  L.  Durrieu,  Comptes  Rendus  251_,  2836  (1960). 

L  Observation  de  la  matifere  vivante  au  moven  d'un  microscope  £lectronique 
fonctionnant  sans  tres  haute  tension."  (The  Observation  of  Living 
Material  by  Means  of  an  Electron  Microscope  Operating  at  Very  Hiqh 
Voltage. )' 


G.  Dupouy,  F.  Perrier,  and  L.  Durrieu,  Comptes  Rendus  254,  3786  (1962). 

L  Observation  des  objets  en  milieu  gazeur.  Application  h  l'gtude  de  la 
contamination  dans  le  microscope  glectronique."  (The  Observation  of 
Objects  in  the  Midst  of  GaAs.  Application  to  the  Study  of  Contamination 
in  the  Electron  Microscope." 

P.  R.  Ward  and  R.  F.  Mitchell,  p  44  in  reference  G16  (25th  EMAG,  1971). 
An  Electron  Microscope  Environmental  Specimen  Chamber." 

and  M-  ^at0h’  p  612  in  reference  G15  (30th  EMSA,  Los  Angeles, 
1972).  Construction  and  Application  of  Environmental  Cell." 

V*  n,a^r1'card''  *  Subjeck,  and  D.  F.  Parsons,  p  180  in  reference  G15 
(30th  EMSA,  Los  Angeles,  1972).  "Hydration  Stages  -  Achievements  and 
Obstacles." 


A.  Fukami  and  K.  Adachi,  J.  Electron  Microscopy  J4_,  112  (1965).  "A  New 
Method  of  Preparation  of  a  Self  Perforated  Micro  Plastic  Grid  and  Its 
Application." 

A.  Fukami,  K.  Adachi,  and  M.  Katoh,  J.  Electron  Microscopy  21_,  99  (1972). 
Micro  Grid  Techniques  and  Their  Contribution  to  Specimen  Preparation 
Techniques  for  High  Precisions  Work." 
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L_(  20-39) — jjejated  Topics: _ Ion  Pimping  and  Residual  Gas  Analysis 


1 20 , 

E21 

E22 , 

E23 . 

E24 . 

E25 . 

E26 . 

E27 . 

E28. 

E29 . 

E30 . 

E31 . 
E32 . 

E33 . 


k.E.  Hartman,  H.  Akahori,  C.  Garrett,  R.  S.  Hartman,  and  P  L  R„mn< 

St-  19H>-  S*  ^Upe"^' 

II.  The  Design  of  a  Gas  Control  Chamber  for  Quantitative  Observations." 

J.  A  Maliakal,  p  78  in  reference  Gil  (27th  EMSA,  St.  Paul,  1969) 
Application  of  Ion  Pumps  for  Electron  Microscopes." 

5*  h?2£?n  aS  R‘  Sl  Hartman«  P  159  in  reference  G2  (6th  Int'l 

Microscope/  '  ReSldua1  Gas  Anal*sis  in  an  Ion  Pumped  Electron 

(6th '  Int^l *  LL  K?ofiSMr’  and  K;  fiauri"’  P  161  in  reference  G2 
Si emens ‘ei ectron°Mi croscope .  *S  1  ^  ^  A"alySiS  in  *"  Auxilia^  PumPed 

R.  E.  Hartman  and  R.  S.  Hartman,  J.  Appl .  Phys  37  3936  (1966) 

Of  two  Electron  Scopes  to 

/iQ™iMar"7S’  C'  J'  Ca1bick‘  and  T.  T.  Sher.g,  t>v.  Sci .  Inst  41  1319 

VibralionlessXeTtfon!"  a"  E'eCtr'0"  MiCr°Scope  f°r  tlean.  Qu1?{. 

C'  J'  Calbick>  and  T.  T.  Sheng,  p  163  in  reference  G1 
for  Clean, 'Qu)e t°Vi bra tionless  Operation?"5  a"  ^'-^on  Microscope 

5e1  lil  .  S1965-*  • 

Boston'iaf9711  !?arbna[J  ’  p  74  in  reference  G14  (29th  EMSA, 

iv.  A'factoMn  RedifSLagT”'0"5  the  EleCtr0"  Hicr°sc°Pe: 

(?6fh  Rn  ?•  Hartnan  and  P-  L-  Ramos»  p  292  in  reference  G12 

(26th  EMSA,  New  Orleans,  1968).  "Residual  Gas  Reactions  in  the  Electron 
Microscope:  I.  Qualitative  Observations  on  the  Water  Gas  Reaction." 

5;  uinErSSk1,  P  ?14Jn  r,eference  G 1 2  (26th  EMSA,  New  Orleans,  1968) 

A  High  Vacuum  Electron  Microscope."  ' 

E.  A  Buvinger,  p  348  in  reference  G12  (26th  EMSA,  New  Orleans  1968) 
"Background  Gas  Reduction  in  the  Electron  Microscope."  ' 

Et.V3:1’-1-/-  Kassner’  and  J-  K>  Maurin,  Phil.  Mag.  21  ,  453  (1970) 

The  Contamination  of  Surfaces  during  High  Energy  Electron  Bombardment.’" 
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E34 .  J.  R.  Parsons  and  C.  W.  Hoelke,  Can.  J.  Phys.  44,  2685  (1966). 
"Electron  Microscopic  Observation  of  Atomic  Planes  and  Atomic 
Positions  in  Germanium." 

E35 .  G.  M.  Greene,  J.  W.  Sprys,  and  C.  D.  Melvin,  p  620  in  reference  G15 
(30th  EMSA,  Los  Angeles,  1972).  "Reduction  of  Contamination  Rates 
in  Electron  Microscopes." 

E36.  P.  R.  Swann  and  N.  J.  Tighe,  in  reference  G17  (5th  European, 

Manchester,  1972).  "Performance  of  a  Differentially  Pumped  Environ¬ 
mental  Cell  in  the  AEI  EM7 . "  See  also  ref.  A27. 
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L40.  A.  W.  Agar  and  J.  H.  Lucas,  p  E2  in  reference  G3  (5th  Int'l,  Philadelphia, 
1962).  "Use  of  a  New  Heating  Stage  for  the  Electron  Microscope." 

E41.  ;.  0.  McPartland,  p  E3  in  reference  G3  (5th  Int'l,  Philadelphia,  1962). 

"A  High  Temperature  Stage  for  Transmission  Electron  Microscopy." 

E42.  N.  Takaliashi,  T.  Takeyama,  K.  Ito,  T.  Ito,  K.  Mihama,  and  M.  Watanabe, 
Japan.  J.  Electronmicroscopy  4_,  16  (156).  "High  Temperature  Furnace 
for  the  Electron  Microscope  (JEM  5)." 

E43.  N.  Sasaki  and  R.  Uyeda,  Rev.  Sci .  Instr.  23,  136  (1952).  "A  Specimen 
Treating  Adaptor  for  the  Electron  Microscope." 

E44 .  I.  Okazaki,  M.  Watanabe,  and  K.  Mihama,  p  93  in  reference  G4  (4th  Int'l, 
Berlin,  1958).  "Improvement  of  the  Specimen  Heating  Device  for  the 
Electron  Microscope." 

E45.  K.  Itoh,  T.  Itoh,  and  M.  Watanabe,  p  658  in  reference  G9  (3rd  Int'l 
London,  1954). 

E46.  M.  J.  Whelan,  p  96  in  reference  G4  (4th  Int'l,  Berlin,  1958).  "A  High 
Temperature  Stage  for  the  Elmiskop  I." 

E47.  W.  Loebe,  0.  Schott  and  F.  Wilko,  p  131  in  reference  G7  (2nd  European, 
Delft,  1960).  "Objektheizungseinvichtung  zum  Elmiskop  I  (Specimen 
Heating  Stage  for  the  Siemens)." 

E48.  R.  M.  Fisher,  P.  R.  Swann,  and  J.  Nutting,  p  131  in  reference  G7  (2nd 
Liropean,  Delft,  1960).  "A  New  Objective  Pole-Piece  and  Specimen 
Heating  Stage  for  the  Elmiskop." 

E41.  U.  Valdre,  J.  Sci.  Instrum.  42,  853  (1965).  "A  Double  Tilting  Heating 
Stage  for  an  Electron  Microscope." 

L50.  U.  Valdre,  p  165  in  reference  G2  (6th  Int'l,  Kyoto,  1966).  "A  Universal 
Specimen  Stage  and  Combined  Cartridges  for  an  Electron  Microscope." 

E51.  U.  Valdre,  II  Nuovo  Cimento  X53B,  157  (1968).  "Combined  Cartridges 
and  Versatile  Specimen  Stage  for  Electron  Microscopy. 

E52.  D.  Jones,  p  662  in  reference  G9  (3rd  Int'l,  London,  1954).  "A  Specimen 
Furnace  for  the  Reflection  Electron  Microscope." 

E53.  I.  W.  Sprys  and  P.  C.  J.  Gallagher,  p  312  in  reference  G12  (26th  EMSA 
New  Orleans,  1968).  "A  New  Hot  Stage  for  the  Philips  EM200  and  its 
Calibration." 

E54.  V.  Valdre,  p  317  in  reference  G17  (5th  European,  Manchester,  1972). 

General  Considerations  on  Specimen  Stages." 
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E.  (60-79)  Related  Topics:  Cinematography 

E60.  V.  A.  Phillips,  Phil.  Mag.  !5,  571  (1960).  "Direct  Observation  of 
Defects  in  Evaporated  Silver." 

E61.  R.  W.  Horne  and  R.  H.  0 1 tew ill,  p  140  in  reference  G4  (4th  Int'l, 

Berlin,  1958).  "Cinematographic  Studies  on  the  Interaction  of 
Electrons  with  Microcrystals  of  Silver  Iodide." 

E62.  K.  F.  Hale,  M.  H.  Brown,  and  D.  L.  Allinson,  p  297  in  reference  G1 
(7th  Int'l,  Grenoble,  1970).  "Resolution  and  Contrast  Experiments 
with  Various  Arrangements  of  Environmental  Cell  Windows  at  High 
Vol tages . " 

E63.  E.  Sugata,  Y.  Nishitani,  S.  Kaneda,  M.  Tateishi,  and  H.  Yoloya,  p  452 
in  reference  G9  (3rd  Int'l,  London,  1954).  "Fundamental  Researches 
for  Observing  Specimens  in  Gas  Layers." 

E64.  F.  D.  Lugton  and  C.  E.  Warble,  Rev.  Sci.  Instr.  4J_,  1793  (1970). 

"High  Resolution  Cinematography  at  High  Temperature  in  the  Electron 
Microscope . " 

E65.  M.  J.  Whelan,  P.  B.  Hirsch,  R.  W.  Horne,  and  W.  Bollmann,  Proc.  Royal 
Soc.  (London)  A240,  524  (1957).  "Dislocations  and  Stacking  Faults  in 
Stainless  Steel ." 

E66.  K.  Anderson,  P.  B.  Kenway,  p  244  in  reference  G12  (26th  EMSA,  New  Orleans, 
1968).  "External  Photography  of  the  Microscope  Image." 

E67.  J.  W.  S.  Hearle,  U.  J.  Clarke,  B.  Lomas,  D.  A.  Reeves,  and  J.  T.  Sparrow, 
p  210  in  reference  G16  (25th  EMAG,  1971).  "A  Simple  Method  of  Recording 
Dynamic  Events  in  the  SEM  Using  a  16  mm  Camera." 

M.  Heerschap  and  R.  DeCat,  p  170  in  reference  G17  (5th  European, 
Manchester,  1972).  "Direct  Filming  of  Transient  Phenomena  with  a 
Closed  TV  Circuit  as  Viewfinder." 
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F.  Rplated  Topics:  Cooling  Stages  (Partial  Listing) 

FI.  J.  A.  Venables  and  D.  J.  Ball,  Proc.  Royal  Society,  London  A322,  331 
(1971).  "Nucleation  and  Growth  of  Rare  Gas  Crystals. 

F2.  J.  A.  Venables,  Phil.  Mag.  21_,  147  (1970).  “Electron  Micrscopy  and 
Diffraction  of  Solid  a^." 

F3.  J.  A.  Venables,  D.  J.  Ball,  and  G.  J .  Thomas,  J.  Phys.  (Sci.  Instr.) 

El,  121  (1968).  "An  Electron  Microscope  Liquid  Helium  Stage  tor  use 

wTth  Accessories." 

F4.  D.  J.  Ball  and  J.  A.  Venables,  p  239  in  reference  G5  (4th  European,  ^ 
Kane,  1968).  "The  Use  of  a  Pressure  Cell  for  Studying  Condensed  Gasis. 

F5_  j.  a.  Venables  and  D.  J.  Ball,  J.  Crystal  Growth,  3,  4,  180  (1969). 
"Microstructure  of  Condensed  Gases." 

F6 .  J.  A.  Venables,  p  175  in  reference  G2  (6th  Int'l,  Kyoto,  1966). 

"A  Goniometer  Stage  for  Use  at  Liquid  Helium  Temperatures. 

F7  j  A  Venables  and  D.  J .  Ball,  p  333  in  reference  G2  (6th  Int'l,  Kyoto, 
1966).  "Defects  in  Rare  Gas  Crystals. 

F8.  U.  Valdre  and  M.  J.  Goringe,  J .  Phys. (Sci.  Instr.)  336  (1970). 

"An  Improved  Liquid  Helium  Stage  for  an  Electron  Microscope. 

F9.  U.  Valdre  and  M.  J.  Goringe,  0.  Sci.  Instr.  42,  268  (1965).  A  Liquid 
Helium  Cooled  Goniometer  Stage  for  an  Electron  Microscope. 

F10  U  Valdre,  p  61  in  reference  G6  (3rd  European,  Prague,  1964)  .  (See 

also  references  E50  and  E51).  "A  Double  Tilting  Liquid  Helium  Cooled 
Object  Stage  for  the  Siemens  Electron  Microscope." 

Fll.  P.  R.  Swann,  in  reference  Gil  (27th  EMSA,  St.  Paul,  1969).  A  Tilting 
Cold  Stage  for  the  AEI  EM802." 

F12.  R.  M.  J.  Cotterill,  p  63  in  reference  G6  (3rd  European,  Prague,  1964). 
"A  Liquid  Helium  Stage  for  the  Siemens  Elmiskop  I." 

F13  G.  R.  Piercy,  R.  W.  Gilbert  and  L.  M.  Howe,  J.  Sci.  Instr.  40,  487 
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APPENDIX  6 


ELECTRON  SCATTERING  BY  GASEOUS  ATOMS 


The  presence  of  a  gas  ph.’.se  near  the  specimen  in  an  electron  microscope 
causes  unwanted  electron  scattering  of  the  primary  beam,  resulting  in  a 
loss  of  resolution.  More  importantly,  there  is  an  intensity  loss  in  the 
transmitted  beam  and,  also,  a  background-intensity  increase  and  partial 
destruction  of  contrast  information  because  of  multiply-scattered  electrons, 
i.e.,  electrons  originally  scattered  out  of  the  primary  beam  but  rescattered 
into  the  transmitted  beam. 

The  beam  intensity  loss  due  to  gas-phase  scattering  has  been  calculated 
and  used  in  determining  the  upper  limits  of  the  gas  pressure  and  gas 
layer  thickness  which  can  be  tolerated  in  an  ^n  situ  CVD  experiment. 

These  two  related  variables  determine  the  dimensions  of  that  portion  of 
the  specimen  chamber  which  contains  the  substrate  and  the  surrounding 
reactant  gases,  and  through  which  the  electron  beam  must  pass. 

The  electron  scattering  for  a  single  molecule  is  a  function  of  the  electron 
scattering  cross-section  a  for  the  molecule,  which  varies  strongly  with  the 
energy  of  the  incident  electrons,  and  the  number  and  distribution  of  electrons 
in  the  molecule.  The  total  cross  section  for  a  volume  of  gas  containing 
N  atoms/ cnr  is 


Q  =  No, 


(1) 


where  Q  has  the  units  cm^/cm'*.  For  a  gas  of  molecular  weight  M  and  pressure 


P  (torr) 


ON 


PaN 


Q  = 


(2) 


M 


RT 
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where  N  *  6.025  x  10  molecules/mole  and  R  =  62, 400  torr-cm/deg  K-mole. 
o 


The  product  Qt  determines  the  amount  of  scattering,  where  t  is  the  thickness 
of  gas  through  which  the  electron  beam  passes.  The  intensity  lost  from  the 
primary  beam,  i.e.,  the  beam  intensity  scattered  outside  the  objective 
aperature,  is 


I 


I  e 
o 


-(Pt)crN 

o 

RT 


(3) 


where  I  is  the  electron  beam  intensity  incident  on  the  gas  region  under 
consideration. 
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The  major  uncertainty  in  Equation  (3)  is  the  proper  value  to  use  for  Q 
(or  a).  Experimental  values  for  Q  or  a  are  known  for  only  a  few  elements 
and  compounds  and  do  not  agree  well  with  theoretical  values.  The  values 
for  SiH^  and  O2  used  for  these  calculations  were  obtained  as  follows. 

Measured  a  values*  (for  65K.V  electrons)  for  thin  solid  films  of  Be,  C, 

SiO,  Cr,  Ge,  Pd,  Pt  and  U  were  plotted  versus  the  number  of  electrons  per 
molecule  (Z)  on  a  log-log  plot.  A  linear  relationship  was  obtained  by 
least-squares  curve  fitting,  as  follows: 

1  Q 

0  =  0.14  x  10  Z.  (4) 

Equation  (4)  can  be  used  to  estimate  0  for  SiH^,  Oj,  and  other  gases  of 
interest  if  the  somewhat  doubtful  assumption  is  made  that  the  electron 
scattering  depends  linearly  on  the  number  of  electrons  present  in  the 
molecule,  independent  of  the  number  of  atomic  nuclei  present.  This  assumption 
is  valid  to  a  first  approximation  only,  for  it  is  unlikely  that  the  scattering 
for  Si(Z  =  14),  for  example,  is  7/3  times  that  of  C(Z=6)  or  14  times  that 
of  ^2=1),  Nevertheless ,  the  assumption  does  permit  determination  of  an 
interrelationship  between  the  two  parameters  1  and  t. 

The  use  of  Equation  (4)  yields  a  =  2.5  x  10_l8cm2  for  SiH,(Z=18)  and 
o  =  2.2  x  10  cm  for  02(Z=16).  These  are  the  best  available  estimates 
of  a  for  these  materials  but  may  still  be  in  error. 

The  intensity  loss  which  can  be  tolerated  in  practice  before  a  noticeable 
resolution  loss  occurs  in  the  electron  microscope  is  not  well  known.  At 
higher  magnifications  the  limit  to  proper  focusing  is  frequently  provided 
by  the  low  intensity  of  the  electron  beam,  so  scattering  losses  become 
more  severe  in  t ->at  range  of  operation.  An  intensity  loss  of  50%  is 
clearly  unacceptable,  and  a  loss  of  10%  is  more  likely  a  reasonable  upper 
limit.  An  intensity  difference  of  1%  is  not  noticeable  in  electron  micro¬ 
graphs,  so  a  loss  of  this  magnitude  is  unimportant.  Thus  the  maximum 
tolerable  intensity  loss  due  to  gas-phase  scattering  appears  to  be  between 
1  and  10%. 


*C.  E,  Hall,  Introduction  to  Electron  Microscopy  (McGraw-Hill,  New  York, 
1953)  ,  p.  246. 
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PRELIMINARY  EM6  MOD  I  El  CAT  IONS  FOR  IN  SITU  CHEMICAL  VAPOR  I)1:I*0SH  ION* 

J.  L.  Kent)’ 

North  Ame: jean  Rockwell  Corporation,  Electronics  Croup,  Anaheim,  CA  92803 


An  AE1  EM6  electron  microscope  is  being  modified  for  iji  situ  chemical 
vapor  deposition.  The  objective  is  to  observe  the  nuclcation  and  growth 
kinetics  and  structure  of  silicon  deposited  by  Sill  gas  pyrolysis  on 
substrates  of  sapphire  (single  crystal  ot-Af,o3).  4 

An  Edwards  DCF 2  thcrmoelectrically  cooled  haffle  has  been  installed  using 
a  simple  adapter  which  permitted  attachment  without  drilling  or  cutting 
of  the  EM(>  frame  or  vacuum  system.  Prior  to  installation  the  measured 
contamination  rate  was  240A/min;  afterward  the  contamination  rate  was 
less  than  lOOA/min. 

An  auxiliary  pumping  system  is  mounted  to  the  left  rear  of  the  column 
and  attached  hy  a  special  stainless  steel  flexible  bellows  to  the  pumping 
port  at  the  rear  ol  the  specimen  chamber.  The  auxiliary  system  has  an 
ultimate  presure  <1  x  Ur”  torr,  but  the  measured  pressure  in  the  specimen 
chamber  was  no  better  than  1  x  10  terr,  due  primarily  to  the  high  out - 
gassing  rate  of  the  beam  deflector  coil  immediately  above  the  specimen 
chamber.  Future  modifications  planned  include  completion  of  the  differ¬ 
ential  pumping  system  by  addition  of  an  impedance  tube  between  the 
specimen  chamber  and  beam  deflector  stage  and  a  permanent  aperture  between 
the  specimen  chamber  and  objective  lens. 

The  bottom  of  the  plate  camera  has  necn  modified  to  accept  a  transmission 
phosphor  screen  of  aluminized  P-11  phosphor  material.  A  1<>  mm  motion 
picture  camera  has  been  mounted  directly  below  this  screen  so  that 
sequences  of  in  situ  growth  may  be  photographed.  Normal  operation  of  the 
microscope  is  not  affected. 

A  CVD  microchamber.  Figure  1,  has  been  constructed  which  sits  in  the 
objective  lens  at  the  normal  4  mm  focal  length  position.  A  new  center 
post,  brazed  into  the  extensively  modified  base  of  a  standard  Al 1  speci¬ 
men  holder,  acts  as  one  electrode.  The  center  post  contains  a  lOOum  dia. 
hole  which  serves  as  a  gas  limiting  aperture.  An  outer  cylindrical 
electrode  is  fastened  by  screws  to  a  fired  lava  insulator,  which  is 
similarly  fastened  to  the  base.  A  gas-tight  seal  is  obtained  by  lapping 
the  insu'ator,  base,  and  ouier  electrode  base  flat.  Once  assembled,  the 
outer  electrode  need  not  be  disassembled  to  load  the  specimen.  A 
removable  disc  containing  a  100pm  dia.  aperture  is  held  tightly  in  place 
by  a  threaded  cap.  The  specimen  is  bold  between  two  resist ively  heated 
200-mesh  stainless  steel  grids  attached  to  the  electrodes  by  two  *80  NM 
screws.  The  holder  is  inverted  for  loading  permitting  the  specimen  and 
removable  aperture  to  be  aligned  with  the  fixed  aperture,  ('.as  is  admitted 
to  the  microchambcr  through  a  flexible  nickel  bellows,  flows  through  the 
annular  space  between  the  electrodes,  and  exhausts  through  the  apertures. 
The  space  between  apertures  is  2  mm,  which  at  ,a  Sill  pressure  of  2.9  torr 
corresponds  to  a  mass  thickness  of  0.55  ug/cm‘,  giving  a  10”.  intensity 
loss  assuming  o  =  2.5  x  10'1  cm  . 


•This  research  sponsored  by  ARPA  under  Order  1585,  monitored  by  IISAM1C0M, 
Redstone  Arsenal,  AL,  under  Contract  DAAII01 -70-C-131 1 . 
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Figure  1.  In  Situ  Chemical  Vapor  Deposition  Microchamber.  A,  center 
electrode  and  gas  limiting  aperture;  B,  lava  insulator; 

C,  outer  electrode;  D,  heating  grids;  E,  sample;  F,  removable 
aperture;  G,  aperture  cap. 
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